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I. Introduction 


The goals of this project have been to provide physical insight into the 
local structure of Io's atmosphere, the manner by which gases escape the gra- 
vitational well of the satellite and produce neutral clouds in the circum- 
planetary space, the nature of east -west and possibly magnetic longitude asym- 
metries in the plasma torus, and the stability or variability of both the 
neutral gas clouds and the plasma torus over a several year time period. 

These goals have been pursued with significant accomplishment achieved in this 
AER/JPL collaborative effort by studying and modeling the spatial morphology, 
intensity, and spacetime variability of the Io sodium cloud as preserved in 
the JPL Table Mountain Io sodium cloud data set, the most complete of the Io 
sodium cloud data sets currently available. This data set documents the 2-D 
spatial morphology of the D-line emission intensities of the cloud on the sky 
plane as a function of both Io geocentric phase angle and the magnetic longi- 
tude of Io in the System III coordinate frame over a time period from 1976 to 
1981. Highly developed and unique models for the Io sodium cloud at AER have 
provided the key tool for extracting physical insights from this data set. 

This method of approach has been proven here and shows great promise for pro- 
viding even more substantial advances in our understanding through continuing 
studies of this massive data set and other earlier but less -comprehensive data 
sets. 

The complete JPL Io sodium cloud data set is composed of 1974-1979 spec- 
tral data, 1981 spectral data, 1976-1979 Region B/C image data, 1981 Region 
B/C image data, and 1981 Region A image data. The complete data set was re- 
viewed and assessed in the first project year (see the 1985 Annual Report). 
From this review, both the quality and quantity of data acquired in 1981 were 
in general much superior to the earlier observations. Because of this fact, 
the primary emphasis in the second year (see the 1986 Annual Report) was di- 
rected toward the 1981 data set and, in particular, to the 1981 Region B/C 
images. Second year efforts included three main areas: (1) data quality re- 

view for the 1981 Region B/C images, (2) selection and data processing of a 
subset of 1981 Region B/C images, and (3) physical enhancement of the Io so- 
dium cloud model in preparation for third-year analysis of the data. A major 
and time-consuming step in the second area was the successful development and 
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application of software at the Multimission Image Processing Laboratory (MIPL) 
of JPL to remove distortions present in the brightness distribution of the 
Region B/C images . 

Research efforts in the past project year are divided into three main 
areas: (1) completion of data processing and calibration for a subset of 34 

representative images of the 1981 Region B/C cloud, (2) identification and 
examination of the basic features and observed changes in the morphological 
characteristics of the sodium cloud, and (3) physical interpretation of some 
of these observed characteristics using the highly- developed sodium cloud 
model. The modeling analysis has led to definite conclusions regarding the 
flux distribution and density profile of sodium in the exosphere of Io, the 
level of the collisional exobase in the atmosphere, the nature of the inter- 
action of the corotating plasma torus and the satellite atmosphere/corona, and 
the presence of both an east-west electric field and an inherent System III 
longitudinal asymmetry in the planetary magnetosphere. This report, which is 
organized into the three main areas noted above, summarizes only these third- 
year accomplishments. 



II. Data Processing and Calibration of 1981 Region B/C Images 


Observational parameters for the 263 Region B/C images of the 1981 JPL 
data set acquired in both the D 2 and emission lines over 14 nights between 
March 25 and June 14 are summarized in numerical form in Table 1 and in 
graphical form in Figure 1. Images each night were recorded with an integra- 
tion time of approximately 10 minutes. Three consecutive images were usually 
added to improve signal to noise. These added images provide a time-averaged 
picture of the sodium cloud over an Io phase angle interval of about 4° and an 
Io System III magnetic longitudinal interval of about 14° and, therefore, rea- 
sonably well approximating instantaneous snapshots of the Io sodium cloud. 
Prior to this project most of these images had been processed only preliminary 
to remove background signals (see 1985 and 1986 Annual reports for further 
discussion) . 

Figure 1 also shows by dots the mid-point locations for 34 representative 
images which have undergone complete data processing and calibration pro- 
cedures as part of this project at the Multi-mission Image Processing Labora- 
tory (MIPL) of JPL. The software development and application of these pro- 
cedures has required a large amount of time and project resources ,as discussed 
in the 1986 Annual Report. The observational start, end, and mid-point, con- 
ditions as well as the Image ID Number for these thirty- four representative 
images are given in Table 2. Additional information for Io's location in the 
D 2 image, the spatial scale of the image, and the relative intensity normali- 
zation factor for each of these images is given in Table 3. Each of these 
images is a combination of three (occasionally only two) of the measured *10 
minute images as is reflected in the start and end times and in the number of 
frames listed in the Image ID Number of Table 2 and Table 3. The final 
absolutely-calibrated thirty-four images in both the D 2 and D-^ sodium emission 
lines are ordered by their UT observation date and presented graphically in 
Appendix I . 

Data processing procedures applied to these thirty-four images were dis- 
cussed in the 1986 Annual Report and have been optimized for the brighter D 2 
images. Absolute calibration is achieved by multiplying each processed image 
array by the relative intensity normalization factor of Table 3 and applying 
the common conversion factor of 13 DN/kR. The absolute conversion factor is 
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1981 Region B/C Images: Observing Chronology 
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Figure 1. Observing Parameters for the 1981 Io Sodium Cloud Image Data for 
Table Mountain Observatory. The solid lines show the complete 
angular coverage for the data as listed in Table 1. The dots and 
squares show the mid-point conditions for the 34 completely 
reduced images listed in Table 2. 



Table 2 

1981 Region B/C luges: Processed by M1PL 


Start Conditions End Conditions Mid-Point Conditions 

lo Magnetic lo Magnetic Io Magnetic 


Date , Phase Longitude Phase 


of 

Observations 

Time 

(UT) 

Angle 

(deg) 

of lot 
(deg) 

Time 

(UT) 

Angle 

(deg) 

6 April 

5:00 

291.1 

303.2 

5:36 

296.3 

28 April 

4:49 

88.0 

213.1 

5:24 

92.9 


7:28 

110.3 

286.9 

7:58 

114.5 


7:39 

111.9 

292.0 

8:07 

115.8 

29 April 

3:48 

283.5 

131.3 

4:19 

287.9 


4:22 

288.3 

147.0 

4:53 

292.7 

4 May 

3:25 

217.6 

216.2 

3:55 

221.8 


4:20 

225.4 

241.6 

4:53 

230.1 


4:56 

230.5 

258.2 

5:33 

235.8 


5:49 

238.1 

282.7 

6:28 

243.6 


7:19 

250.9 

324.3 

7:52 

255.6 


9:06 

266.1 

13.7 

9:50 

272.4 

5 May 

3:59 

65.9 

178.9 

4:30 

70.3 


4:32 

70.5 

194.2 

5: 15 

76.6 


6:34 

87.7 

250.9 

7:10 

92.7 


8:03 

100.2 

292.2 

8:35 

104.7 

6 May 

5:36 

283.6 

170.4 

6:10 

288.5 


6:11 

288.6 

186.6 

6:43 

293.2 

12 May 

3:28 

46.4 

153.5 

3:55 

50.2 


4:10 

52.3 

173.0 

4:40 

56.5 


5:05 

60.0 

198.5 

6:28 

71.6 


6:30 

71.9 

238.0 

7:09 

77.4 


8:28 

88.5 

292.7 

9:07 

93.9 

13 May 

3:28 

250.1 

100.3 

4:03 

255.1 


4:43 

260.8 

134.9 

5:04 

263.7 


5:42 

269.1 

162.2 

6:07 

272.7 


6:09 

273.0 

174.7 

6:47 

278.4 


6:49 

278.7 

193.2 

7:23 

283.5 


7:25 

283.8 

209.8 

7:58 

288.5 


8:00 

288.7 

226.0 

8:35 

293.7 


8:25 

292.3 

237.6 

9:00 

297.2 

6 June 

4:19 

100.2 

292.4 

4:53 

105.0 

14 June 

3:47 

284.0 

212.6 

4:28 

289.8 


4:17 

288.2 

226.5 

5:05 

295.0 


t System III (1965) 

* images processed are redundant 


Longitude 
of lot 
(deg) 

Phase 

Angle 

(deg) 

Longitude 
of lot 
(deg) 

Image ID 
Number 

(Tape/f rames) 

319.8 

293.7 

311.5 

SIP 410/13-15 

229.4 

90.5 

221.3 

SIP 415/27-29 

300.8 

112.4 

293.9 

*SIP 415/43-45 

305.0 

113.9 

298.5 

*SIP 415/44-46 

145.6 

285.7 

138.5 

SIP 416/4-6 

161.4 

290.5 

154.2 

SIP 416/7-9 

230.0 

219.7 

223.1 

SIP 417/8-10 

256.8 

227.8 

249.2 

SIP 417/13-15 

275.3 

233.2 

266.8 

SIP 417/16-18 

300.7 

240.9 

291.7 

SIP 417/20-22 

339.5 

253.3 

331.9 

SIP 417/27-29 

34.1 

269.3 

23.9 

SIP 417/36-39 

193.3 

68.1 

186.1 

SIP 418/10-12 

214.2 

73.6 

204.2 

SIP 418/13-16 

267.6 

90.2 

259.3 

SIP 418/24-26 

307.0 

102.5 

299.6 

SIP 418/31-33 

186.1 

286.1 

178.3 

SIP 419/19-21 

201.4 

290.9 

194.0 

SIP 419/22-24 

166.0 

48.3 

159.8 

SIP 420/6-8 

186.9 

54.4 

180.0 

SIP 420/10-12 

237.0 

65.8 

217.8 

SIP 420/15-20 

256.1 

74.7 

247.1 

SIP 420/21-23 

310.8 

91.2 

301.8 

SIP 420/30-32 

116.4 

252.6 

108.4 

SIP 421/21-23 

144.6 

262.3 

139.8 

SIP 421/27-28 

173.8 

270.9 

168.0 

SIP 421/32-33 

192.3 

275.7 

183.5 

SIP 421/34-36 

208.9 

281.1 

201.1 

SIP 421/37-39 

225.1 

286.2 

217.5 

SIP 421/40-42 

242.2 

291.2 

234.1 

SIP 421/43-45 

253.8 

294.8 

245.7 

SIP 421/45-47 

308.1 

102.6 

300.3 

SIP 424/10-12 

231.5 

286.9 

222.1 

SIP 425/6-9 

248.7 

291.6 

237.6 

SIP 425/9-12 
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Table 3 

1981 Region B/C Images: Pixel and Calibration Information 


Date 

of 

Observation 

Start 

Time 

(UT) 

End 

Time 

(UT) 

Image 

Number 

Tape/frames 

Io 1 8 Location in 
the D 2 Image 
(Line) (Sample) 

Spatial 

Scale 

(km/pixel) 

Intensity 

Normalization 

Factor 

6 April 

5:00 

5:36 

SIP 410/13-15 

113 

150 

2572 

0.48t 

28 April 

4:49 

5:24 

SIP 415/27-29 

128 

155 

2680 

0.80 


7:28 

7:58 

*SIP 415/43-45 

130 

154 

2680 

0.98 


7:39 

8:07 

*SIP 415/44-46 

131 

154 

2680 

1. 59t 

29 April 

3:48 

4:19 

SIP 416/4-6 

128 

146 

2680 

0.93 


4:22 

4:53 

SIP 416/7-9 

128 

145 

2680 

0.87 

4 May 

3:25 

3:55 

SIP 417/8-10 

140 

145 

2707 

1.04 


4:20 

4:53 

SIP 417/13-15 

141 

145 

2707 

0.85 


4:56 

5:33 

SIP 417/16-18 

141 

145 

2707 

0.78 


5:49 

6:28 

SIP 417/20-22 

142 

144 

2707 

0.86 


7:19 

7:52 

SIP 417/27-29 

144 

143 

2707 

0.92 


9:06 

9:50 

SIP 417/36-39 

145 

143 

2707 

1.04 

5 May 

3:59 

4:30 

SIP 418/10-12 

143 

152 

2707 

0.97 


4:32 

5:15 

SIP 418/13-16 

143 

152 

2707 

0.94 


6:34 

7:10 

SIP 418/24-26 

145 

152 

2707 

0.80 


8:03 

8:35 

SIP 418/31-33 

146 

152 

2707 

1.00 

6 May 

5:36 

6:10 

SIP 419/19-21 

143 

143 

2734 

1. 24t 


6:11 

6:43 

SIP 419/22-24 

143 

143 

2734 

0.93 

12 May 

3:28 

3:55 

SIP 420/6-8 

144 

152 

2761 

1.21 


4:10 

4:40 

SIP 420/10-12 

146 

153 

2761 

0.99 


5:05 

6:28 

SIP 420/15-20 

147 

153 

2761 

0.87 


6:30 

7:09 

SIP 420/21-23 

149 

153 

2761 

0.78 


8:28 

9:07 

SIP 420/30-32 

150 

153 

2761 

0.87 

13 May 

3:28 

4:03 

SIP 421/21-23 

149 

145 

2761 

0.86 


4:43 

5:04 

SIP 421/27-28 

150 

144 

2761 

0.89 


5:42 

6:07 

SIP 421/32-33 

151 

144 

2761 

0.81 


6:09 

6:47 

SIP 421/34-36 

151 

143 

2761 

0.80 


6:49 

7:23 

SIP 421/37-39 

151 

143 

2761 

0.89 


7:25 

7:58 

SIP 421/40-42 

151 

143 

2761 

0.91 


8:00 

8:35 

SIP 421/43-45 

152 

144 

2761 

0.92 


8:25 

9:00 

SIP 421/45-47 

152 

144 

2761 

1.00 

6 June 

4:19 

4:53 

SIP 424/10-12 

165 

159 

2951 

0.97 

14 June 

3:47 

4:28 

SIP 425/6-9 

159 

149 

3032 

1.18 


4:17 

5:05 

SIP 425/9-12 

159 

149 

3032 

0.98 


*images processed are redundant 

tvalue in error; new value not available at this time 
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estimated to be accurate to within -±10%. The image number SIP 418/31-33 was 
chosen as the standard upon which the absolute conversion factor was estab- 
lished as indicated in Table 3 by the intensity normalization factor of 
1.00. The intensity normalization factor contains information to properly 
correct each image to the standard image for different exposure times, dif- 
ferent air masses, different earth- Jupiter distances, and differences in the 
response of the instrument over the observational period. 

In addition to these representative images, these are a large number of 
partially processed images as noted above that can be used to qualitatively 
verify and extend the characteristics of the cloud over a larger range of ob- 
servational parameter space. The representative images, however, provide the 
best images for quantitative modeling analysis purposes, for comparing images 
acquired with similar observational parameters but at different times, and for 
studying the character of the cloud with a minimum amount of instrumental and 
background distortions. 
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III. Examination of the Observational Data 


3 . 1 Correlation Studies 

In the 1986 Annual Report, five correlation studies for the 1981 Region 
B/C image data were defined and their potential value evaluated. These five 
correlation studies are defined in Table 4 and provide a grouping of these 
images based upon different classifications of their Io geocentric phase 

angle, <f > , and their Io System III magnetic longitude angle The five 

studies are important (1) in assessing the stability of the sodium cloud over 

the two and one-half month observational time interval, and (2) in identifying 

and studying the effects of solar radiation pressure, the atom source charac- 
teristics, and the plasma torus sink characteristics on the morphology of the 
sodium cloud. These studies together with the data quality reviews undertaken 
in the first and second project years were used in selecting the thirty- four 
representative images that were processed and calibrated by MIPL. 

The five correlation studies of Table 4 have been applied in the third 
project year to these 34 fully-reduced and representative 1981 Region B/C 
images. The representative images, property ordered and compared for each of 
the five studies, are presented, respectively, in the five appendices I, II, 
III, IV, and V. The results of these studies are discussed below and form the 
basis for a number of modeling studies, some of which have been undertaken in 
this project and are discussed in section IV. 

3.2 Stability of the Sodium Cloud 

It is important in examining and comparing the 1981 Region B/C images to 
determine first if the sodium cloud was relatively stable over the -12 week 
observing window for observations summarized in Table 1 or more importantly 
over the -10 week interval for the selected observations in Appendix I. Since 
the cloud morphology is time variable and depends upon both the Io geocentric 
phase angle (-42.5 hour period) and the Io System III longitude angle (-13 
hour period), the most appropriate comparison is for images that have a recur- 
rence of these two angle pairs (-32 day period). Fortunately, there are three 
such occurences in the 1981 Region B/C data set (see Figure 1). Image pairs 
from two of these occurrences (one east and two west of Jupiter) are compared 
(study 2 of Table 4) in Appendix II. Comparisons made from the east image 
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Table 4 


1981 Region B/C Images: Five Correlation Studies^ 


Consecutive Image Study: 

H t). *(t) 

Cloud Stability Study: 

~ (^2> ^2^ 

image sequences 

1 separated by 32 days 

East/West Image Study: 


*E * % 

only one image -sequence pair 

« * E + 180° 

no image pairs 

*E " *W 

many comparison images 


IV. System III Variability Study: ~ 4> 2 < *1 ** #2 

V. Io Phase Angle Variability Study: * $2’ *1 “ ^2 


— Io Geocentric Phase Angle 
® - Io System III Magnetic Longitude Angle 
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pair (image set 2. A of Appendix II) show that the sodium cloud (apart from 
some recognized scatter in the data because of signal to noise effects) on 
May 5 and June 6 are very similar both in absolute intensity and in morpho- 
logical structure. The same is also true for the two west image pairs (image 
sets 2.B and 2.C of Appendix II) of May 13 and June 14 for which a System III 
dependent "bite -out" of the cloud south of Io (to be discussed later) is re- 
produced in all images. 

In order to test more extensively the stability of the sodium cloud for 
data in Appendix I, the next best comparison of images beyond the one dis- 
cussed above is that for similar Io geocentric phase angles but different Io 
System III longitude angles (study 4 of Table 4) . This comparison is pre- 
sented in Appendix IV. To address the question of stability in these images, 
it is necessary to compare as many images as possbile for different observing 
dates in order to identify and eliminate System III changes .from considera- 
tion. Comparison of three image pairs (image sets 4. A, 4.B and 4.C) on 
March 5 and March 12 all show similar structure, spatial extent and bright- 
ness, although a System III correlated asymmetry in the sodium distribution 
north and south of Io is suggested. This will be discussed further in 
Section 3.3. An east image (image set 4.C) acquired on April 28 .has a similar 
spatial structure and brightness, but the spatial extent of the forward cloud 
and the trailing north directional feature are not as pronounced. This does 
not seem to be due to any observing/reduction problem since it is also pre- 
sented in other April 28 images (SIP 415/43-45, SIP 415/44-46) at larger east 
phase angles. 

For the sodium cloud west of Jupiter, more extensive image comparisons 
for similar Io phase angles are possible and are presented in Appendix IV. 
Inspection of image sets 4.E and 4.F shows similar structure, extent, and 
brightness of the forward cloud when the expected changing nature of the di- 
rectional feature in the trailing cloud is considered. The image set 4.G 
extends the earlier discussed image pair 2.B (May 13, June 14) by adding 
April 29 and May 6 images. Similarly, image set 4.H extends the earlier dis- 
cussed image pair 2.C (May 13, June 14) by adding April 29, May 6, another 
May 13, and April 6 images. In addition to a clearly delineated System III 
correlated south "bite-out", the images compare very well in both image 
sets. For image set 4.G, the forward cloud of the April 29 image is not as 
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extended and the sodium not as bright nearer Io. For image set 4.H, the for- 
ward cloud of the April 29 image is of normal length but the cloud is still 
not as bright near Io, while the April 6 image is generally weaker and may be 
attributed to data of lesser quality. 

In summary, the spatial structure of the sodium cloud appears quite 
stable over the -10 week period of data in Appendix I. There is only a hint 
of weaker sodium brightness in the cloud in April images. The basic structure 
of the cloud remains the same with the identification of stable features and 
several time -dependent periodic signatures readily apparent. These features 
and time -dependent signatures may therefore be safely classified and further 
studied by additional comparisons of images in Appendix I. The main features 
and time -dependent signatures of the sodium cloud are discussed below. 

3.3 Features and Periodic Signatures of the Sodium Cloud 

The main features and periodic signatures of the sodium cloud are listed 
in Table 5. All of these features and signatures have been either discovered 
or confirmed in the JPL Table Mountain Io sodium cloud data set. The features 
and signatures of Table 5 as well as some less stable signatures are discussed 
separately below. 

East-West Intensity Asymmetry 

The east-west intensity asymmetry near Io was discovered by Bergstralh et 
al. (1975, 1977) from the earliest single slit spectrographic observations ac- 
quired for the JPL Table Mountain Io sodium cloud data set in 1974 and 1975. 
These observations showed that the sodium cloud in a 3 x 8 arcsecond slit cen- 
tered on Io was (on the average) brighter by about 25% when the satellite was 

east of Jupiter than when it was west of Jupiter. Similar slit spectrographic 
observations were acquired for the JPL data set in 1976-1981 and were reviewed 
in the 1985 Annual Report. These observations in addition to 1981 Region A 
image data are most valuable in extending the results of Bergstralh et al. 

The spectra and image data have, however, not yet been properly reduced for 
this purpose. 

The spatial region near Io observed by Bergstralh et al. is essentially 
confined to sodium in the Lagrange sphere (~ 10^ km in radius) of the satel- 
lite. This region is not observed in the 1981 Region B/C images of Appendix I 
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Table 5 


Features and Time-Varying Periodic Signatures of the Sodium Cloud 

1. East-West Intenstity Asymmetry near Io 

2. Predominance of a Bright Forward Cloud 

3. East-West Orbital Asymmetry of the Forward Cloud 

4. South Bite-Out in the West Cloud 

5. Directional Feature (Trailing Cloud) 
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because of the larger circular occulting mask (4.5 arcsec in radius and cor- 
responding to a distance varying between 1.46 x 10^ and 1.71 x 10^ km) cen- 
tered on Io. A small increase of order 25% in the cloud brightness east of 
Jupiter is, however, not readily distinguishable in the 1981 Region B/C cloud 
at distances from Io of the 1 to 2 kR contour levels when visual comparisons 
are made between the east-west image sets in Appendix III and when proper ac- 
count is taken for different Io phase angle and Io System III longitude 
effects. 

To test more quantitatively for a brightness enhancement at distances 
from Io larger than the occulting mask (i.e. >17,1000 km radius) but smaller 
than the 1 to 2 kR contour levels, the average D 2 brightness of the cloud in 
an annulus between 19,000 km and 21,000 km was computed for 31 of the 34 
images in Appendix I for which the intensity calibration is established (see 
Table 3). The results are shown in Figure 2 as a function of the Io geocent- 
ric phase angle. From Figure 2, it is immediately apparent that the east 
cloud is generally brighter just beyond the occulting mask, although there is 
scatter in the data. For the 12 east images in Figure 2, the average annulus 
brightness is 3.15 kR while for the 19 west images, the annulus brightness is 
2.24. The east cloud annulus is therefore, on average, -40% brighter than the 
west cloud annulus. This confirms the east-west asymmetry discovered by 
Bergstralh et al. (1975, 1977). It also indicates that the east-west bright- 
ness ratio at a distance from Io greater than that sampled by Bergstralh et 
al. is larger than their average 1.25 value and argues that the intensity 
asymmetric mechanism is primarily sink related in contrast to source related. 

Predominant Bright Forward Cloud 

The existence of a brighter sodium cloud inside of Io's orbit and ahead 
of the satellite's location can be traced back to some of the earliest slit 
observations of the D-line emissions acquired in 1973 and 1974 (i.e. pre- and 
post-Pioneer 10) by Macy and Trafton (1975a, b) and Munch and Bergstralh 
(1977). Later two-dimensional image data acquired for the sodium cloud in 
1976 and 1977 (Murcray 1978; Matson et al. 1978; Murcray and Goody 1978) first 
revealed the presence of a bright and predominantly forward sodium cloud. The 
large number of images in the 1981 portion of the JPL Table Mountain data set 
provide the first opportunity to refine and examine the spatial character of 
the forward cloud more quantitatively as discussed in the next subsection. 
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Figure 2. Sodium Cloud D 2 Brightness in an Annulus about lo. The average brightness in an annulus 

centered on lo with inner radius of 19,000 km and outer radius of 21,000 km is shown for thirty- 
one of the fully-reduced 1981 Region B/C images in Appendix !• 



East-West Orbital Asymmetry of the Forward Cloud 

Examination of the fully-reduced 1981 images in Appendix I shows that the 
forward cloud is always present. The appearance of the forward cloud on the 
sky plane, however, changes as Io moves about Jupiter. These changes are con- 
sistent with changes recognized earlier (Smyth and McElroy 1978; Matson et al. 
1978) in 1976 and 1977 images acquired for the JPL data set (Matson et al. 
1978) and acquired independently at Harvard University (Murcray 1978; Murcray 
and Goody 1978). The much larger number of 1981 images in the JPL data set, 
however, allow one to track these changes in more detail as a function of the 
satellite phase angle. In particular, the forward cloud on the sky plane can 
be seen to be approximately parallel to the line of sight at a critical east 
phase angle of -65-70° and a critical west phase angle of -235° from the 
images of Appendix I. This was first discovered by Goldberg et al. (1978) and 
later quantified for the 1981 JPL image data by Goldberg, Garneau, and LaVoie 
(1984). The fact that the east and west critical phase angles do not differ 
by 180° and furthermore that the forward cloud west of Jupiter is generally 
significantly longer than the forward cloud east of Jupiter (see images in 
Appendices I, III and IV) directly indicate an east-west orbital asymmetry in 
the cloud. 

South Bite-Out in the West Cloud 

A zone of deficient sodium southwest of Io was noted in the 13 May 1981 
Region B/C images in the brief review of the JPL Table Mountain data set con- 
ducted by Goldberg, Garneau, and LaVoie (1984). This "bite-out" of a south 
portion of the west cloud can be seen to develop in time (see images in 
Appendix I) as a function of the Io System III longitude angle. The 
n bite-out" is absent for an Io System III longitude angle of 168° (image 
SIP 421/32-33), is noticeable at an angle of 183.5° (image SIP 421/34-36), is 
increasing as the angle increases beyond 200° , and is still present near the 
end of the observational sequence where the angle is 245.7° (image 
SIP 421/45-47). 

Examination of additional images in the 1981 Region B/C data set shows 
that this south "bite-out" of the west cloud is a normal and periodic feature 
of the Io sodium cloud for at least a seven week interval. From Figure 1, it 
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is immediately evident that images on April 29, May 6, and June 14 have 
similar values for their Io System III longitude and Io geocentric phase 
angle-pairs. Comparisons of the fully-reduced images on these dates in Appen- 
dix I confirm this behavior. Examination of raw data for images on April 29 
with Io System III longitude angles larger than 154.2° (image SIP 416/7-9 in 
Appendix I) also reveals the same behavior. 

The south "bite-out” of the west cloud is also a relatively stable fea- 
ture of the sodium cloud over a several year time frame. The south "bite-out" 
is consistent with two 1977 images, and one 1979 (Voyager 1 encounter) image 
of the JPL Table Mountain data set published earlier (Goldberg et al. 1980) 
even though the integration times for these sets are -1-2 hours (see the 1985 
Annual Report). The south "bite-out" of the west cloud is also present in the 
very high quality (but unpublished) sodium cloud data acquired on 13 June 1983 
and 16 August 1984 by Morgan (1984). The physical meaning of the "bite-out" 
will be explored and explained in section IV. 

Evidence for South Deficiency in the East Cloud 

It is natural to ask if there is a sodium deficiency south of Io in the 
east cloud for the same satellite System III longitude range for which the 
"bite-out" occurs in the west cloud. From Figure 1, observations obtained on 
April 28, May 5, and May 12 (seven-day intervals and one day earlier than the 
west cloud images on April 29, May 6, and May 13 discussed above) are very 
well suited for such an investigation. 

The three fully-processed east images for April 28 in Appendix I have a 
sodium cloud that is brighter north of Io. The image (SIP 415/27-29) at a 
satellite System III angle of 221.3° is much more asymmetric north to south 
than the two similar images (SIP 415/43-45 and SIP 415/44-46) at larger angles 
near 300° . 

The images in Appendix I for May 5 show a fairly symmetric cloud north 
and south of Io for a satellite System III angle of 204.2° (image 
SIP 418/13-16) or less, a cloud enhanced north of the satellite for an angle 
of 259.2° (image SIP 418/24-26), and a cloud symmetric north and south of the 
satellite for an angle of 299.6° (image SIP 418/31-33). Examination of a 
number of additional images on May 5 that have not been fully reduced rein- 
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forces this behavior and indicates that the cloud is symmetric for an Io 
System III longitude angle of 210° , that it becomes brighter north of Io for 
larger angles, and that the asymmetry then decreases and is absent by an angle 
of 294* . The spatial morphology of the south deficient zone of sodium in the 
east cloud is, however, characterized by a steeper gradient in the brightness 
for the whole sodium cloud south of Io rather than a sharp and inclined bound- 
ary for the deficient zone as is present in the south bite-out of the west 
cloud. Finally, the image of June 6 (see Appendix II) exhibits only a slight 
north brightness enhancement when compared with the May 5 image (SIP 418/31-3) 
at near identical Io phase and Io System III longitude angles . 

The images on May 12 do not, however, show the south zone of deficient 
sodium that is indicated in the May 5 images and suggested by the April 28 
images. Examination of May 12 images in Appendix I and additional images that 
have not been fully reduced shows that the cloud brightness is essentially 
symmetric north and south of Io for satellite System III longitude angles 
ranging from -210° to 302°. 

The different north- south behavior of the east cloud images on May 5 and 
May 12 for the same range of Io System III longitude angles (-180° -300°) is 
surprising in the light of the consistency and repeatability of the south 
"bite-out" signature that occur in the west cloud on May 6 and May 13 over 
this same range of Io System III longitude angles. The long-term stability of 
the "bite -out" in the west cloud and the temporal or more erratic behavior of 
the deficient sodium south of Io in the east cloud raise some very interesting 
questions for Io and the magnetosphere that are discussed in Section IV. 

Directional Features 

The presence in the trailing cloud of an elongated feature having a 
changing north/south inclination relative to the satellite plane was disco- 
vered in 1980 (January 5 - April 22) and 1981 (March 10 - April 13) image data 
of Pilcher et al. (1984). The time varying inclination of the directional 
features was shown to be periodic and correlated with the Io System III longi- 
tude angle. In the data of Pilcher et al., the directional feature was in- 
clined north of the satellite orbit plane for Io System III (1965) longitude 
angles between -170° and 360° , and south of the orbit plane for longitudes 
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outside of this range. The south inclination of the feature was observed to 
decrease to zero as the Io System III longitude angle approached -150-160°, 
and the north inclination of the directional feature was observed to increase 
for larger angles. 

To investigate the behavior of the directional features in the JPL data 
set, the 1981 Region B/C images of Appendix I have been arranged in Appendix V 
(Study 5 of Table 4) in ascending order of the Io System III (1965) longitude 
angle for both the east cloud and west cloud separately. Examination of the 
images in Appendix V shows that the inclination of the directional feature is 
indeed the same as that discovered in the 1980 and 1981 image data of Pilcher 
et al. (1984). This north/south correlation of the inclination of the direc- 
tional feature with the Io System III (1965) longitude angle for images in 
Appendix V is shown in Figure 3. 

As noted in the 1985 Annual Report, the correlation in Figure 3 is also 
reflected (although not at that time discovered) in 1976-1979 sodium cloud 
images of the JPL data set published by Goldberg et al. (1980). In particu- 
lar, comparison of an image on 5 March 1979 (Voyager 1 encounter day) and an 
image on 5 December 1977 having the same average Io geocentric phase angles 
(256°) showed that while their forward clouds were similar in length and 
brightness, their trailing clouds had different brightness distributions north 
and south of the satellite orbit plane. The trailing cloud of the March 5 
image (centered at an average Io System III angle of 34.8°) was distributed 
south of the satellite orbit plane, while the trailing cloud of the December 5 
image (centered at an average Io System III angle of 142°) was distributed 
symmetrically about the orbit plane. This is the expected behavior based on 
Figure 3 and implies that the directional feature is a permanent feature of 
the sodium cloud. 
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IV. Interpretation of the Observational Data 


4 . 1 Mode 1 Deve lopment 

Overview 

In order to analyze effectively the main features and time -varying pe- 
riodic signatures of the sodium cloud in Table 5, it is necessary to have a 
model that contains to a reasonable level of accuracy a description of the dy- 
namics , sink processes, and the satellite source mechanism for sodium atoms in 
the Io- planetary -magnetosphere system. The orbital dynamics of sodium atoms 
in the cloud model are determined by the gravitational fields of both Io and 
Jupiter and by the acceleration of solar radiation pressure produced as atoms 
undergo solar resonance scattering in the D*^ and D 2 emission lines. All three 
factors were already implemented in the AER sodium cloud model at the be- 
ginning of this project. The sink for sodium atoms is determined by their in- 
teraction with the plasma torus which is dominated by electron impact ioniza- 
tion. Just prior to the beginning of this project, the spacetime dependent 
sodium sink produced by the effective oscillation of the plasma torus about 
the satellite orbit plane was explicitly included in the model. The plasma 
torus was assumed in this description to be longitudinally symmetric, and its 
oscillatory motion was based upon a tilled and planet- centered magnetic dipole 
field with a corotating heavy-ion plasma. A significant number of model im- 
provements for the plasma torus sink have been implemented in this project and 
are discussed below. The satellite source for the sodium cloud has a most 
probable speed for the atom flux distribution that has been estimated to be 
—3 km sec but the shape of initial velocity dispersion of the flux dis- 
tribution has remained somewhat uncertain. Particular emphasis in this pro- 
ject has been directed toward more carefully determination of the underlying 
physical nature of this source and its velocity dispersion as discussed below. 

Lifetime of Sodium in the Plasma Torus 

The description of the plasma torus in the sodium cloud model to deter- 
mine the spacetime varying sodium lifetime has been improved. In the first 
project year, improvements were made in the electron temperature and density 
by incorporating new results acquired from the analysis of both the Voyager 1 
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PLS electron data and ion data. These improvements were discussed in the 1985 
Annual Report. In the second project year, the lifetime description of sodium 
in the planetary magnetosphere was improved by including the offset dipole na- 
ture of the magnetic field as well as the effects of a constant east-west 
electric field (Barbosa and Kivelson 1983; Ip and Goertz 1983; Goertz and Ip 
1984). These improvements were discussed in the 1986 Annual Report. 

A model documentation paper describing a general mathematical formalism 
for the neutral clouds of Io and including a number of plasma-neutral cloud 
interactions in the torus was completed in the second year and is included in 
Appendix VI. This paper (Smyth and Combi 1988a) is scheduled to be published 
in the 1988 April issue of The Astrophvsical Jou rnal Supplement . The general 
mathematical formalism was applied to the sodium cloud in a second paper 
(Smyth and Combi 1988b) that was completed in the third project year and that 
is scheduled to be published in the 1988 May 15 issue of The AstronhvsiraT 
J ournal . This second paper is attached in Appendix VII. The implementation 
of these improvements in the lifetime description of sodium was driven by the 
presence of physically based signatures observed in the cloud data and the 
plasma torus. A paper (Smyth and Combi 1987) correlating east-west asymmet- 
ries in the Jovian magnetosphere and the Io sodium cloud that illustrates some 
of these physical signatures was published in Geophysical Research Letters and 
is attached in Appendix VIII. 

The lifetime description for sodium was further refined in the third 
project year by including an inherent System III asymmetry in the plasma torus 
in the longitude range of -180* to 300* . This addition was required by the 
presence of the south "bite-out" signature in the west sodium cloud and will 
be discussed in Section 4.2. A physically corresponding but anticorrelated 
System III brightness asymmetry in the plasma torus S + (6716A° and 6731A*) ion 
emission has been observed from ground-based telescope since 1976 and was 
later discovered in the S'*- 1 - (9532A) ion emission. Pertinent S + spectral data 
for this asymmetry have been acquired from ground-based observations in 1976 
(Trafton 1980; Trauger, Munch, and Roesler 1980), in 1977 and 1978 (Trafton 
1980; Pilcher and Morgan 1980); in 1979 (Trafton 1980; Morgan and Pilcher 
1982), and in 1981 (Morgan 1985a,b). Pertinent image data for S + and/or S ++ 
emissions have more recently been obtained from ground-based observations in 
1980 (Roesler et al. 1982), in 1981 (Trauger 1984; Pilcher and Morgan 1985; 
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Oliversen et al. 1987), in 1982 (Trauger 1984), in 1983 (Trauger 1984; Pilcher 
and Morgan 1985) and in 1984 (Trauger, private communication) . 

A comparison of these plasma observations indicates that two separate 
longitudinally asymmetric components appear to exist in the torus. The first 
component is centered at -180° -200° with an angular width of -90° or more, and 
is a permanent feature of the torus. When compared to the normal brightness 
of the plasma torus at the diametrically opposite magnetic longitude angle, 
the first component is typically 2-5 times brighter (and at extremes, 10 times 
brighter) with significant brightness changes occuring on a time scale of 
months to years. The second component is centered at -280° -290°, has a compa- 
rable angular width to the first component, and is sometimes present and some- 
times absent in the torus. The brightness of the second component may be com- 
parable to or less than the brightness of the first component. For the ob- 
serving time interval of the 1981 Region B/C images in Appendix I, the first 
and second longitudinally asymmetric components were similar in angular size 
and brightness, and both of these components had an enhanced brightness of -4 
and a plasma density enhancement of ~2 when compared to the plasma torus at 
the diametrically opposite longitude location (Morgan 1985a, b; Morgan 1987, 
private communication) . 

Initial Flux Distribution for Sodium Atoms at Io 

Significant efforts in the third project year have been expended in de- 
termining the character of the flux distribution for the upward-moving sodium 
atoms ejected from Io that is consistent with the general morphology and 
brightness of the sodium cloud. This work is included in the paper entitled 
"A General Model for Io's Neutral Gas Clouds. I. Application to the Sodium 
Cloud” in Appendix VII (see particularly Appendix D of this paper) . A 
raodified-sputtering flux distribution to describe the incomplete collisional 
cascade processes between plasma ions and the atmosphere of Io is formulated 
and shown in the paper to provide the proper type of nonthermal sodium source 
required to simulate the morphological characteristics of the cloud. 
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4.2 Modeling of Sodium Cloud Features and Time Varying Signatures 

The modeling analysis of this project has been undertaken in two steps. 
The first step was to verify the analytic usefulness of the model by success- 
fully establishing an overall consistent and physical explanation for the most 
basic morphological characteristics of the cloud. The second step was to test 
and to refine this physical picture by using the model to analyze specific 
1981 Region B/C images of the JPL data set. The main features and time- 
varying periodic signatures of the data set are summarized in Table 5. 

Basic Morphological Characteristics of the Cloud 

The most basic morphological characteristics of the cloud are (iy the 
predominance of a bright forward cloud with its east-west orbital asymmetry, 

(2) the presence of a less bright trailing cloud with its north-south alter- 
nating directional feature, and (3) an east-west intensity asymmetry near Io 
such that the D-line emissions are brighter, on average, when the satellite is 
east of Jupiter. The east-west orbital asymmetry exhibits itself in an east 
forward cloud which is less extended arid is generally only comparable in 
length to its directional feature, and a west forward cloud which is signi- 
ficantly longer than the east forward cloud and is more inclined toward 
Jupiter (relative to a tangent line to the satellite orbit at Io's location) 
by -10-15°. 

The underlying physical explanation for the nature of the forward cloud 
and the east-west intensity asymmetry of the cloud have been explored exten- 
sively in this project. The results of these modeling efforts are presented 
in papers (Smyth and Combi, 1988b and 1987) attached in Appendices VII and 
VIII. The significant conclusions and physical insights gained from these 
modeling efforts for Io's atmosphere and for the plasma torus are briefly sum- 
marized later in Section 4.3. 

The north- south alternating nature of the directional feature in the JPL 
data set was shown in Section 3.3 (see Figure 3) to be correlated with the 
System III longitude angle of Io. This correlation is identical to that 
determined earlier by Pilcher et al. (1984). For the image data of Pilcher et 
al., the directional features were analyzed with the 1983 sodium cloud model 
of AER. The Io System III longitude behavior of these directional features 
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was adequately explained by Pilcher et al. in terms of a source of high- 
velocity (~20 km sec”^) sodium and a spacetime dependent sodium sink produced 
near Io by the effective geometric and periodic oscillatory motion of the 
plasma torus. The modeling of the directional feature in the JPL data set 
need not therefore be explicitly considered further in this project. The 
velocity dispersion for the sodium flux distribution at Io necessary to pro- 
vide the directional feature, the forward cloud, and the east-west intensity 
asymmetry in a consistent manner is discussed in Appendix VII. 

South Bite-Out of the West Cloud 

The Io System III longitude evolution of the south bite-out feature of 
the west cloud is illustrated in Figure 4 by four Region B/C images of 
Appendix I obtained on May 13, 1981. This bite -out is not produced by model 
computations illustrated in Figure 5 using the best sodium lifetime descrip- 
tion (case C in Appendix VII) of Smyth and Combi (1988b) and an isotropic 
sodium source with a flux of 2 x 10® atoms sec'^ (referenced to the satellite 
surface area) ejected monoenergetically (2.6 km sec'^) from the satellite exo- 
base assumed to have a radius of 2600 km. The lifetime description is based 
upon a tilted and offset magnetic dipole field in the presence of a constant 
east-west electric field of 2.8 mV/m for a heavy- ion magnetospheric plasma 
with no inherent System III longitudinal dependence. The less extended 
trailing cloud beyond Io's orbit in the model results of Figure 5 is under- 
standable due to the absence of higher velocity components in the adopted 
monoenergetic source. The adopted monoenergetic source is chosen near the 
most probable speed of the flux distribution and is quite adequate for our 
present purposes of investigation. 

Examination of 1981 plasma torus data (Morgan 1985a, b; Pilcher and Morgan 
1985; Oliversen et al. 1987) as discussed earlier in Section 4.1 indicated an 
S and S emission enhancement in the torus near Io's orbit in the System III 
longitudinal region of the -180° -300°. The brightness enhancement may be at- 
tributed primarily to an increase in the electron density of the so-called S + 
"ribbon feature” (Trauger 1984) or "field-aligned feature" of Pilcher, Fertel 
and Morgan (1985; Pilcher and Morgan 1985; Morgan and Pilcher 1982) that is 
located inside of Io's L- shell and to a S ++ plasma extension of this feature 
beyond Io's orbit. The longitudinal asymmetry in the torus appears to have 
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Figure 4. 1981 May 13 Region B/C Image Data* Four D 2 emission images of the sodium cloud as measured on 

the sky plane from the Table Mountain Observatory are shown in contour plot format and in proper 
spatial relationship to Jupiter, Io's location, and Io's circular orbit which appears as a thin 

ellipse. Contour brightness levels in kiloRayleighs, from outside to inside, are spaced as 
follows: 0.2, 0.5, 1.0, 2.0, 5.0, 10.0, 20.0. 
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two components, the first one centered at a System III longitude angle A^j^ 
(-180-200° ) and the second one centered at a{|£ (-280-290°). The lifetime of 
sodium through electron impact ionization then has an inherent longitudinal 
dependence . 

An explicit functional dependence of the sodium lifetime r on the atom's 
System III longitude angle, Ajjj, and radial position r in the plasma torus 
coordinate frame was therefore adopted in the sodium cloud model. This new 
lifetime description is given by 



and was suggested by more recent plasma torus models of Morgan (1987, private 
communication). Here r Q — i/ Q ^ is the lifetime description for sodium deter- 
mined by case C of Smyth and Combi (1988b) . The parameters aj determine the 
relative longitudinal strength of the two asymmetric plasma torus components 
while the parameters Oj determine their relative angular width. The radial 
dependence R(r) of the asymmetry is centered at r - r Q and has an inner scale 
length <7 inner that differs from its outer scale length % uter . During the 
1981 JPL observations of the sodium cloud, best estimated values for these 
parameters are given as follows 
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Using the parameter values given by (3), the electron density is enhanced in 
the active longitudinal region by a factor of two (as in the best plasma torus 
models of Morgan 1985b) over the Voyager 1 value it has at the diametrically 
opposite longitudinal angle in the plasma torus. The sodium lifetime is cor- 
respondingly reduced by the same factor. 

New model calculations for the four observations in Figure 4 were per- 
formed using the asymmetric longitudinal sodium lifetime description given by 
the expressions (1), (2), and (3), but using the same satellite source adopted 
in Figure 5. These results are shown in Figure 6 and clearly simulate the ba- 
sic character of the Io System III longitude evolution of the south bit-out 
feature. Comparison of model results in Figure 5 and Figure 6 also shows that 
the increased sink active in Figure 6 has reduced the size of the forward 
cloud so that it is now only a little longer than in the data of Figure 4. 

The 1 kR and 2 kR level contours of the forward cloud in Figure 6 are, how- 
ever, located at a similar distance from Io as in the data in Figure 4. Re- 
ducing only the active flux region of the model exobase to a band source 
(which excludes a portion of the leading and trailing hemispheres, used in the 
paper in Appendix VII) and adopting the asymmetric longitudinal sodium life- 
time description, monoenergetic model calculations for the first and third 
images of Figure 4 are shown in Figure 7. Comparison of appropriate images in 
Figure 4 and Figure 7 shows that the forward cloud lengths at the 0.2 kR level 
are similar in length, but that the model results for the forward cloud closer 
to Io are not bright enough by approximately a factor of 2. This suggests 
that a more correct specification of the source would lie somewhat immedately 
between the band and isotropic ejection condition. This is similar to the 
conclusion reached in the paper in Appendix VII from modeling the east-west 
intensity asymmetry measured by Bergstralh et al. (1975, 1977). The appro- 
priate sodium source may then correspond to a flux of -2 x 10^ cm"^ sec’^ for 
isotropic ejection or a slightly larger flux for band ejection which both lead 
to a total source strength of -1 x 10^ atoms sec"^*. 
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West Cloud Model Calculations for an Isotropic Source an 
Torus Sink* Same description as Figure 5 with the excep 
plasma torus sink for sodium (see text) has now been ass 
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Figure 7. West Cloud Model Calculations for a Band Source and a Longitudinally Asymmetric Plasma Torus 

Sink. Model calculations for the D 2 emission brightness of two of the four images in Figure 4 
are shown with the same adopted contour levels. The plasma torus sink and sodium source are the 

same as in Figure 6 with the exception that the active exobase area has now been reduced to a 
band source* 


East Cloud Analysis 


Two model calculations for the sodium cloud at an Io geocentric phase 
angle of 90.5° and an Io System III angle of 221.3° assuming the longitudinal- 
ly asymmetric plasma torus sink given by expressions (1) , (2), and (3) and the 
two previously adopted source conditions with a flux of 2 x 10® atoms cm’^ 
sec ^ are shown in Figure 8. The Io geocentric phase and Io System III 
longitude angles for the model results are identical to the April 28 image 
SIP 415/27-29, and differ only appreciably for the latter angle from the May 5 
image SIP 418/24-26 (i.e. 90.2°, 259.3°) and the May 12 image SIP 420/30-32 
(i.e. 91.2°, 301.8°). Comparison of the two model results in Figure 8 
illustrates how the band source reduces both the length and brightness of the 
forward cloud (as also noted in the paper in Appendix VII) substantially from 
that of an isotropic source. Comparison of the band source results in 
Figure 8 with the three images noted above, however, indicates that the 
forward cloud predicted by the model, although only slightly longer, is too 
bright at the 1 to 2 kR level by at least a factor of two. This is in 
contrast to the model results for the west cloud at comparable Io System III 
angles in Figure 6 and Figure 7, where the isotropic source matches these 
observed 1-2 kR brightness levels and the band source produces a cloud that is 
too dim at these observed brightness levels by a factor of —2. 

The ability of the cloud model for the band source to reasonably match- 
the west cloud shape but underestimate its brightness, and to reasonably match 
the length of the east forward cloud but overestimate its brightness is most 
interesting. This likely indicates that the angular distribution of the so- 
dium source is oriented more in the vector direction of motion of the satel- 
lite in the east than in the west. Such an east-west asymmetry in the source 
would not appreciably alter the brightness of the cloud near Io (where it is 
constrained by the east- west intensity asymmetry observations) but would 
change the orbital trajectories of sodium atoms at larger distances from the 
satellite so that not as many atoms would be able to populate the east forward 
cloud and more atoms would be able to populate the west forward cloud. The 
comparison of the isotropic and band source in Figure 8 is a simple example 
where the angular distribution has been changed, but in this case the same re- 
duced angular distribution (i.e. the band source) has been assumed for all 
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Figure 8. East Cloud Model Calculations for a Longitudinally Asymmetric Plasma Torus Sink. Model calcula- 
tions for the D 2 emission brightness of the 1981 April 28 image (SIP 415/27-29) are shown on the 
sky plane with the same contour levels adopted in Figure 4. The plasma torus sink for sodium is 
the same as in Figure 6. For both the isotropic and band sources, monoenergetic (2.6 km sec”^) 

ejection of sodium was assumed with the same exobase radius and uniform flux as adopted in 
Figure 5. r 



east and west Io geocentric phase angles. The east-west asymmetry in the an- 
gular distribution of the source that is suggested might be produced by the 
interaction of the plasma torus ions with Io's local atmosphere if the atmo- 
sphere were enhanced above the sunlit hemisphere of the surface. The investi- 
gation of this very interesting possibility, however, requires additional sup- 
port for research beyond that available in this project. 

The model results of Figure 8 show an east cloud that exhibits little or 
no asymmetric distribution of sodium north and south of Io, in contrast to the 
south bite-out feature present in the model results for the west cloud 
(Figure 6 and Figure 7). Additional model results (not shown) for the May 5 
east image (SIP 418/24-26) and the May 12 east image (SIP 420/30-32) that 
cover the Io System III longitude angle range of the asymmetric sink substan- 
tiate the results of Figure 8. The reason there is no deficient sodium signa- 
ture south of Io in the east cloud model is because the east-west electric 
field has displaced Io radially Inward from the location of the asymmetric 
sink and thereby has literally removed the sodium from the spatial region of 
enhanced ionization. This is readily apparent by studying Figure 9 where the 
sodium lifetime (case C of Smyth and Combi 1988b) in the plasma torus at east- 
ern and western elongation is shown in relationship to Io's motion (dark oval) 
during one Io System III longitude period. The Io System III longitude angle 
on the oval in the centrifugal equator plane at the smaller radial displace- 
ment is 110° and increases as the satellite moves upward along the oval and 
reaches 200° at the top of the oval. The longitudinally asymmetric plasma 
torus sink for sodium is centered on the centrifugal plane at a radial dis- 
tance of 5.65 Rj near the middle of the oval at eastern elongation and some- 
what to the left of the oval at western elongation. The small inner radial 
scale height (a- nner - 0.125 Rj) of the inherent System III enhancement causes 
the enhancement to decay so rapidly that it is unable to reach Io's location 
in the 110° -200° sector in the east, while the outer radial scale height 
(a ou ter * 1.0 Rj) of the inherent System III enhancement allows the enhance- 
ment to be felt at almost full force at Io's location in the 110° -290° sector 
in the west. 

The north- south symmetry of the east cloud model results in Figure 8 is 
similar to the Region B/C images observed on March 12, but not the images ob- 
served on March 5 and April 28 which show some deficiency in the sodium 
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brightness south of Io as discussed earlier in Section 3.2. The apparent time 
variability or erratic behavior of the east cloud and the stability of the 
west cloud regarding a zone of deficient sodium south of Io are likely to be 
related to a difference (or time variability) in the strength of the east-west 
electric field east and west of Jupiter. It is clear from Figure 9 that the 
west cloud is in a rather stable relationship with the location of the inher- 
ent System III plasma torus enhancement for small modulations of the electric 
field strength west of Jupiter. This is not, however, true for the east cloud 
because of the very small scale height cr£ nner and the rather close radial 
proximity of Io's orbit to the System III enhancement in the plasma torus. 
Small changes in the east-west electric field east and west of Jupiter can 
therefore modulate the radial inward displacement of Io relative to the loca- 
tion of the longitudinally asymmetric plasma torus sink and hence modulate the 
sodium population south of the satellite accordingly. 

In the sodium lifetime description of Figure 9 (case C of Smyth and Combi 
1988b), the east-west field has been assumed simply to be constant. Studies 
of Voyager EUV ion emission data (Sandel 1985), however, suggest that the 
electric field may be time variable, but that it also has a preferred time- 
average value directed from east- to-west , as assumed in Figure 9. More recent 
studies of the plasma ion emission data (Sandel and Dessler 1987) show that 
these emissions exhibit a clear and persistent dual periodicity (System III, 
and a new System IV differing by almost exactly 3% in period). Sandel and 
Dessler propose that (1) the convection electric field across Jupiter's 
magneto tail is limited primarily to the west side of the planet, where the 
overflow of plasma from the Io torus occurs, and (2) that this electric field 

is stronger and also modulated at the System IV period west of Jupiter, but is 

weaker and shows little System IV modulation east of Jupiter. This causes the 
plasma torus west of Jupiter (as is evident in ion emission data) to be 
shifted on average closer to the dipole offset of the magnetic field by 
~0.3 Rj (similar to the case C in Figure 9) but shifted in a different manner 
than by equal radial contributions by a constant electric field east and west 
of Jupiter. The assumption of a constant electric field as adopted in 
Figure 9 would then need to be modified to properly include one shift con- 
tributed by the electric field west of Jupiter and a different shift con- 
tributed by the electric field east of Jupiter. The presence of an unequal 
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electric field strength east and west of Jupiter would then change the abso- 
lute location of the dark ovals in Figure 9 but not their time -average rela- 
tive separation. Small changes in the absolute location of the east dark oval 
could then produce different modulations of the sodium abundance south of Io 
and thereby provide an explanation for the apparent erratic behavior of the 
east image data. The investigation of these physical processes is, however, 
beyond the scope of the present project and should be undertaken in future 
research efforts for Io. 

4.3 Summary of Physical Insights 

Physical insights that have been gained in this project are briefly sum- 
marized below. These advances in our understanding are divided for discussion 
purposes into three topics: (1) the local structure of Io's atmosphere, (2) 

the gas escape mechanism at Io, and (3) the properties and stability of 
Jupiter's magnetosphere. 

Local Structure of Io's Atmosphere 

From modeling analysis of the basic morphological characteristics of the 
sodium cloud (Smyth and Combi 1988b, in Appendix VII), a number of conclusions 
were formed. The collisional exobase for the atmosphere is somewhat below 
2600 km in radius (i.e. somewhat below an altitude of 785 km). Such heights 
for the collisional exobase are nearer the "minimum" SO 2 atmosphere for Io at 
the Voyager 1 encounter required by photochemical models of Kumar (1984). 

Above the collisional exobase, the exosphere (or corona) of the satellite is 
dominated by a nonthermal upward flux distribution with a sodium density of 
-3 x 10 cm’ at a radius of 2600 km and a radial sodium density profile in 
agreement with recent measurements of Schneider et al. (1987). The nonthermal 
sodium atoms escape from the Lagrange sphere of the satellite (~10^ km in 
radius) and populate the sodium cloud. The nonthermal energy is supplied by 
ion-neutral interactions (described in the following subsection) which occur 
near the collisional exobase and in the corona. The nonthermal collisional 
processes will penetrate energetically somewhat below the classical collision- 
al exobase level and may modify significantly the gas density and simple dif- 
fusive equilibrium (vertical stratification) assumptions used in photochemical 
models for Io (Kumar 1984, 1985). In this manner, the simple photochemical 
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picture could be replaced in the upper atmosphere by a nonthermal heating 
layer and an effective global atmospheric blowoff process. 

From analysis of Smyth and Combi (1988b) and new results presented in 
this report, additional information about the local atmosphere is suggested. 
The local atmosphere is likely collisionally thick above the satellite surface 
because of the nature of the speed of dispersion of the sodium flux distribu- 
tion. The atmosphere is also likely enhanced on the sunlit hemisphere so as 
to produce a different angular distribution for the sodium flux distribution 
east and west of Jupiter as suggested by image data. Additional analysis is 
required to investigate more fully this very interesting latter point. 

Gas Escape Mechanism 

The gas escape mechanism at Io is driven by the corotating plasma through 
ion-neutral momentum transfer. This momentum transfer most likely occurs in a 
collisional thick atmosphere. In the atmosphere, ion-neutral interactions are 
characterized by an incomplete collisional cascade process described reason- 
ably well by a nonthermal modif ied-sputtering flux distribution that has been 
explicitly developed (see Smyth and Combi 1988b). The nonthermal flux dis- 
tribution dominates a thermal (1000° -2000°K) flux disbribution by roughly an 
order of magnitude in the exosphere (or corona) , and a thermal flux distribu- 
tion by itself is incompatible with the sodium cloud data. The most probable 
speed of the nonthermal flux distribution for sodium atoms is -3 km sec"^. 

The speed disperson of the nonthermal flux distribution is more appropriately 
described by a Thomas -Fermi differential cross section with an incomplete col- 
lisional cascade process than by a classical sputtering distributing for a 
complete collisional cascade process. 

The angular distribution of the nonthermal flux distribution is likely 
intermediate between an isotropic and band source. There also appears to be 
an east-west asymmetry in the angular distribution of the sodium source such 
that the distribution is peaked more in the direction of Io's orbital motion 
in the east than in the west. This accounts for the suggestion above that the 
atmosphere may well be enhanced more above the sunlit hemisphere. A more 
exact understanding of the angular distribution of the flux distribution and 
its east-west asymmetries will require additional modeling of image data. The 
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total source strength of the upward directed sodium for the nonthermal flux 
distribution is in the range of -1-3 x 10 4 atoms sec’ and depends upon the 
details of the angular distribution of the sodium source. 

Magnetosphere of Jupiter 

A 7° tilted and corotating heavy- ion plasma torus with a planetary offset 
dipole magnetic field and Voyager-like radial profiles from the plasma density 
and temperatures is essential in understanding the basic morphology of the so- 
dium cloud. The predominance and length of the forward sodium cloud are de- 
termined naturally within this plasma torus description by means of the re- 
sulting spactime lifetime structure produced for sodium atoms. The north- 
south alternating nature of the directional feature in the trailing cloud is 
also naturally explained by this spacetime lifetime behavior and higher speed 
(-15-20 km sec’l) sodium atoms in the nonthermal flux distributions. The ba- 
sic structure of the plasma torus can be judged to be reasonably stable for 
more than ten years by the consistent presence of these identified sodium sig- 
natures in the cloud data over this time period. 

An east-west electric field with a time-average value of -2.8 mV/m is 
present in the magnetosphere „ This electric field produces the east-west in- 
tensity asymmetry in the sodium cloud near Io (Smyth and Combi 1987, 1988b) 
and also the anti-correlated east-west intensity asymmetry observed both in 
optical and ultraviolet ion emissions in the torus (see Smyth and Combi 
1987). Both sodium data and ion emission data suggest that the strength of 
this east-west electric field east and west of Jupiter may be different and 
may undergo different time variabilities. 

The plasma torus contains an inherent System III longitudinally asym- 
metric enhancement in the active sector (-180° -320° ) . This has been noted in 

ground based observations of the torus ion emissions. This System III asym- 

metric enhancement has been shown to produce a periodic south "bite -out" in 
the west sodium cloud which has been present in the various sodium data sets 
for at least seven years. A zone of deficient sodium south of Io in the east 

cloud is sometimes observed and sometimes not observed. It appears that this 

time variable or erratic behavior may be linked to a time variability in the 
east-west electric field as discussed above in Section 4.2. 
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4 . 4 Future Studies 


As noted in the introduction, the complete JPL Io sodium cloud data set 
is composed of 1974-1979 spectral data, 1981 spectral data, 1976-1979 Region 
B/C image data, 1981 Region B/C image data, and 1981 Region A image data. 
Studies undertaken in this project have considered only a subset of the 1981 
Region B/C image data. At this time, a number of very important studies for 
the 1981 Region B/C image data still remain to be considered. Studies for the 
other components of the data set have yet to be even initiated. Some import- 
ant aspects of future studies are briefly described below. 

Additional studies of the 1981 Region B/C image data could be most easily 
undertaken by extending and refining the analysis already begun in this proj- 
ect. For the local structure of Io's atmosphere, the radial density structure 
of sodium may be readily probed beyond the occulting mask and questions as to 
the enhancement of the gas density on the sunlit hemisphere straightforwardly 
investigated by comparison and analysis of east and west images. The effects 
of unequal east-west electric fields east and west of Jupiter may likewise be 
pursued by further comparison and. analysis of east images acquired on differ- 
ent days. This latter analysis will lead to a better understanding of the 
structure, stability, and time variability of the planetary magnetosphere and 
their connection to the Iogenic plasma source. 

The analysis of the 1981 Region B/C images may also be readily applied to 
the 1976-1979 Region B/C images. This will allow us to extend the time line 
and to test for the stability or time variability of various morphological 
features of the cloud and their corresponding physical properties and struc- 
tures in the plasma torus. The 1976-1979 Region B/C data set includes two 
images on the Voyager 1 encounter day as well as a number of images around 
this encounter time. Study of these images would be most valuable in linking 
the whole sodium cloud data set to Voyager 1 encounter conditions determined 
by the spacecraft. 

Analysis of the 1974-1979 spectral data, 1981 spectral data, and the 1981 
Region A image data would also be valuable for refining our understanding of 
Io's local atmosphere and its stability or time variability. Studies of these 
data can be used to extend the 1974-1975 measurements for the east-west 
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asymmetry near Io of Bergstralh et al. (1975, 1977) to the 1981 time frame. 
Radial profiles for the sodium density within the Lagrange sphere and beyond 
could be extracted from these data and analysis applied. These studies would 
also be most valuable in understanding the stability of lo's atmosphere, the 
plasma source, and the magnetosphere over a seven year period. 
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Appendix I 


Fully Reduced 1981 Sodium Cloud Images: 


Study 1 



Appendix I 


For each observation, the D 2 image (upper contour plot) is shown in 
proper relationship and scale to Jupiter, Io, and Io's orbit (as seen by an 
Earth observer) while this information is suppressed for the simultaneously 
measured but dimmer image (lower contour plot). Eastern elongation for Io 
is to the left, western elongation is to the right, and north is up. Contour 
brightness levels, from outside to inside, have increasing values that are 
spaced as follows: 0.1, 0.2, 0.5, 1, 2, 5, 10 and 20 kiloRayleighs . In all 
cases the sodium very near Io is occultated by a circular mask (9 arc sec in 
diameter) centered on the satellite. This mask accounts for the decrease in 
the brightness of the cloud very near Io. 

For each observation, the Image Number (Tape/Frames), UT date, and mid- 
point condition for the Io geocentric phase angle and the Io System III longi- 
tude angle (from Table 2) are given. The number in parentheses following the 
UT date identifies the first, second, third, etc. image fully reduced on that 
day as ordered in Table 2 and indicated by dots in Figure 1. 



Appendix I 

1981 Region B/C Images: Processed by MIPL 


Start Conditions End Conditions Mid-Point Conditions 


Date 

of 

Observations 

Time 

(UT) 

lo 

Phase 

Angle 

(deg) 

Magnetic 
Longitude 
of lot 
(deg) 

Time 

(UT) 

Io 

Phase 

Angle 

( de g) 

6 April 

5:00 

291.1 

303.2 

5:36 

296.3 

28 April 

4:49 

88.0 

213.1 

5:24 

92.9 


7:28 

110.3 

286.9 

7:58 

114.5 


7:39 

111.9 

292.0 

8:07 

115.8 

29 April 

3:48 

283.5 

131.3 

4:19 

287.9 


4:22 

288.3 

147.0 

4:53 

292.7 

4 May 

3:25 

217.6 

216.2 

3:55 

221.8 


4:20 

225.4 

241.6 

4:53 

230.1 


4:56 

230.5 

258.2 

5:33 

235.8 


5:49 

238.1 

282.7 

6:28 

243.6 


7:19 

250.9 

324.3 

7:52 

255.6 


9:06 

266.1 

13.7 

9:50 

272.4 

5 May 

3:59 

65.9 

178.9 

4:30 

70.3 


4:32 

70.5 

194.2 

5:15 

76.6 


6:34 

87.7 

250.9 

7:10 

92.7 


8:03 

100.2 

292.2 

8:35 

104.7 

6 May 

5:36 

283.6 

170.4 

6:10 

288.5 


6:11 

288.6 

186.6 

6:43 

293.2 

12 May 

3:28 

46.4 

153.5 

3:55 

50.2 


4:10 

52.3 

173.0 

4:40 

56.5 


5:05 

60.0 

198.5 

6:28 

71.6 


6:30 

71.9 

238.0 

7:09 

77.4 


8:28 

88.5 

292.7 

9:07 

93.9 

13 May 

3:28 

250.1 

100.3 

4:03 

255.1 


4:43 

260.8 

134.9 

5:04 

263.7 


5:42 

269.1 

162.2 

6:07 

272.7 


6:09 

273.0 

174.7 

6:47 

278.4 


6:49 

278.7 

193.2 

7:23 

283.5 


7:25 

283.8 

209.8 

7:58 

288.5 


8:00 

288.7 

226.0 

8:35 

293.7 


8:25 

292.3 

237.6 

9:00 

297.2 

6 June 

4:19 

100.2 

292.4 

4:53 

105.0 

14 June 
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1981 Region B/C Images: 

East-West Comparisons on Consecutive Days 
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Appendix II 


Fully Reduced 1981 Sodium Cloud Images: 


Study 2 



Appendix II 


1981 Region B/C Images with 

Recurrence of Io Geocentric Phase and Io System III Longitude Angles 


Mid-Point Condition 
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Appendix III 


Fully Reduced 1981 Sodium Cloud Images: 


Study 3 


Appendix III 


1981 Region B/C Images: 

Comparison for Diametrically Opposite Io Phase Angles 


Mid-Point Condition 
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Image ID 

Phase 

Phase 

Phase 

System III 


UT 

Number 

Angle 

+ 180° 

Angle 

Longitude 

Image 

Date 

( Tape /Frames) 

I deg) 

.(deg) 

(deg) 

(deg) 

Set 

12 May (1) 

SIP 420/6-8 

48.3 

228.3 

- 

159.8 

3. A 

4 May (2) 

SIP 417/13-15 

- 

- 

227.8 

249.2 

3. A 

12 May (2) 

SIP 420/10-12 

54.4 

234.4 . 

* 

180.0 

3.B 

4 May (3) 

SIP 417/16-18 

- 

- 

233.2 

266.8 

3 . B 

5 May (1) 

SIP 418/10-12 

68.1 

248.1 

_ 

186.1 

3 . C 

12 May (3) 

SIP 420/15-20 

65.8 

245.8 

- 

217.8 

3 . C 

4 May (4) 

SIP 417/20-22 

- 

- 

240.9 

291.7 

3 . C 

4 May (5) 

SIP 417/27-29 

- . 


253.3 

331.9 

3.C 

5 May (2) 

SIP 418/13-16 

73.6 

253.6 

_ 

204.2 

3.D 

12 May (4) 

SIP 420/21-23 

74.7 

254.7 

- 

247.1 

3 . D 

13 May (1) 

SIP 421/21-23 

- 

- 

252.6 

108.4 

3 . D 

4 May (5) 

SIP 417/27-29 
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353.3 

331.9 

3 . D 

28 Apr (1) 

SIP 415/27-29 

90.5 

270.5 
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221.3 

3 . E 

5 May (3) 

SIP 418/24-26 

90.2 

270.2 . 
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259.3 

3 . E 

12 May (5) 

SIP 420/30-32 

91.2 

271.2 

' - 

301.8 

3 . E 

4 May (6) 

SIP 417/36-39 
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23.9 

3 . E 

13 May (3) 
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168.0 
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13 May (4) 

SIP 421/34-36 
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275.7 

183.5 

3 . E 

5 May (4) 

SIP 418/31-33 

102.5 

282.5 
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299.6 

3 . F 

6 Jun (1) 

SIP 424/10-12 

102.6 

282.6 
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300.3 

3 . F 

29 Apr (1) 

SIP 416/4-6 
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285.7 

138.5 

3 . F 

6 May (1) 

SIP 419/19-21 
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286.1 

178.3 

3 . F 

13 May (6) 

SIP 421/40-42 
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286.2 
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3 . F 

14 Jun (1) 

SIP 425/6-9 
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3 . F 

28 Apr (2) 

SIP 415/43-45 

112.4 

292.4 
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3.G 

28 Apr (3) 
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- 
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3 . G 

6 May (2) 
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SIP417/27-29 Image Set: 

UT Date: May 4, 1981 (5) 
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UT Date: May 12, 1981 
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UT Date: April 28, 1981 (1) 
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UT Date: May 5, 1981 
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UT Date: May 13, 1981 
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285 . 7 ° 

SIP416/4-6 Image Set: 

UT Date: April 29, 1981 (1) 
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UT Date: May 13, 1981 
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SIP425/6-9 Image Set: 

UT Date: June 14, 1981 (1) 
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112.4* (292.4°) 

Image Set: 3.G SIP415/43-45 

UT Date: April 28, 1981 (2) 
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UT Date: April 28, 1981 (3) 
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SIP419/22-24 Image Set: 

UT Date: May 6, 1981 (2) 
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SIP421/45-47 Image Set: 

UT Date: May 13, 1981 (8) 
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UT Date: April 6, 1981 (1) 
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Appendix IV 


Fully Reduced 1981 Sodium Cloud Images: 


Study 4 



Appendix IV 


1981 Region B/C East Images: 
Comparison for Similar Phase Angles 


Mid-Point Condition 



Image ID 

Io Phase 

Io System III 


UT 

Number 

Angle 

Longitude 

Image 

Date 

( Tape /Frames) 

(deg) 

(deg) 

Set 

12 May (1) 

SIP 420/6-8 

48.3 

159.8 

- 

12 May (2) 

SIP 420/10-12 

54.4 

180.0 

- 

5 May (1) 

SIP 418/10-12 

68.1 

186.1 

4. A 

12 May (3) 

SIP 420/15-20 

65.8* 

217.8 

4. A 

5 May (2) 

SIP 418/13-16 

73.6 

204.2 

4.B 

12 May (4) 

SIP 420/21-23 

74.7 

247.1 

4 . B 

28 Apr (1) 

SIP 415/27-29 

90.5 

221.3 

4.C 

5 May (3) 

SIP 418/24-26 

90.2 

259.3 

4 . C 

12 May (5) 

SIP 420/30-32 

91.2 

301.8 

4 . C 

5 May (4) 

SIP 418/31-33 

102.5 

299.6 

4.Dt 

6 Jun (1) 

SIP 424/10-12 

102.6 

300.3 

4.Dt 

28 Apr (2) 

SIP 415/43-45 

112.4 

293.9 

_ 

28 Apr (3) 

SIP 415/44-46 

113.9 

298.5 
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critical east phase angle 
•same as Image Set 2. A 
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Appendix IV 


1981 Region B/C West Images: 
Comparison for Similar Phase Angles 




Mid- Point 

Condition 



Image ID 

Io Phase 

Io System III 


UT 

Number 

Angle 

Longitude 

Image 

Date 

( Tape /Frames') 

(dec) 

(deg) 

Set 

4 May (1) 

SIP 417/8-10 

219.7 

223.1 

_ 

4 May (2) 

SIP 417/13-15 

227.8 

249.2 

_ 

4 May ( 3 ) 

SIP 417/16-18 

233 . 2* 

266.8 

_ 

4 May (4) 

SIP 417/20-22 

240.9 

291.7 

- 

13 May (1) 

SIP 421/21-23 

252.6 

108.4 

4. E 

4 May (5) 

SIP 417/27-29 

253.3 

331.9 

4 . E 

13 May (2) 

SIP 421/27-28 

262.3 

139.8 

- 

4 May (6) 

SIP 417/36-39 

269.3 

23.9 

4.F 

13 May (3) 

SIP 421/32-33 

270.9 

168.0 

4 . F 

13 May (4) 

SIP 421/34-36 

275.7 

183.5 

. 

13 May (5) 

SIP 421/37-39 

281.1 

201.1 

- 

29 Apr (1) 

SIP 416/4-6 

285.7 

138.5 

4.G 

6 May (1) 

SIP 419/19-21 

286.1 

178.3 

4.G 

13 May (6) 

SIP 421/40-42 

286.2 

217.5 

4.Gt 

14 Jun (1) 

SIP 425/6-9 

286.9 

222.1 

4 . Gf 

29 Apr (2) 

SIP 416/7-9 

290.5 

154.2 

4.H 

6 May (2) 

• SIP 419/22-24 

290.9 

194.0 

4.H 

13 May (7) 

SIP 421/43-45 

291.2 

234.1 

4.H § 

14 Jun (2) 

SIP 425/9-12 

291.6 

237.6 

4.H § 

13 May (8) 

SIP 421/45-47 

294.8 

245.7 

4.H 

6 Apr (1) 

SIP 410/13-15 

293.7 

311.5 

4.H 

critical west 
Tsame as image 
S same as image 

phase angle * 
set 2.B 
set 2.C 
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SIP421/32-33 Image Set: 

UT Date: May 13, 1981 (3) 
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SIP419/22-24 Ima e e Set: 4 H 

UT Date: May 6, 1981 (2) 
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Appendix V 


Fully Reduced 1981 Sodium Cloud Images: 


Study 5 



Appendix V 


1981 Region B/C East Images: 

Comparison for Similar Io System III Longitude Angles 


Mid- Point Condition 



Image ID 

Io System III 

Io Phase 


UT 

Number 

Longitude 

Angle 

Image 

Date 

(Tape /Frames) 

(dec) 


Set 

12 May (1) 

SIP 420/6-8 

159.8 

48.3 

_ 

12 May (2) 

SIP 420/10-12 

180.0 

54.4 

- 

5 May (1) 

SIP 418/10-12 

186.1 

68.1 

- 

5 May (2) 

SIP 418/13-16 

204.2 

73.6 

- 

12 May (3) 

SIP 420/15-20 

217.8 

65.8 

5. A 

28 Apr (1) 

SIP 415/27-29 

221.3 

90.5 

5. A 

12 May (4) 

SIP 420/21-23 

247.1 

74.7 

- 

5 May (3) 

SIP 418/24-26 

259.3 

90.2 

- 

28 Apr (2) 

SIP 415/43-45 

293.9 

112.4 

- 

28 Apr (3) 

SIP 415/44-46 

298.5 

113.9 

5.B 

5 May (4) 

SIP 418/31-33 

299.6 

102.5 

5 . B 

6 Jun (1) 

SIP 424/10-12 

300.3 

102.6 

5 . B 

12 May (5) 

SIP 420/30-32 

301.8 

91.2 

5 . B 
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Appendix V 


1981 Region B/C West Images: 

Comparison for Similar Io System III Longitude Angles 


Mid-Point Condition 



Image ID 

Io System III 

Io Phase 


UT 

Number 

Longitude 

Angle 

Image 

Date 

(Tape/Frames) 

(deg) 

(deg) 

Set 

4 May (6) 

SIP 417/36-39 

23.9 

269.3 


13 May (1) 

SIP 421/21-23 

108.4 

252.6 

- 

29 Apr (1) 

SIP 416/4-6 

138.5 

285.7 

5 . C 

13 May (2) 

SIP 421/27-28 

139.8 

262.3 

5 . C 

29 Apr (2) 

SIP 416/7-9 

154.2 

290.5 

_ 

13 May (3) 

SIP 421/32-33 

168.0 

270.9 

- 

6 May (1) 

SIP 419/19-21 

178.3 

286.1 

- 

13 May (4) 

SIP 421/34-36 

183.5 

275.7 

- 

6 May (2) 

SIP 419/22-24 

194.0 

290.9 

- 

13 May (5) 

SIP 421/37-39 

201.1 

281.1 

- 

13 May (6) 

SIP 421/40-42 

217.5 

286.2 

5 . D 

14 Jun (1) 

SIP 425/6-9 

222.1 

286.9 

5 . D 

4 May (1) 

SIP 417/8-10 

223.1 

219.7 

5 . D 

13 May (7) 

SIP 421/43-45 

234.1 

291.2 

5 . E 

14 Jun (2) 

SIP 425/9-12 

237.6 

291.6 

5 . E 

13 May (8) 

SIP 421/45-47 

245.7 

294.8 

5 . F 

4 May (2) 

SIP 417/13-15 

249.2 

227.8 

5 . F 

4 May (3) 

SIP 417/16-18 

266.8 

233.2 

_ 

4 May (4) 

SIP 417/20-22 

291.7 

240.9 

- 

6 Apr (1) 

SIP 410/13-15 

311.5 

293.7 

- 

4 May (5) 

SIP 417/27-29 

331.9 

253.3 
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ABSTRACT 


A general mathematical formalism for calculating the physical properties 
of any of Io's neutral gas clouds (Na, K, 0, S, S0 2 ) is presented. The 
dynamical effects of both the gravitational fields of Io and Jupiter and solar 
radiation pressure are included, in addition to the many complex space- and 
time -dependent interactions that occur between the neutral clouds and the 
plasma torus. Earlier models have not included these complex space- and time- 
dependent interactions. The importance of this new model in studying both 
the plasma conditions prevalent in the inner planetary magnetosphere and 
the nature of Io's local atmosphere is discussed. A numerical method for 
evaluating the physical properties of the neutral clouds using the new model 
has been developed and is described. The mathematical formalism and numerical 
method are applied to a study of the Io sodium cloud in a companion paper. 
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1. Introduction 

During the past two decades, Io, the innermost of the four Galilean 
satellites of Jupiter, and its profound influence on the planetary magneto- 
sphere have been one of the most exciting and intensely studied subjects in 
the solar system. The first evidence that Jupiter and Io constituted a com- 
plex interactive system was introduced by Bigg (1964), who discovered that 
some of the decametric radio emission of Jupiter was directly influenced by 
the satellite's orbital and magnetic location about the planet. Since then, 
a large information base has been accumulated from ground-based observations, 
rocket measurements, Earth-orbiting satellites, and four spacecraft encounters 
with Jupiter: Pioneer 10 (4 December 1973), Pioneer 11 (2 December 1974), 

Voyager 1 (5 March 1979), and Voyager 2 (9 July 1979). 

From ensuing studies of this information base , a picture for the Io- 
Jupiter system is emerging and slowly coming into focus. Io has a local atmo- 
sphere, although it is not well characterized or structurally understood. Io 
has extended neutral atmospheres or clouds that exist large distances beyond 
the gravitational grasp of the satellite and fill partial or somewhat incom- 
plete toroidal -shaped volumes around Jupiter. These clouds occupy an overlap- 
ping volume with the Io plasma torus - the densest part of the planetary 
magnetosphere - and interact with it through collisional processes including 
ionization and charge exchange reactions. These interactions not only shape 
the spatial morphology of the clouds but also supply heavy ions to the torus 
and affect the composition, structure, and energy budget of the plasma torus 
and its extension into the larger magnetosphere. The primary heavy ions that 
compose the corotating (or nearly corotating) plasma torus are known to be 
oxygen (0 + , 0 ++ ) and sulfur (S + , S ++ , S +++ ) , although the relative abundances 
of these ions throughout the torus volume (even for a two-dimensional 
longitudinally averaged picture) are still not well established. Io and its 
local atmosphere also interact with the corotating plasma torus and planetary 
magnetic field in ways that are yet to be fully understood. Some of these 
interactions almost certainly provide ions and energy directly to the plasma 
torus as well as drive escape mechanisms for gases in the local atmosphere. 

For purposes of this paper, a brief review of Io's local and extended 
atmospheres and their relevant interactions with the plasma torus will be 
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given in section 2. This will provide necessary information for understanding 
the evolution and motivation for developing the new general neutral gas cloud 
model presented in sections 3 and 4. A more complete discussion of Io may be 
found in several excellent review articles (Nash et al. 1986; Brown, Pilcher 
and Strobel 1983; Fanale et al. 1982, 1977; Kumar and Hunten 1982; Pilcher and 
Strobel 1982; Brown and Yung 1976). The paper concludes with a brief summary 
in Section 5. The new model is applied to the sodium cloud in a companion 
paper (Smyth and Combi 1988). In a series of future papers, the new model 
will also be applied to other neutral clouds of Io (K, 0, S, S0 2 ). 
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2. Historical Perspective 

Our understanding of Io's local atmosphere, although far from complete, 
has grown steadily since 1964, when Binder and Cruikshank (1964) reported an 
anomalous brightening of the satellite's surface as it emerged from eclipse. 
This posteclipse brightening suggested a dynamic and volatile atmosphere 
for Io, for which in 1971 an upper limit surface pressure of approximately 
IO' 7 bars was deduced by Smith and Smith (1972) from occultation of the star 
£ Scorpii by the satellite. In the summer of 1972, the discovery by Brown 
(1974) of sodium optical emission from Io's vicinity and the subsequent obser- 
vations in 1973 by Trafton, Parkinson, and Macy (1974) of an extended sodium 
cloud well beyond the gravitational reach of Io confirmed the presence of a 
local atmosphere. This discovery also sparked a large number of Earth-based 
observations of the sodium cloud that have continued over the intervening 16 
years. An ionosphere for Io discovered from analyses of Pioneer 10 S-band 
radio occultation data (Kliore et al. 1974) provided additional information to 
characterize the satellite's local atmosphere and its possible interactions 
with the Jovian magnetosphere (Cloutier et al. 1978). A quantum leap forward 
in understanding Io's local atmosphere was obtained by the Voyager 1 space- 
craft's discovery of active volcanic plumes (Morabito et al. 1979) and 
detection of SC >2 gas (~10~ 7 bars surface pressure) over one of many hot spots 
on Io (Pearl et al. 1979). Although SC >2 now appears to be the major consti- 
tuent of the atmosphere (Kumar 1979), the vertical structure, global extent, 
and surface density variations of the atmosphere still remain uncertain even 
after much theoretical analysis (Kumar 1980, 1982, 1984, 1985; Matson and Nash 
1983; Summers 1984; Ingersoll, Summers, and Schlipf 1985; Lunine and Stevenson 
1985; Summers et al. 1986). Insightful information for the density of sodium 
in Io's local atmosphere is just beginning to emerge, based upon analysis of 
novel Earth-based observations recently obtained by Schneider et al. (1985, 
1987) of absorption signatures at the sodium D-line wavelengths of reflected 
sunlight from Europa as it was eclipsed by Io. 

The discovery of the Io sodium cloud marked the first of four extended 
atmospheres that have been detected to date for the satellite as summarized in 
Table 1. The three more recently detected neutral clouds for Io are potassium 
in 1975 and atomic oxygen in 1980, both from ground-based observations of 
optical lines, and atomic sulfur in 1981 from rocket-borne instruments 
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measuring ultraviolet emissions. A fifth cloud of neutral SO 2 is expected, 
because of the presence of a small amount of SC> 2 + ions measured in the plasma 
torus (Bagenal and Sullivan 1981; Bagenal 1985), but has not been detected 
directly. For sodium, more than 1000 observations of its spatial brightness 
have been made in the last 16 years, with most of the earlier measurements 
providing only slit- averaged intensity values and most of the more recent 
measurements providing increasingly higher -quality two-dimensional image 
data. For the dimmer potassium cloud (about one- twentieth the brightness of 
the sodium cloud) , about 200 slit-averaged observations have been made almost 
entirely by Trafton (1975, 1981), who presented a survey of this gas cloud 
that exhibited spatial and temporal variations similar to those documented for 
the Io sodium cloud. For the much dimmer atomic oxygen and atomic sulfur 
clouds, only a very few slit-averaged detection observations have been 
acquired. 

The presence of neutral clouds for Io provides a largely unexplored but 
rich potential for studying both the plasma conditions prevalent in the inner 
planetary magnetosphere and the nature of Io's local atmosphere. Sodium is by 
far the best neutral cloud to use as a probe for these purposes, because it is 
the brightest and most extensively observed and is (together with potassium) 
sensitive to spatial and temporal changes in the plasma torus because of its 
small electron impact ionization lifetime (~ 1-2 hours) near Io's orbit. 

Early observations acquired in 1973 and 1974 (i.e., pre- and post-Pioneer 10) 
by Macy and Trafton (1975 a,b) and Munch and Bergstralh (1977), for example, 
first suggested the presence of a magnetospheric sink for sodium near Io's 
orbit that was both radially asymmetric and oscillatory. The radial asymmetry 
was suggested by the presence of a predominantly forward sodium cloud, and the 
oscillatory nature was suggested by time dependent changes in the sodium 
brightness north and south of the satellite that were correlated with the 
satellite's magnetic latitude (Trafton and Macy 1975; Trafton 1980). As will 
be shown in the companion paper, the predominantly forward cloud is indeed 
naturally produced by the radial asymmetry of the plasma torus for sodium 
ejected from Io at relatively slow speeds (~ 3 km sec"^) as first suggested 
from these early observations. These slow speed atoms may be energized by 
magnetospheric ions both directly through scattering of neutrals from the 
local atmosphere (or surface, if regions of low density exist) as discussed in 
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a number of papers (Sieveka and Johnson 1984; Cheng 1984; Lanzerotti and Brown 
1983; Brown, Pilcher, and Strobel 1983; Ip 1982; Matson, Johnson, and Fanale 
1974), and indirectly through multicollisional processes resulting in atmos- 
pheric sputtering (Summers, Yung, and Haff 1983; Haff, Watson, and Yung 1981; 
McGrath and Johnson 1987) or blowoff (Hunten 1985). More recent observations 
of the sodium cloud have shown, for example, (1) well-defined spatial zones 
that become temporally deficient of atoms, perhaps as a result of locally 
(likely System III) enhanced ionization (Goldberg, Garneau, and LaVoie 1984), 
(2) directional features produced by high-velocity (- 20 km sec' 1 ) sodium 
driven elastically by the magnetospheric wind from the local atmosphere at 
near right angles to Io's orbit and preferentially from the equatorial regions 
(Pilcher et al. 1984), and (3) very fast (< 60 km sec' 1 ) narrowly collimated 
jets of sodium driven directly ahead of the satellite by charge exchange reac- 
tions with the plasma torus (Trauger 1984, 1985). Although the oxygen and 
sulfur clouds are only bright enough to be barely detected, physical insight 
acquired for the bright sodium cloud may be adopted where appropriate in 0, S 
and S0 2 cloud models. These atomic and molecular cloud models so constructed 
then become enhanced tools themselves to determine ion loading rates , plasma 
mass loading rates, and ion energy input rates to the plasma torus as well as 
to specify the character and source rate of neutral gases escaping from Io's 
local atmosphere (Smyth and Shemansky 1983). 

Models for the neutral clouds of Io are essential in probing and re- 
covering physical information from spatial brightness observations of Io's 
extended atmospheres and relating this information to conditions in Io's local 
atmosphere and the planetary magnetosphere. The development of cloud models 
for Io was historically initiated for the sodium cloud in the early 1970 's and 
has evolved during the last decade with progressively more complete descrip- 
tions for the source, lifetime, and orbital dynamics of sodium atoms in the 
Jupiter system. Only in the last few years have models for the Io oxygen and 
sulfur clouds (Smyth and Shemansky 1983) and the Io potassium cloud (Smyth and 
Combi 1984) been developed. A general model, appropriate to any of these neu- 
tral gas clouds of Io, has now been formulated and is presented in Sections 3 
and 4. 
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3. Mathematical Description for the Io Neutral Gas Clouds 
a ) Mathematical Approach 

The basic mathematical approach is based upon the conservation in 
phase space of the one-particle distribution function ftf.I.t) for a gas 
species. In brief, if at an initial time t - ^ the distribution function can 
be specified, then the distribution function at a later time may, under cer- 
tain circumstances of interest here, be simply related to its initial value by 
properly mapping its time evolution in phase space and Including processes 
that modify its initial value during this time evolution. This conservation 
in phase space follows directly from the well-known Liouville theorem for an 

N-particle system upon integrating over the six phase-space coordinates of N-l 
of these particles. 

The general approach may be formally understood as follows. The one - 
particle distribution function in kinetic theory has the physical meaning that 


f (v, x , t) dx dv 


(i) 


represents the (ensemble-average) number of atoms (or molecules) at time t in 

a spatial volume element dj about position J, with velocities A about 5 . The 

evolution of the one-particle distribution function is governed by the 
kinetic equation 


8f - 

at + v 


af f 

dx m 


3v } 


( 2 ) 


where F, which may in general depend on (v.x.t), is the force experienced by 
gas atoms (or molecules) due to macroscopic internal and external sources 
the mass of a gas^atom (or molecule). Furthermore, (Sf/St) is the 
rate of change of f (v.x.t) due to processes such as collisions that may 
add or remove the (ensemble-average) number of atoms (or molecules) from an 
element of phase space. For binary collisions in a gas, for example, one 
obtains the Boltzmann equation when the right-hand side of the kinetic 
equation is represented by the Boltzmann collision integral (a nonlinear and 
integral function of the distribution function). 
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The time evolution or trajectory of a point in phase space (v,x) - 
(v(t),x(t)) that is initially (t-t^) at a point ~ ( v (t b > ,x( t fe ) ) is 

determined by solving the equations of motion for a gas atom (or molecule) of 
mass m 


v(t) - ^ 

F(v(t),x(t),t) _ dv(t) 
m dt 


(3. a) 


(3.b) 


for this initial value problem. These equations of motion are by construction 
the characteristic equations of the kinetic equation, so that the left-hand 
side of the kinetic equation is the total time derivative of the distribution 
function. The kinetic equation thus reduces to 


df ( Sf 
dt " ( St } * 


(4) 


There exists a class of physical problems to which the lo gas clouds may 
be related in which the gas density is sufficiently low that gas atoms (or 
molecules) along their phase space trajectories do not collide with each other 
but may collisionally interact with other material (such as ions and electrons 
in a rarefied plasma) or photons (i.e., photoionization or photodissociation 
processes) so as to change the initial ensemble -average number of atoms; In 
this case, the right-hand side of the kinetic equation is linear in the 
distribution function and may be expressed in the general form 


< ft £ - s , (5) 

where v and S are independent of f but may depend in general upon (v,x,t). 
Here v is an effective loss rate and S is an effective source rate for the 
ensemble -average number of atoms in an element of phase space. Using ex- 
pression (5) in equation (4), the kinetic equation can be explicitly 
integrated along the phase space trajectory to determine the distribution 


function at time t: 



10 


f(v(t) ,x(t) ,t) - f(v b ,^,t b ) P(t,t b ) 

t 

+J dt' S(v(t'),x(t'),t') P(t.t') . 

% 

Here 


( 6 ) 


PCt^.t^) - exp 


-J 2 dt' i/(v(t'),x(t'),t') 
C 1 


(7) 


is the probability that an ensemble -average number of gas atoms introduced at 
a time t^ exist at a later time t 2 because of the effective loss rate u. The 
distribution function at a time t > t b is then given by the initial distribu- 
tion function reduced by this probability plus all new source contributions to 
the initial distribution function that have occurred in the time interval t b 
to t with a proper decay loss included for these contributions. 


The description of the gas may therefore be constructed in terms of a 
large number of phase space trajectories along which the ensemble -average 
number of atoms evolve in time based upon a set of initial conditions for the 
distribution function which properly simulate the boundary conditions of a 
particular physical problem. More specifically, the description of the gas 
can be fundamentally reduced to describing the time evolution of an ensemble 
packet of gas along its phase space trajectory. An ensemble packet is defined 
as the infinitesimal volume element containing the ensemble -average number of 
atoms that initially began its trip along the trajectory in the time interval 
t b to (t b + dt) . The physics of an ensemble packet appropriate for the Io 
neutral gas clouds is described in the next subsection. 


The basic principle of the general mathematical approach described above 
has been employed in planetary science to develop models for a number of 
studies involving collisionless or nearly collisionless atmospheres . The 
basic principle has been, in fact, used with increasing regularity in the last 
two decades. These studies have included coraetary atmospheres, planetary 
atmospheres, and satellite atmospheres for which major references are given in 
Table 2. A more general review of exospheric theories is given by Fahr and 
Shizgal (1983). The studies for Io's neutral gas clouds are perhaps the most 
complex subset of these exospheric problems appropriate to the linear inter- 
action (5). This complexity results because of the inherent velocity, space, 



11 


and time dependence of both the force F in the equations of motion and the 
effective loss rate v of neutrals in the plasma torus and also because of the 
form of the source term S which is described below. 

For the Io neutral gas cloud, the source term S in equation (5) is a 
spatially distributed source that describes the creation of secondary cloud 
atoms produced by interactions of the raagnetospheric plasma with only the 
primary gas cloud atoms (or molecules) initially ejected from the satellite. 
For atomic gas clouds, the distributed source is caused by charge exchange 
reactions between the extended cloud atoms and the corotating plasma torus 
ions. Elastic collisions between these ions and extended cloud atoms may also 
produce a contribution to the source term in equation (5) but this contri- 
bution may be significantly less important over the relatively short lifetime 
of the cloud. Primary cloud atoms have velocities relative to Io that are 
generally slow compared to the corotational ion speed (-57 km sec and, when 
removed by the atom- ion charge exchange reactions, create fast secondary cloud 
atoms that have a rather tightly confined velocity dispersion centered on the 
corotational speed. The fast secondary cloud atoms are thus created in a 
velocity volume of phase space that is generally discontinuous with the velo- 
city volume of phase space containing the slow primary cloud atoms. For 
molecular gas clouds (e.g. , SO 2 or SO), electron impact dissociative reactions 
very rapidly remove the primary cloud molecules (Shemansky 1984) and create 
other secondary cloud atoms (or molecules) that initially may have additional 
amounts of kinetic energy as dictated by the energetics of the reaction. 
Because of the rapid dissociative lifetime of molecules in the plasma torus , 
these dissociative sources are likely to be heavily confined near Io. A gas 
cloud of Io for the j*'* 1 neutral species may therefore be generally described 
by a primary and a secondary population of cloud atoms (or molecules) by 
formally dividing the one-particle distribution function (1) into two parts. 

f _ f (primary) + f (secondary) ^ 

j * j J 

The kinetic equation (4) with the right hand side given by expression (5) 
may thus be formally subdivided for the j species into equations for the 
primary and secondary atom (or molecule) populations 
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df (primary) 

l _ „ € (primary) 

dt jj 

(secondary) 

— !• - -^.f. (seconda,; y ) + s. 

dt J J J 

where the secondary source term may be written in the form 


(9. a) 


(9.b) 


s =sc f (P rimar y) / Q . 

j l c jk f k ' (9 - c) 

where C.^ describes the creation rate of secondary cloud atoms (or molecules) 
of the j species from primary cloud atoms (or molecules) of the k th 
species. The source term (9.c) therefore represents the secondary cloud atoms 
(or molecules) produced by all primary cloud atoms (or molecules). More 
generally, a source term similar to expression (9.c) may likewise be intro- 
duced if necessary to include third generation extended source terms produced 
by dissociation of secondary cloud molecules (e.g. SO -► S + 0) . Limiting the 
present development to second generation source terms, the description of a 
neutral cloud is then obtained by first solving equation (9. a) for the 
properties of the primary cloud atom (or molecule) population, 

(primary) (primary) 

fj(v(t),x(t),t) - fj(v b ,x b ,t b ) Pj(t,t b ) (10. a) 


and then, second, determining Sj and solving expression (9.b) for the 
properties of the secondary cloud atom (or molecule) population, 

(secondary) 

f. (v(t) ,x(t) ,t) - / dt" S (v(t M ) ,x(t") ,t")P. (t,t") . (10. b) 

J t J J 

s 

(secondary) 

where fjCVsiXgitg) “ 0 by selecting the initial point of the secondary cloud 
phase space trajectory ( v s » x s »tg) at a point on the primary cloud phase space 
trajectory where the secondary cloud distribution originates and is hence zero 
by definition. 

In this paper, a general model for the primary cloud atoms (or molecules) 
is explicitly developed and will include expressions to calculate the relevant 
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plasma input rates supplied to the plasma torus by the primary cloud. A model 
for the secondary cloud atoms (or molecules) is also given but is more briefly 
outlined. For simplicity, the general model developed will adopt the language 
appropriate to a gas cloud of atoms but, with slight generalization, will 
equally well apply to a molecular gas cloud. For a gas cloud of atoms, the 
primary cloud population is much brighter than the secondary cloud population 
in the general vicinity of Io's orbit, but the secondary cloud population ex- 
tends throughout and beyond the planetary magnetosphere and provides a back- 
ground density and cloud brightness into which the primary cloud population 
merges . 

b) Physics of an Ensemble Packet 

A neutral gas cloud of Io may be described for the primary atom 
population discussed above by following the spacetime history of neutral atoms 
ejected from the surface of the satellite (in low atmospheric density limit) 
or from the exosphere of the satellite (in the presence of a substantial 
atmosphere). Without loss of generality, a spherical exobase of radius Rp 

Cj 

centered on Io is adopted. For a differential surface area element 
2 

dfl ( - R £ sin$ d$d0) of the exobase centered on the angular direction 
^ ™ (©*$) $ let the neutral flux (number of atoms ejected per unit area, 
per unit initial velocity interval, per unit time) in the initial velocity 
interval tj to w + dw and at the absolute time t^ be denoted by ^(w,0;t^). 

The rate at which these atoms are ejected from the surface element is then 

^ A 

Vf t b ) dfl dw. (11a) 


In terms of the one-particle distribution function of the preceding subsec- 
tion, the flux <f> is defined by 

A A (primary) 

<Ku, 0, t b ) - w • 0 f(w, x b ,t b ) (11. b) 

where x b ~ (R £ , 0, $) is the spherical polar coordinate location of the surface 
area element dfl. 

To describe the neutral cloud and its impact on the planetary magneto- 
sphere, it is sufficient to follow in time an ensemble packet of atoms for 
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each neutral species initially ejected in an absolute time interval t b to 
fc b + > w ^ er e t' is defined at the flight time since the packet was initial 

ly ejected from the exobase at time t b . For one such neutral species, the 
motion of a packet is depicted in Figure 1. A number of physical quantities 
will be defined for the packet in this section. The complete description of 
a neutral cloud for a particular physical quantity is obtained by integrating 
the packet contribution for this quantity over its complete flight time and 
over all initial exobase conditions. To develop this description, the loca- 
tion of the center of the packet at flight time t' may be denoted by 
r(t',|w,0) in a coordinate system fixed relative to the satellite. In the 
case of a sodium or potassium cloud, the additional force of solar radiation 
pressure experienced by the packet atoms as they resonantly scatter sunlight 
requires a more complex solution for the packet location (see Smyth 1979, 

1983) that depends in addition on the present value of the absolute time, 
t b + t', and hence must be denoted by r(t' , t b | 3 , 2 ) . This absolute time de- 
pendence is introduced because the magnitude of the solar radiation force is 
modulated by the Doppler motion that occurs between the Sun and the packet 
atoms in the satellite frame as Io executes circular motion about Jupiter and 
as Jupiter moves in its elliptical orbit about the Sun. 

As an ensemble packet of atoms moves through the circumplanetary space, 
it may experience ionization, charge exchange, excitation, and even elastic 
collisions because of its interactions with the Io plasma torus. Neutral- 
neutral collisions are generally not important (except for recapture near Io 
and Jupiter where they are easily included) because the accumulated column 
density determined by the relatively short lifetime of cloud atoms is small. 
Electron impact ionization, photoionization, and charge exchange reactions 
remove atoms from the packet and hence determine its decay lifetime. Ioniza- 
tion processes also supply new ions to the planetary magnetosphere. Charge 
exchange processes modify the ion composition and energy of the Io plasma 
torus and also provide a flux of high-speed neutrals that (unless captured in 
the middle magnetosphere) escape from the planetary system. Electron impact 
excitation processes for neutral oxygen and sulfur atoms provide a photon 
emission mechanism for these two neutral gas clouds. Elastic ion-neutral 
collisions alter the trajectory of atoms in the packet and, if important, 
require that this momentum transfer information be properly incorporated in 
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_ A 

determining the equivalent of r(t',t^|w,0) for the secondary atom population 
as described in the previous subsection. 

To describe the interactions of the ensemble packet of atoms with the 

plasma torus, it is necessary at an absolute time t to map a location R in the 

-4 _+ A 

satellite fixed coordinate system where the packet location r(t' , t b |u,0) is 

■4 

specified to a location ft in the magnetic (or, more correctly and generally, 
what shall be called the plasma) coordinate frame where the plasma properties 
are specified. The general transformation for this mapping, which is a 
function of the absolute time t, will be denoted as follows: 

ft - T(t)*R + D(t) . (12) 

The tensor part of the transformation, T(t), includes the relative rotational 
motion and the tilt angle between the satellite coordinate frame and the 
plasma torus coordinate frame. The vector part of the transformation, D(t), 
includes the more complex mapping required between the satellite frame and the 
plasma frame for a dipole offset of the origin of the magnetic and planetary 
spin coordinate frames as well as for the additional shift of the center of 
the plasma torus produced by an east-west electric field such as proposed by 
Barbosa and Kivelson (1983) and Ip and Goertz (1983). The location of a 
packet at flight time t' in the plasma coordinate system is then given as 

?( t',t b |« f O) - T(t b + t') • r(t',t b |-w f S) + D(t b +t') . (13) 

? A 

The velocity of the packet v(t' , t b |ta,fl) at flight time t' in the plasma co- 
ordinate system is then determined directly by time differentiation of 
equation (13). 

In the remainder of this subsection, the ensemble packet contributions 

for a number of neutral cloud physical quantities of interest (e.g. density, 

photon emission rate, neutral loss rate, and various plasma input rates to 

the magnetosphere) are formulated. To facilitate this formulation, a hat 

" A " notation is adopted to indicate that a physical quantity is described 

H. A 

in terms of its phase space representation ( t ' , t b | u,fl) along its phase space 

trajectory. Thus, for a general physical quantity Q, its ensemble packet 

A A A 

contribution will be denoted by Q — Q (t',t |w,0). For an observer in a 
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coordinate frame (r^.t^), the spacetime description of the physical 
quantity Q(r^,t o ) is then determined, as discussed in a later section, 
from Q by evaluating its ensemble average < Q > . 

The number of atoms in an ensemble packet at flight time t' depend upon 
the loss processes suffered by the packet since its beginning at time t h . For 
an ensemble packet at a location r having a velocity v in the plasma coordi- 
nate frame, the lifetime of cloud atoms is denoted by r(r,v). The velocity 
dependence is required in the lifetime to calculate properly the cloud atom 
charge exchange reaction rates with the plasma torus ions. The loss collision 
frequency for the packet is then given by 


A _ A 

v(f ,t b |w,fi) 


[r(r(t' ,t b |w,Q),v(t' .tblw.O))]’ 1 , 


( 14 ) 


and the probability that the ensemble packet of atoms at flight time t' exists 
is defined as follows: 


P(t',t b |w,fl) - exp [-/ i/(t" ,t b |w,n) dt"]. 

0 

The number of atoms in the packet at flight time t' is then given by 


A . A 

n (t' ,t b |w,0) dfl dw dt' 

where 


( 15 ) 


( 16 ) 


A 

n 


(t',t b |w,n) - P(t* ,t b |w,n)^(3,fi,t b ) 


( 17 ) 


is the number of atoms per unit surface element dfl at the initial surface 
A 

location O, per unit initial velocity volume dw at the initial velocity w, and 
per unit time interval dt' at the flight time t' . 

The number of photons emitted by an ensemble packet at flight time t' 
depends upon its photon excitation mechanism. The two important mechanisms, 
electron impact excitation in the plasma torus and solar resonance scattering, 
are described by an atom excitation rate (i.e., photons per atom per unit 
time) denoted by J A (r, v g ) for an emission line with wavelength A. The posi- 
tion vector r of the packet in the plasma coordinate frame is required in 
calculating the electron impact excitation rate for oxygen and sulfur atoms. 
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The radial speed v g of the packet relative to the Sun is required in calcu- 
lating the solar resonance scattering excitation rate for sodium, potassium, 
oxygen, and sulfur. The photon emission rate (i.e., number of photons emitted 
per unit time) at wavelength A for the packet at flight time t' is then given 
by 


A A _ 

I A (t' , t b |t*>,0) dfl du dt' (18) 


where 


A A 

i A (t',t b |w,n) 


J A [r(t' ,t b |w,n) , v 8 (t',t b |«,0)] ft(t',t b |u;,Q). (19) 


The rate at which neutral atoms are lost from the packet at flight time 
t' is determined from the time derivative of expression (17) and is given by 

A 


fa P(t',t b |w,n) 


at ' 


^(w,n, t b ) dfl dtj dt' , 


( 20 ) 


which may be written in terms of a loss rate 


A -> A 

L(t ' ,t b |S,Q) dfl dtj dt' 


( 21 ) 


where 


L(t',t b |w,S) - i(t',t h \u 9 Q) S(t',t b |w > S). (22) 

Loss processes contributing to the collision frequency u in equation (22) 
occur for the first four types of neutral reactions summarized in Table 3, 
where the neutral cloud atoms are denoted by the general chemical symbol Y. 

(It should be noted in passing that to generalize the collision 
frequency u for molecular gas clouds, dissociative reactions would from this 
point forward have to be properly included.) It should be emphasized that the 
expression (21) is not , however, the net rate at which ions are deposited into 
the magnetosphere by the packet, since charge exchange reactions of type C in 
Table 3, while providing from the packet a new ion (Y + ) , simultaneously remove 
an old ion (X + ) from the magnetosphere by production of a fast escaping 
neutral . 
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In general, to determine the impact of the neutral cloud on the magneto- 
sphere, it is necessary to assess simultaneously the contributions of the 
packets of several neutral clouds on the plasma torus because of reactions of 
type C in Table 3. It is of interest to determine these contributions not 
only for the net ion loading rate N but also for the plasma mass loading rate 
M, and for the ion energy input rate E. For a mixture of N neutral gas 
species (j-1 , 2 , . . . ,N) , the instantaneous packet contributions to any one of 
these three input rates can be formulated (see Smyth and Shemansky 1983) and 
will be denoted by R(t' , t b | w,ft) , where 


A N 
R - S 

j-1 


A A 
[£ .L. 
J J 


Yi ] 


(23) 


Here Lj is the loss rate of the packet for the j th species defined in a 

parallel manner to equation (22) , and P- is the production rate of the packet 
th ^ 

for the j species as determined from the reactions of type C in Table 3. 

_ , , A A 

lor the three cases of interest, and are defined in Table 4, where V R is 
the relative velocity between the neutral atom Y and the plasma coordinate 
frame at the time the neutral is lost from the packet, and V T is the thermal 
velocity of the X ion at the time it is converted to a fast neutral because 
of reactions of type C in Table 3 . 


Combining the ionization reactions A and B in Table 3 into a single term, 

the collision frequency 0.(t\t |w,0) for the j th species may be divided into 

J b 

a contribution from ionization (I) , charge exchange (CE) , and charge transfer 
(CT) as follows: 


A . A(I) A(CE) a(CT) 

J J J j 


(24) 


The charge exchange collision frequency may be further divided into terms 
A(CE) A N A(k) A 

v j (t',t b |w,0)- ES. (t',t b |w,Q) (25) 

k-1 

that describe the individual charge exchange reactions. Here Sj^ is the loss 
rate for the j species neutral packet because of the charge exchange reac- 
tion that produces a fast k th species neutral (i.e., it is the product of the 
ion density for the ionized k-species times the reaction rate constant for the 
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charge exchange reaction). Symmetric charge exchange (k-j) and nonsymmetric 
charge exchange (k*j) are included. The production rate Pj (t' , t b | w,fl) for all 

fast neutrals of the j th species that are produced by all neutral packets is 
then given by 


A 

P. 

J 


N 

2 

k-1 


s (j) A. 

s k ”k 


(26) 


so that expression (23) may be rewritten as the sum of two different terms 


i - j, . j, <y;» tj 


A A(J) a 

’j S k V 


(27) 


The contribution of the charge exchange reactions £o expression (27) is given 
by the second term. This second term is zer^ for N, nonzero for £ because of 

the nonsymmetric reactions, and nonzero for E because of every charge exchange 
reaction. 


A 

In the case of N, the number of neutrals lost from the j 1 "* 1 species 
packet, Lj is divided between the net ion loading rate, N. , and the production 
of new fast neutrals, , that escape from the system because of 

reactions of type C (Table 3) : 

A 


A 

L, 



a(CE)a 
+ *J n j 


(28) 


as expected on physical grounds. The fraction of L, that contributes to t. 
is time dependent because of the oscillation of the plasma torus about the^ 
satellite plane and other jnore complex motions that are contained in the 
transformation (13). For M and E, the separate j th species contributions are 
coupled as described by equation (27). The overall mass and energy 
conservation is reflected in equation (23), where the first term is the 
contribution to the plasma torus by the neutrals lost from the packets, and 
the second term is the contribution produced by new neutrals being created by 
the packet-plasma reactions of type C in Table 3. 

A 

The cumulative value of the instantaneous rate R supplied to the plasma 
torus in a time period T by all of the neutral species packets, each of which 
has a current flight time t' , may also equally well be formulated. For the 
absolute time interval from t Q -T to t Q , this cumulative value is given by 
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[J R (t\ t -t'|w,Q) dt"] dn dw dt' 


(29) 


where 


0 

t' -T 


if t' < T 1 

if t' > T J 


(30) 


A 

The quantity (29) is useful in that, even though the instantaneous rate R 
integrated over all flight time is time dependent, its cumulative value (29) 


averaged over all flight time t' is not time dependent, ff the source and sink 

processes for the neutral cloud are periodic with the period T. This will be 

discussed further in a later section. It should be noted here in passing, 

however, that in addition to the instantaneous contribution (27) or the 

A A A 

cumulative contributions (29) of the neutral packets to N, M and E, the plasma 
torus may also contribute directly through reactions of type E, F, and G in 
Table 3. These additional processes are not in the realm of the neutral cloud 
description and may be treated independently. 


Physics of an Ensemble Subpacket 

A neutral gas cloud of Io may be described for the secondary atom popu- 
lation by following the spacetime history of all neutral atoms produced by the 
primary atom population. To accomplish this explicitly, the secondary atom 
population will be described in terms of an ensemble subpacket in a way 
analogous to that in which the primary atom population is described in terms 
ensemble packet. The relationship between the packet and the subpacket 
is illustrated in Figure 2. For a neutral gas subpacket of the j*'* 1 species, 
the number of atoms created with ihitial velocity £ to £+d£ along the primary 
packet trajectory r(t' , t^ |w,Q) in the time interval t' to t'+dt' is given by 

A A -► A -► 

Sj(t' >t b |w,fi,£) to d ” dt ' d £ 


( 31 ) 
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where 


A a -+ 

Sj<t',t b |S r n f o 


Pj(t',t b |S,n) gj t' ,t b ) 


(32) 


-* A 

Here, gj (f|w,fi, t' , t b ) is the local initial velocity distribution for all 
secondary atoms created at time t' along this packet trajectory and is 
normalized to unity, 


- A -*• 

J gj(£|w,0,t' ,t b ) d£ - 1 , (33) 

A 

and Pj is given by equation (26) for the case of charge exchange creation of 
fast secondary neutrals. For a secondary population created by elastic ion- 
atom or dissociative electron-molecule collisions, the expression for Pj must 
be appropriately generalized. The creation rate (32) in the subpacket 

formalism occupies a parallel position to the flux (11. b) in the packet 
formalism. 


The location of the center of the subpacket, originally created at 
a flight time t"-=0 with an initial velocity £ at the primary location 
r (t',t b |w,Q) and corresponding to the primary trajectory time t b +t' , will 
be denoted for t">0 by r(t",£; t',t b |w,S). In the plasma coordinate system, 
the location of the subpacket at a flight time t" is therefore given by 


r (t M ,£;t' ,t b |w,Q) - T(t b +t'+t") • r(t" ; t',t b |u,Q) + D(t b +t'+t") (34) 


in an exactly parallel manner to the packet location (13). The velocity 

of the subpacket in the plasma coordinate system will be denoted by 
-► A J 

v (t",£; t',t b |w,0) and is determined by the t" time derivative of expression 
(34). 

The ensemble subpacket contributions for different neutral cloud physical 
quantities of interest can be formulated as in the previous section for the 
packet ensemble. To facilitate this formulation, a double hat "a" notation is 
adopted to indicate that a physical subpacket quantity is described in terms 
of its phase space representation (t",e; t',t b |w,n) along its phase space 
trajectory. The loss collision frequency for the subpacket of the j 1 "* 1 species 
is given by 
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t' ,t b |w,fi) - t' ,t b |w,n), v(t",£; t'.tblw.O))]* 1 (35) 


parallel to the packet expression (14), and the probability that the ensemble 
subpacket of atoms at flight time t ' * exist is defined by 


t',t b |w,fl) = exp [-J i/j (t'",f; t',t b |w,n) dt'"]. 


(36) 


The number of atoms in the subpacket of the j** 1 species at time t" is then 
given by 


ft 


n j ( t ' ' , £ ; t',t b |S, 0 ) dn dw dt' d£ dt" 


(37) 


where 


n-j (t",J; t' ,t b |w,n) « SjCt'.thiw.o,?) j(t",J; f 9 t h \» 9 a). (38) 


In a similar manner to the packet quantities defined in the previous section, 

A -* A 

the subpacket photon emission rate density Ix(t'',£; t # ,t, |w,Q), the loss rate 
A -► A 

^ensity L(t'',f; t ' , t^fu, 0) , and the plasma input rate densities 

R(t'',£; t',tk|w,0) may be defined parallel to the expressions (18), (21), 
and (27), respectively. 


Kinetic Theory Distribution Function 


The kinetic theory distribution function for the neutral gas cloud of the 
j species as noted earlier may be divided into contributions from the 
primary and secondary atom populations 

(primary) ( secondary) 


(V ,X ,t ) 
J o o o 


JOOO JOOO 


(39) 


Here, v Q and x Q are the velocity and spatial coordinates of phase space for an 
observation coordinate frame, and t Q is the absolute time in the observation 
frame. The primary and secondary contributions to equation (39) can be 
directly expressed, respectively, in terms of the primary packet density 
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A A 

n and the secondary subpacket density n defined in the previous two 


subsections : 
(primary) 


f (v x t ) - J d£l / dw J dt' n (t'.t -t |w,0) fi(r ,r.) S(v r ) , (40) 
^ exobase velocity flight J 


A . £ 


space time 

(secondary) 

f.(x o ,v o ,t o ) - J dfl J dw f dt' / d{ /dt" ^j(t" |w,0)5(r o ,rj)6(v o , 

** ° exobase velocity flight velocity flight 

space time space time (4 

A 


In equations (40) and (41), r and r are abbreviated notations for the 
ensemble packet spatial location and velocity 


r - r(t ' ,t 0 -t' |w,0) , r - [— ?(t' ,t b |t5,0)] t t# 
a bo 


(42) 


A A 

A A 

A 


and r and r are abbreviated notations for the ensemble subpacket location and 
velocity: 


k 


r*r(t",£ ;t' ,t Q -t' |w,G) , r [r(t",£ ; t' , t Q - 1 ' |5, O) ] . (43) 


The quantities $(r o ,rj) and 5(v 0 ,rj), respectively, are the spatial volume 
projection operator and the velocity volume projection operator defined as 
follows : 


s( V r> - 


f 1 if r < R < r + Ar 
o o o o 


otherwise 


(44) 


S( v r) - 


r i 

0 


if v < V < v + Av 
o o o o 


otherwise 


(45) 


Here, R q is the location of the ensemble packet (or subpacket) in the 


observation coordinate frame defined and mapped by an appropriately defined 
transformation of the form 


R - T (t ) • R 
o o o 


( 46 ) 


K*1 » 
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4 ^ 

where R * r (or r) and V Q - R Q is the velocity of the ensemble packet (or 
subpacket) in the observation coordinate frame. In addition, the quantities 
A 3r o and A 3 v q are, respectively, the spatial and velocity volume elements 
establishing the volume resolution element in the phase space of the observa- 
tion coordinate system and correspond to the spatial volume Ar Q about r Q and 
the velocity volume Av Q about v Q . 


Spatial Properties of the Neutral Cloud 


The spatial properties of the neutral cloud (e.g., the three-dimensional 
number density or the two-dimensional brightness of the cloud on the sky 
plane) as well as the physical ion input rates to the planetary magnetosphere 
can be determined directly from the kinetic theory distribution function (39) 
or, alternatively, from the expressions derived from the ensemble packet and 
subpacket in the previous sections. In the former case, one multiplies the 
distribution function by a physical quantity of interest, which may in general 
be a function of (r ,v f t ) f and integrates the product over the velocity 
coordinate space v . _ In the latter case, one defines an ensemble average for 
the packet and subpacket physical quantity of interest, producing a result 
identical with that obtained by the former method because of the nature of the 
spatial and velocity volume projection operators (44) and (45)*. The ensemble 
average approach will be explicitly developed below and is the description 
adopted for the numerical evaluation ’of the neutral cloud model discussed in 
the next section. 


A A '\ 

A A • 

If an ensemble packet quantity (e.g., n, Ix,N) and its subpacket quantity 

A A A _ . ■ 

A A • A A 

(e.g., n,I^,N) are denoted respectively by Q(t',t b |w,Q) and 

A A 

t # ,t^|u,Q) f and if the corresponding ensemble -average physical 

quantity in the observer's coordinate frame at r and at an absolute time t 
- o ° 

is denoted by Q(r o ,t Q ) f then 




A 

Q 


A 


> + < Q > 


(47) 


where 

< Q > - J dfl J dw dt' Q(t' ,t 0 -t' |w,ft)6(r 0 ,r) 

exobase velocity flight 
space time 


( 48 ) 
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> 


- J dfl J dw f dt' J d£ f dt" 



k 


k 


Q(t",£; t' ,t Q -t' |w,0)5(r o ,r), 


(49 


where r and r are given by equation (42) and (43) , and where the volume 
projection operator 6(r^, r) is defined by equation (44). To project an 
ensemble average physical quantity Q(r Q , t Q ) onto the observer's viewing 
plane, as is appropriate for the gas density to obtain a column density or for 
the volume photon emission rates to obtain an image of the cloud brightness, 
an integration along the line of sight in the observer's coordinate frame is 
required in addition to those specified in equations (48) and (49). For the 
simple case of no parallax with the z Q -axis along the line of sight, the 
column integrated quantity is given by 


Q(x ,y t ) 
^ o J o o 


r 


Q(r , t ) dz . 
o o o 


(50) 


In this manner, the basic properties of the cloud and also its interactions 
with the magnetosphere may be calculated. 

It should be noted that if the ensemble packet quantity and the packet 
trajectory are periodic with period T in their initial ejection time t^, 
that is , if 


A A A . A 

Q(t',t+T|w,n) = Q(t',t[w,0) 


(51a) 


r(t' ,t+T|w,Q) - r(t' ,t|w,S), 


(51b) 


it then follows that the corresponding ensemble average physical quantities of 
the neutral cloud (47), although time dependent in general, are also periodic 
with period T 

Q(r ,t + T) - Q(r ,t ) . ( 52 ) 

xv o o x o o 

In general, the ensemble -average value of the cumulative value [see eqs. (29) 

A 

and (30)] of any packet quantity Q over a time t Q -T to t Q , 
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t'A 


< [J Q(t",t -t|w,0) dt"] > , 


(53) 


is exactly equal to the ensemble average value of the instantaneous packet 
A 

quantity Q integrated over the time interval 


A 

S < Q > dt; , 


(54) 


to ' 1 


as can be readily verified. A similar result also holds for the subpacket 

quantity. If equations (51) hold, the quantities (53) and (54) are also 

explicitly independent of the time t_, as would, of course, be expected purely 

A A ° 

on physical grounds. For Q - L, this time independence is clearly illustrated 
in the simple case of a constant flux, 


^(w.H,t b ) - Q) , 


(55) 


by the obvious result 


S 


t +T 


< L > dt' 


all 

space 


T f dfl f dw *(w,0) 

exobase velocity 
space 


(56) 


developed directly from expressions (20) and (48). 

The general periodic relationship (52) and the equality of expressions 
(53) and (54) are useful in studying the impact of the neutral clouds on the 
magnetosphere. A period T of approximately 13 hours exists between Io and its 
location in the oscillating/rotating magnetosphere if east-west torus 
differences are assumed small or zero. If, on the other hand, east-west torus 
differences cannot be ignored, a period T of approximately 32 days exists for 
recurrence of the same east-west position and System III magnetic longitude of 
Io. To illustrate the impact of the neutral clouds, the packet expression 
(28) is summed over the various neutral cloud species and its ensemble average 
is taken to obtain 


A ^ N A(CE)A 

<L>-<N> + <Si/| tii >. 

J-l 


(57) 
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This expression divides the instantaneous total packet neutral loss rate 
A 

of the cloud < L > into its instantaneous net ion loading rate for the 
A 

magnetosphere < N > and its instantaneous production rate of fast neutrals 
A f CE) A 

< “j >• This partitioning is time- as well as space - dependent . 

The latter term of equation (57) is, of course, the rate at which subpacket 
fast neutrals are produced. If the packet periodicity (51) holds, the period 
averaged value of each term in equation (57) is independent of time, so that 
the average fractions of neutrals contributing to the net ion loading rate and 
to the production of fast neutrals are given, respectively, by 


o 


N 
C S 

J-l 


(0 a> + 

J 


a(CT). A _ . 

V 1 >v / 


N 

< 2 
j-l 


A A 
v ,n.> 

j j 


(58a) 


and 


N N 

-(fast),-* v . „ A(CE) A , _ A A 

f v '(r ) - < 2 v) 'n> / < S ^.n > 

° j-l J J j-l J J 


(58b) 


where the angular brackets denote the packet ensemble average and the bar 
indicates periodic accumulation, that is. 


A t _+T A 

< Q > - J < Q > dt^ 

t o 


; dt; , 

t O 


(58c) 


and where by construction 


f m a > + f <fast >(? > 

o o 


1. 


(58d) 
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4 . Numerical Approach 

The ensemble average of a packet and subpacket quantity given by 
equations (47) - (49) and its column integrated value given by equation (50) 
may be evaluated numerically. For the packet description, this is accomp- 
lished by dividing the exobase into many small area elements. On each area 
element, the initial velocity distribution as well as the absolute flux are 
specified independently. For each initial condition specified on an area 
element, a trajectory r(t ' , t^ |w,Q) for an ensemble packet is calculated for a 

sufficiently long flight time t 9 to determine accurately the cloud proper- 

A 

ties. Along each trajectory, the packet quantities Q of interest are 

calculated and accumulated in their appropriate spatial volume (or area of the 

sky plane) in such a way as to determine the packet ensemble average (48) or 

its contribution to the column- integrated ensemble average (50). For the 

A 

subpacket description, the source density rate S for the subpacket as defined 

by expression (31) is calculated in the packet description along the packet 

trajectory and has the analogous role that the flux (11a) has in the packet 

A 

description. The source density rate S is numerically divided into many small 
packet flight time increments along the packet trajectory and from each 
increment the subpacket trajectories r(t",£; t',t^|w,Q) are calculated for 
sufficiently long flight time t" to determine accurately the cloud 

k 

properties. Along each trajectory, the subpacket quantities Q of interest are 
calculated and accumulated in the appropriate spatial volume element of the 
observation coordinate frame. 

The packet trajectories r(t' , t^ |t*>,0) and subpacket trajectories 
r (t" t # ,t^|w,0) may be determined by solving the standard (Smyth and 
McElroy 1977) or modified (Smyth 1983) circular-restricted three-body 
equations of motion for each prescribed initial condition. The modified 
equations are required for the sodium and potassium cloud because of the 
additional effects of solar radiation pressure (Smyth 1979, 1983) on these 
atoms. The equations of motion are solved numerically by applying a straight- 
forward Runge Kutta technique. The time steps employed in the numerical 
solution must be selected small enough not only to determine an accurate 
trajectory but also to ensure that the space and time resolution of the 
ensemble average quantities (47) , which are quite sensitive to the oscilla- 
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tional/rotational motion of the plasma torus, are calculated to the level of 
desired accuracy. 

A 

For a particular neutral cloud, calculation of the packet quantity Q, the 
subpacket quantity Q, and their ensemble averages (48) and (49) requires first 
that the relevant reactions and their rates (symbolically indicated by A-D in 
Table 3) and also the excitation rates (18) be identified and defined, and, 
second, that the spacetime properties of the plasma torus be specified along 
the packet and subpacket trajectories. For the sodium cloud, this information 
is presented in the companion paper (Smyth and Combi 1988). 
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5 . Summary 

A brief review of Io's local atmosphere and its extended gas clouds as 
well as their interactions with the planetary magnetosphere was given in 
section 2. The Io neutral clouds were seen to occupy a central position in 
understanding the coupled planet-satellite-magnetosphere system. Analysis of 
Earth-based, rocket, or Earth- orbiting satellite observations of neutral cloud 
emissions provides a direct, very promising, and still largely untapped 
approach for increasing our knowledge and understanding of complex physical 
phenomena in the system. Physically realistic models for the neutral clouds 
are essential in this analysis. These models must include the complex space- 
time interactions that occur between the neutral clouds and the plasma torus 
which have hitherto not been incorporated. The general Io neutral cloud model 
developed here includes these interactions and provides a general and solid 
base for future analysis of cloud data. 

A mathematical description for the general neutral cloud model was pre- 
sented for Io in section 3. For a given neutral gas species, the description 
is divided into a primary cloud population that is initially ejected from the 
satellite exobase and a secondary cloud population that is created in a spa- 
tially extended volume by the primary cloud population through various time- 
dependent processes. For simplicity, the mathematical description for the 
primary and secondary gas populations is explicitly developed for an atomic 
cloud, although it may be readily generalized for a molecular cloud. The 
description of the primary gas cloud may be fundamentally reduced to describ- 
es the time evolution of an ensemble packet of gas along its phase space 
trajectory. The description of the secondary gas cloud may likewise be 
reduced to describing the time evolution of an ensemble subpacket along its 
phase space trajectory. The physics of an ensemble packet and ensemble 
subpacket are presented and are also directly related to the kinetic theory 
distribution function for the gas cloud. The spatial properties of the 
neutral cloud as well as various ion input rates to the planetary magneto- 
sphere are obtained from their ensemble packet and subpacket values by 
evaluating their ensemble averages. A straightforward numerical approach for 
evaluating the model is outlined with appropriate references that define and 
present a method of solution for the equations of motion of the phase space 
trajectories for the ensemble packet and subpacket. 
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As noted earlier, the general model is applied to the Io sodium cloud in 
a companion paper (Smyth and Combi 1988). The application of the general 
model to other neutral gas clouds of Io (e.g., K, 0, S, S0 2 , 
considered in a series of future papers. 


etc.) will be 
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Table 2 

Some Studies for Exospheric Problems in Planetary Science 


Subject 


Reference 


1. Cometary Atmospheres Eddington 1910; Haser 1957, 1966; Keller 

1973; Keller and Thomas 1975; Keller and 
Meier 1976; Combi and Delsemme 1980; Combi 
1980; Festou 1981a, b;. Combi and Smyth 1985; 
Combi, Stewart and Smyth 1986; Combi and 
Smyth 1988a, b 


2. Planetary Atmospheres 

• Mercury Hartle 1971; Hartle, Ogilvie and Wu 1973; 

Hodges 1974; Hartle, Curtis and Thomas 1975; 
Smith et al. 1978; Curtis and Hartle 1978; 
Hodges 1980; Ip 1986; Smyth 1986 

• Venus Penner 1977; Hodges and Tinsley 1981, 1982; 

Rodriguez, Prather and McElroy 1984; Hodges 
and Tinsley 1986 

• Earth Chamberlain 1963; Hartle 1973; Bertaux and 

Blamont 1973; Fahr and Paul 1976; Penner 
1977; Hodges, Rohrbaugh and Tinsley 1981; 
Rodriguez 1983; Bishop 1985 


3. Satellite Atmospheres 
• Moon 


Hodges 1973; Hartle and Thomas 1974; Hodges 
1980 


• Titan Fang, Smyth and McElroy 1976; Smyth and 

McElroy 1977; Smyth 1981; Shemansky, Smith, 
Smyth and Combi 1984; Hilton and Hunten 
1985; Ip 1985; Hilton 1987 

• Io Carlson, Matson and Johnson 1975; Fang, 

Smyth and McElroy 1976; Smyth and McElroy 
1977, 1978; Matson et al. 1978; Carlson et 
al. 1978; Smyth 1979; Goldberg et al. 1980; 
Macy and Trafton 1980; Smyth 1983; Smyth and 
Shemansky 1983; Sieveka 1983; Pilcher et al. 
1984; Smyth and Combi 1983, 1984, 1985; 
Sieveka and Johnson 1985; McGrath and 
Johnson 1987; Smyth and Combi 1988. 
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Figure 1. 


Figure 2. 


FIGURE CAPTIONS 

Motion of a Neutral Cloud Packet. The mathematical description for 
the primary population of the neutral cloud based upon the dynamic 
evolution of neutral cloud packets is depicted. The one neutral 

t 

cloud packet shown on its trajectory at a flight time t > 0 is 
ejected from a satellite exosphere area element dfl with a velocity 
w at flight time t' - 0. 

Motion of a Neutral Cloud Subpacket. The mathematical description 
of the secondary population of the neutral cloud based upon the 
dynamic evolution of neutral cloud subpackets is depicted. The one 
neutral cloud subpacket shown on its trajectory at flight time t" > 
0 is ejected at flight time t" - 0 with initial velocity f from 
the neutral packet. The location of the neutral packet at the 
later time (t" >0) is also indicated. 
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Abstract 


The general model for Io's neutral gas cloud developed in Paper I is 
applied to the sodium cloud. It includes realistic spacetime dependent sinks 
for sodium in the plasma torus, sources for sodium that describe both thermal 
and non thermal processes in the satellite exosphere/corona, and forces on 
sodium atoms exerted by the planetary and satellite masses and by solar radia- 
tion pressure. The sodium cloud model provides a natural and consistent 
physical explanation for the basic morphological character of the cloud and 
its east-west asymmetries as observed from ground-based telescopes. The 
emerging physical explanation is understood in terms of a sodium sink based 
upon Voyager 1 data for a 7° tilted corotating plasma torus with an offset- 
dipole planetary magnetic field in the presence of an east-west electric field 
of -2.8 mV/m, and a nonthermal sodium source peaked at -3 km sec' 1 describing 
an incomplete collisional cascade process driven by the interaction of the co- 
rotating plasma torus with Io. The active exobase area for the nonthermal so- 
dium source is intermediate between a spherically symmetric- source 
26 1 

(~2.7 x 10 atoms sec" 1 ) and a band source (-1.7 x 10 26 atoms sec' 1 ), where 
the latter excludes a portion of the leading and trailing hemispheres. The 
nonthermal source produces a sodium density (in agreement with recent measure - 
ments) of -3.0 x 10 atoms cm' J at the model- adopted exobase of 2600 km radius 
and implies a collisional exobase that is somewhat below that level. At the 
model-adopted exobase, the flux of the nonthermal sodium source dominates by 
roughly an order of magnitude or more the flux of a 1000-2000°K thermal sodium 
source. By itself such a thermal source is incompatible with the sodium cloud 
data. The nature of the nonthermal processes in the exosphere/corona and 
their implications for the other gas clouds of Io are discussed. 
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1. Introduction 

Studies for the four (Na, K, 0, S) known neutral gas clouds of Io coupled 
with observational data provide an excellent remote sensing tool for recover- 
ing physical information regarding the local atmosphere/corona of the satel- 
lite, the interactions of this atmosphere/corona with the corotating plasma 
torus, and the processes that shape the properties of the inner magnetosphere 
of Jupiter. Recognizing this potential, a general model, appropriate to any 
of the gas clouds of Io, was developed by Smyth and Combi (1988, hereafter 
Paper I). This general model is applied here to the sodium cloud with physi- 
cal insights gained that decisively demonstrate the power of this approach. 

The sodium cloud is a particularly good first candidate for this applica- 
tion, for two major reasons. First of all, the sodium cloud is significantly 
2 3 

brighter (—10 to 10 ) than the other neutral clouds and has been observed 
extensively during most of the last 16 years, with increasingly btf^ker image 
data obtained in the last decade. This not only allows typical characteris- 
tics of the cloud to be identified unambiguously, but also provides a long- 
term data base that, upon understanding its physical basis, gives a direct and 
unique means of studying the time history of the Io-planetary-magnetosphere 
system since about 1972. Second, the sodium cloud has interactions with the 
plasma torus that are less complex than those for the atomic potassium, 
oxygen, and sulfur clouds. The excitation mechanism for the sodium cloud (and 
the potassium cloud) is solar resonance scattering, which is independent of 
the plasma torus properties, in contrast to the primary electron impact exci- 
tation mechanism for the oxygen and sulfur clouds. The sodium cloud, further- 
more, has a very small cross section for charge exchange reactions with the 
major oxygen and sulfur ions of the torus but a very short lifetime that de- 
pends on the plasma torus electrons through electron impact ionization. Simi- 
lar results hold for the dimmer Io potassium cloud, although charge exchange 
cross sections are considerably larger. In comparison, elastic collisions, 
charge exchange reactions, and electron impact ionization are all important in 
determining the lifetime of the oxygen and sulfur clouds. Because of these 
properties, the study of the sodium cloud allows the different plasma-neutral 
gas processes to be more easily separated and thereby allows a greater level 
of confidence and certainty to be gained in the analysis of all the observa- 
tional cloud data for Io. 
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The studies of the sodium cloud in this paper are restricted to ex- 
plaining some of the typical or most basic characteristics of the cloud drawn 
from published and unpublished observational data. The purpose of this is to 
verify the analytic usefulness of the model and to provide an overall consist- 
ent and physical understanding within which the answers to more refined ques- 
tions may be pursued in future studies of specific observational data. In 
keeping with this purpose, for example, three successively more complex des- 
criptions of the plasma torus were considered and were sufficient to explain 
the basic characteristic of the cloud. A yet more complex description for a 
plasma torus with properties inherently dependent upon the System III magnetic 
longitude angle was not included, although observational data for the plasma 
torus suggest that such features are present. Such System III asymmetries in 
the plasma torus would produce some time -dependent perturbations in the basic 
morphology of the sodium cloud. The study of such modifications will be pur- 
sued in future model investigations of specific sodium cloud observations. 

This paper is organized into five remaining sections. In section 2, past 
models for the sodium cloud are reviewed and the motivation for developing 
better models is discussed. The general model of Paper I is adapted to the 
sodium cloud in section 3 , and the interactions of sodium in the plasma torus 
are discussed. In section 4, the model is used to investigate and provide a 
consistent physical explanation both for the most basic morphological charac- 
teristics of the spatially extended sodium cloud and for the east-west inten- 
sity asymmetry observed very near Io. This leads to definite conclusions re- 
garding the flux distribution and density of sodium in the exosphere of 
1 ° > level of the collisional exobase in the atmosphere, the nature of 

the interaction of the corotating plasma torus and the satellite atmosphere/ 
corona, and the presence of an east-west: electric field in the planetary 
system. These results and some of their implications are discussed in 
section 5. The paper is briefly summarized in section 6. The Appendices con- 
tain mathematical expressions to define the plasma torus descriptions de- 
veloped in section 3 and the sodium atom flux distributions adopted in 
section 4. 
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2. Review of Past Models 

Models for the lo sodium cloud have evolved over the past decade by 
gradual refinements in the description of the source, lifetime, and orbital 
dynamics of sodium atoms in the Jovian environment as summarized in Table 1 . 
These refinements have been introduced in response to improved cloud data and 
also to new data obtained for Io and the plasma torus. All of the cloud 
models in Table 1 assume that sodium atoms are initially ejected radially 
(with specified source conditions) from Io's surface or local atmosphere (non- 
radial ejection is also considered by Pilcher et al. 1984) and then travel un- 
hindered along trajectories in the circumplanetary space until they are lost 
from the cloud by some ionizing sink mechanism. The density and brightness of 
the cloud are determined by appropriately following many sodium atoms along 
their trajectories, with proper account taken of ionization and excitation 
phenomena. For most models, these trajectories are determined by solving in 
three dimensions the circular-restricted three-body equations of motion for 
the Jupiter- Io- atom system. 

The results of early modeling of the sodium cloud were to identify 
electron impact ionization as the dominant sink mechanism (Carlson, Matson, 
and Johnson 1975) , and to explore the implication of this finite ionization 
lifetime (assumed to be spatially uniform) on the cloud morphology (Carlson, 
Matson, and Johnson 1975; Fang, Smyth, and McElroy 1976; Smyth and McElroy 

1977) , based upon limited slit-average data showing the absence of a symmetric 
cloud about Jupiter. Upon the acquisition of two-dimensional image data for 
the sodium cloud in 1976 and 1977 (Murcray 1978; Matson et al. 1978; Murcray 
and Goody 1978) , the character of a bright and predominantly forward cloud of 
sodium (i.e., a cloud leading Io in its orbit) became evident. In modeling 
efforts (Matson et al. 1978; Smyth and McElroy 1978), the presence of the pre- 
dominantly forward cloud was reproduced by restricting the surface from which 
atoms were ejected from Io to the inner hemisphere (or within ±30° of it) and 
by limiting the forward-cloud length by selecting a suitably short spatially 
uniform sodium lifetime (-20-30 hr). The forward cloud was shown to depend 
primarily on the low speed components (-3 km sec"^) of the initial velocity 
distribution with a source rate of about 2 x 10 sec (Smyth and McElroy 

1978) , while a higher velocity tail (<15 km sec“^) was also shown to be 
important for populating the cloud north and south of Io (Trafton and Macy 
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1975; Macy and Trafton 1980; Trafton 1980b) and in reproducing the asymmetric 
character of observed sodium line profiles (Smyth and McElroy 1977; Carlson et 
al. 1978; Macy and Trafton 1980). East-west shape asymmetries, first observed 
in the forward cloud by Goldberg et al. (1978) when images of diametrically 
opposite phase angles were compared, were shown to arise naturally upon incor- 
porating in the orbital dynamics the acceleration of sodium atoms by solar ra- 
diation pressure (Smyth 1979, 1983). East-west intensity asymmetries observed 
very near Io (Bergstralh et al. 1975, 1977) were also shown (Smyth 1983) to be 
caused at least in part by the action of solar radiation pressure on sodium 
atoms ejected at slow velocities (<2.2 km sec~^) from Io. More recently, 
atmospheric shielding of the sunlit surface by S (>2 vapor to vary the effi- 
ciency at which sodium may be sputtered from the surface by corotating ions as 
a function of the satellite geocentric phase angle has been suggested (Thomas 
1986) as an alternative mechanism for the east-west intensity asymmetry. 
Asymmetries in the plasma torus that may more completely explain this inten- 
sity asymmetry will be discussed in sections 3 and 4. 

Until 1979, a fundamental deficiency in models for the sodium cloud was 
the lack of any spatial description of the magnetospheric plasma and hence the 
ionization lifetime of sodium in the circumplanetary space of Jupiter. The 
encounter of Voyager 1 with Jupiter provided this missing information and in- 
dicated the presence of an extremely radially asymmetric plasma torus, spa- 
tially including Io's orbit and oscillating about its orbital plane, that 
yielded a sodium lifetime as small as 1-2 hours near Io. Incorporating a ra- 
dial, one -dimensional (i.e. f vertically time -averaged) spatial description of 
the sodium lifetime based on early Voyager data to model a forward cloud image 
obtained on the Voyager encounter day by ground-based observations, Goldberg 
et al . (1980) obtained a reasonable match for a sodium source rate of about 
3 x 10 atoms sec by restricting gas ejection to one hemisphere. This 
source rate is about an order of magnitude larger than that determined earlier 
by Smyth and McElroy (1978) and reflects the increased escape of sodium re- 
quired to supply the forward cloud when sodium ionization by the plasma torus 
is explicitly included in the model. Incorporating a more realistic two- 
dimensional oscillating description of the plasma, Smyth and Combi (1983) pre- 
sented model calculations for the forward sodium cloud suggesting an isotropic 
source rate of 2 x 10 atoms sec and demonstrating, as suggested earlier by 
Smyth (1979, 1983), that the predominantly forward cloud was actually caused 
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by the extremely radial asymmetric lifetime sink of sodium in the plasma torus 
and not necessarily by a hemispherical sodium source at Io. This subject will 
be explored further in section 4. The sodium cloud model of Smyth and Combi 
(1983) was later used (Pilcher et al. 1984) to explain the time -dependent 
north/south oscillatory nature of less -bright sodium directional features in 
terms of the oscillating sodium sink and a high-velocity (-20 km sec"'*’) sodium 
source of about 1 x 10 atoms sec’ produced by magnetospheric wind-driven 
escape, preferentially from the equatorial region of Io. The relationship of 
such a sodium source to the longitudinal and latitudinal distribution of dif- 
ferent volcanic plumes (McEwen and Soderblom 1983) , thermal emission from hot 
spots (Johnson et al. 1984; McEwen et al. 1985), and surface composition 
(Nelson et al . 1980) poses interesting questions which are now beginning to be 
probed by the neutral cloud models. 

The advances anticipated in physical insight by accurately describing the 
spacetime interactions of the plasma torus with Io and its extended neutral 
cloud are aptly illustrated by the sodium cloud model of Smyth and Combi 
(1983). In general, electron impact ionization, ion-neutral charge exchange 
reactions, and even elastic ion-neutral collisions (Brown, Pilcher, and 
Strobel 1983) between the corotating magnetospheric plasma and lo's local and 
extended atmospheres are processes that must be evaluated and included where 
necessary in developing each new cloud model. For sodium, improved versions 
of the model of Smyth and Combi (1983) are developed in section 3 which in- 
clude even more complete descriptions for the plasma torus as illustrated in 
section 4. Applications of these more complete models to the large JPL Io 
sodium cloud data base recently reviewed by Goldberg, Garneau, and LaVoie 
(1984) and to the complementary and novel Io sodium cloud data-cube images 
(Doppler resolved, two-dimensional sky plane images) recently reported by 
Trauger (1984, 1985) provide excellent opportunities for answering fundamental 
questions in future studies. 
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3. Model Description 

The general mathematical formalism developed in Paper I to describe the 
neutral gas clouds of Io is here explicitly applied to the Io sodium cloud. 

For a given neutral gas, the formalism divides the cloud into a description 
for a primary population and a description for a secondary population. The 
description for the primary population is based upon following the spacetime 
history of an ensemble packet of atoms initially ejected from the satellite 
exobase. The description for the secondary population is similarly based upon 
following the spacetime history of ensemble subpackets that are created in the 
spatially extended circumplanetary volume about Jupiter from primary cloud 
atoms or molecules. The spacetime history of an ensemble packet (or sub- 
packet) and the creation of ensemble subpackets depend upon the physical pro- 
cesses by which the neutral gas interacts with the plasma torus. 

a • Interactions in the Plasma Torus 

The interactions of sodium atoms with the plasma torus are sum- 
marized in terms of chemical reactions and their minimum reaction times in 
Table 2. There are four basic processes: ionization by solar photons and 

electrons (R1 and R2) , charge exchange (R3 - R5) , charge transfer (R6 and R7) , 
and elastic ion-neutral collisions (R8). 

Ionization is the most important process for removing sodium atoms di- 
rectly from an ensemble packet or subpacket. Electron impact ionization (R2) 
has the smallest reaction time in Table 2 by more than an order of magnitude 
and effectively determines the lifetime of the cloud. This lifetime for so- 
dium in the plasma torus, however, varies significantly both in space and in 
time. For an assumed longitudinally symmetric plasma torus tilted by 7 
degrees relative to Io's orbital plane, the electron impact ionization life- 
time of sodium is shown in Figure 1 for plasma properties to be described 
later that are appropriate to the Voyager 1 encounter epoch. A minimum life- 
time of between 1 and 2 hours occurs on the centrifugal equator between 5.7 
and 6.1 planetary radii (Rj) from Jupiter (i.e., ~± 0.2 Rj from the place 
where Io's orbit intersects the centrifugal equator plane at -110° and 290° 
System III magnetic longitude). Inside of a radial distance of about 5.4 Rj 
on the plasma centrifugal symmetry plane, the sodium lifetime increases ex- 
tremely rapidly because of a very steep decrease primarily in the electron 
temperature . 
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Charge exchange reactions (R3-R5) and charge transfer reactions (R6 and 
R7), in contrast to ionization processes, remove only a small amount of sodium 
from the cloud. Charge exchange reactions remove sodium from the primary 
cloud (an ensemble packet) by creating ensemble subpackets of fast neutrals 
(~60 km sec”^) which effectively escape from the Jupiter system unless they 
are trapped in the outer magnetosphere. Because of extremely small cross sec- 
tions for 0 + and S + with Na, only fast sodium atoms are effectively pro- 
duced. Charge transfer reactions remove sodium directly from the primary 
cloud by creating additional ions. Both charge exchange and charge transfer 
reactions have their minimum lifetimes between about 5.0 Rj and 6.0 Rj on the 
centrifugal equator where the corotating ions have their maximum densities. 
Inside a radial distance of about 5.4, where the electron impact ionization 
lifetime increases rapidly, charge transfer reactions are, however, still not 
effective, since the dominant ion species are then S + and 0 + . In this inner 
region, charge exchange reactions with Na + are also not effective, since the 
ion densities are contained by the decreasing ion temperature to be very 
closely (~± 0.1 Rj) confined to the centrifugal equator, thereby minimizing 
their spatial overlap with the neutral cloud during the nominal 15-20 hours of 
flight time of sodium atoms in the forward cloud (see Figure 2). Consequent- 
ly, electron impact ionization is the only lifetime process that will be con- 
sidered in this paper for calculating the loss of sodium atoms from an en- 
semble packet. Furthermore, because the ensemble subpackets created by reso- 
nance charge exchange (R5) have a relatively low production rate and move so 
rapidly in the circumplanetary space that they contribute only minimally to 
the bright forward (Region B) and near- trailing lo (close Region C) portion of 
the cloud, they may, for purposes of this paper, also be safely neglected. 
Those charge exchange reactions that occur directly with the local atmosphere 
within the Lagrange sphere (~ 6 satellite radii in radius) about Io and that 
are important sources for the fast component of the fainter sodium observed 
remote from Io by Brown and Schneider (1981) and more recently near Io by 
Schneider, Hunten, and Brown (1984) and Trauger (1984, 1985) may be properly 
included in the source characteristics of the ensemble packets. 

The fourth process in Table 2, ion-neutral elastic collisions (R8) , as 
noted by Brown, Pilcher, and Strobel (1983), may be either "hard" (resulting 
from the repulsive core of the interatomic potential) or "soft" (resulting 
from the attractive force induced by the passing ionic charge). The "hard" 
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c °H^ s i° ns result in significant momentum transfer to the sodium atoms and 
cause a simple but small loss of material from the primary cloud that creates 
ensemble subpackets and are, in this way, somewhat similar in effect to the 
resonance charge exchange reaction (R5) . "Soft" collisions will result in 
much smaller momentum transfer that severs the kinematic connection between 
the initial escape conditions at Io and the morphological structure of the 
near (Region B, close Region C) cloud and thereby also create ensemble sub- 
packets . The attractive potential for the soft sodium collisions has recently 
been shown (Johnson 1985) to be much smaller than earlier estimates (Brown, 
Pilcher, and Strobel 1983) so that only hard collisions need be considered for 
sodium. Since the "hard" collision reaction time for R8 listed in Table 2 is 
similar in magnitude and spatial character to that of reaction R5 and is more 
than an order of magnitude larger than reaction R2 , the creation of ensemble 
subpackets by reaction R8 can also be safely neglected for purposes of this 
paper. The explicit development of the ensemble subpacket formalism for the 
secondary population of sodium atoms will therefore not be considered further 
in this paper. 

b. Primary Cloud Description 

The description of the primary population of sodium atoms is based 
upon following in flight time t # an ensemble packet of atoms ejected initially 
(t'-*0) at an absolute time t^ with a velocity w from a surface area element dfl 
centered upon the angular direction 0 of the satellite exobase (see Figure 1 
in Paper I). The location of the center of the packet at time t'>0 in a coor- 
dinate system fixed relative to the satellite is denoted by r(t' , t b |w,n) . The 
packet trajectories r(t' , t^ |w,Q) for sodium atoms are computed by numerically 
solving the modified circular-restricted three-body equations of motion (see 
Smyth 1983). These equations of motion include both the gravitational field 
Jupiter and Io and the effects of solar radiation acceleration. 

To describe the interactions of the ensemble packet with the plasma 
torus, a transformation at an absolute^ time t of a location R in the satellite 
fixed coordinate system to a location R in the plasma torus coordinate system 
is required and denoted by 

•H- 


R - T(t) • R + D(t) . 


( 1 ) 
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This transformation is explicitly defined, for purposes of this paper, in 

44 - 

Appendix A. The tensor part of the transformation, T(t), explicitly includes 
the tilt angle of the plasma torus and the absolute time dependent angular mo- 
tion of both the plasma torus and Io about Jupiter. The vector part of the 
transformation D(t) explicitly includes the effects of the magnetic dipole 
offset distance and also the stretched-offset distance of the plasma torus be- 
cause of the presence of an east-west field in the Jupiter system. The loca- 
tion r(t',t b |w,Q) in the plasma torus coordinate frame of the packet at flight 
time t' is therefore given by 

?( t\t b |w f S) «T(t b +t') • r(t',t b |w,0) +D(t b +t'). (2) 


The probability that the ensemble packet of atoms at flight time t ' 
exists is 


P(t',t b |w,0) - exp [-J i/(t" ,t b |w,n) dt"] 


(3) 


A A 

where v(t' , t b | td,fl) is the packet collision frequency which depends upon the 
plasma torus properties and is divided into a contribution from ionization 
(I), charge exchange (CE), and charge transfer (CT) , 


A 

v 


A( I ) A A(CE) A A(CT) 

i/ + i/ +i/ 


(4) 


as discussed in Paper I. As was seen in the above discussion of sodium atom 
interactions with the plasma torus, 0^^ for electron impact ionization is the 
only term that needs to be retained in this paper in evaluating expression 
(4). Three different descriptions for varying in spatial complexity will 

be developed in the next subsection. 

The spatial properties of the neutral cloud are obtained by properly in- 
tegrating the contributions of an ensemble packet over its complete flight 
time and over all initial and exobase conditions. Expressions for the contri- 
butions of the packet to various cloud properties (e.g., number density, 

brightness, ion input rate to the plasma torus) are defined in Paper I. In 

A A 

general, if such a packet contribution is denoted by Q( t ' , t b | w,fl) , the corres- 
ponding spatial property of the neutral cloud in the observer's coordinate 
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-t ^ _► 

frame at an absolute time t Q is denoted by Q(r Q ,t ) and is given by the 

ensemble average < Q >. The ensemble average < Q > is defined by equation 
(48) in Paper I and the transformation T Q (t) required in this definition is 
given, for purposes of this paper, in Appendix B. 

For model calculations presented in the next section, the only sodium 
cloud property to be calculated is its D-line emission brightness on the sky 
plane as viewed from Earth. The excitation mechanism for the sodium cloud 
emission in the D x line ( X 1 * 5895.92 A) and D 2 line (A 2 - 5889.95 A) is solar 
resonance scattering. The rate at which a sodium atom absorbs (and then 
instantaneously emits) photons in the D 1 line (i-1) and D 2 line (i*=2) is given 
by 


2 (nF v > 

\ - V” < rz > -isr V 2 ' 

l e i 


(5) 


where R is the Sun-Jupiter distance in astronomical units, is the oscil- 

lator strength of the i line, nF is the solar continuum- level photon flux 

i 2 

between the and D 2 lines at R-l, e^/m e c is the classical electron radius, 


(ttF ) = 6.80 x 10' 5 — F. 

V . c A ' 


A 2 


( 6 ) 


(see Allen 1976), and 7 ^(v) is the fraction of this solar continuum flux 
available to the sodium atom as it is Doppler shifted out of the bottom of 
the solar D-line Fraunhofer absorption feature by the instantaneous radial 
velocity v of the sodium atom relative to the Sun. Further discussion of 
the emission rate (5) is given by Smyth (1979, 1983). The total photon emis- 
sion rate for the packet is given by using equation (5) in equation (19) of 
Paper I. The D-line emission brightness on the sky plane is then determined 
assuming no parallax by the column projection (equation 50 in Paper I) of 
the ensemble average (equation 48 in Paper I) of this emission rate for the 
packet. 

c . Sodium Lifetime Description 

The ensemble packet collision frequency (4) and lifetime r for 
sodium are effectively determined by electron impact ionization and are 
therefore given simply by 

A A ( 1 ) - - 

v = i/ - n K(T ) - r 
e e 


( 7 ) 
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Here n e and are the electron density and temperature in the plasma coor- 
dinate frame evaluated at the physical location of the ensemble packet, and 
K(T^) is the electron impact ionization rate coefficient shown in Figure 3 as 
a function of electron temperature. For model results to be presented in the 
next section, it is sufficient to describe the electron density n e and temper- 
ature for three different steady- state physical situations: (1) Case A, a 

rotating, Jupiter-centered plasma torus tilted 7° with respect to the orbit 
plane of Io, (2) Case B, a rotating, offset-dipole-centered plasma torus 
tilted 7° with respct to the orbit plane of Io, and (3) Case C, a rotating, 
offset-dipole-centered plasma torus tilted 7° with respect to the satellite 
orbit plane in the presence of an east-west electric field. 

The electron density n e and temperature for Case A and Case B may be 
expressed as a function of only the position vector, R, in the plasma torus 
coordinate frame. The collision frequency u and lifetime r for sodium in the 
plasma torus coordinate frame are then given by 

KR) - n (R) K(T (R)) - (R) (8) 

e e 

so that equation (7) becomes 

u - r' 1 - u(i). (9) 

— ► 

Here r is* defined by equation (2), where D(t + t') is zero for Case A and 

b *► 

nonzero for Case B, and the transformation tensor T includes both the seven- 
degree tilt and the relative rotational motion between the satellite and 
plasma coordinate frames. In particular, for purposes of this paper, the 
electron description in Case A and Case B will be further limited to Voyager 1 
encounter conditions and the plasma torus will be assumed longitudinally sym- 
metric. The resulting two-dimensional lifetime r (/?,£) is given in Figure 1, 
where p is the radial along, and z is normal to the centrifugal plasma 
plane. The radial coordinate p in Figure 1 is centered on Jupiter in Case A 
and on the dipole-offset location in Case B. The lifetime of sodium at Io's 
location for both Case A and Case B is shown in Figure 4 as a function of the 
System III longitude angle of Io and are seen to be similar in magnitude and 
to behave approximately like a phase- shifted sine curve with minimum values 
occurring when the satellite intersects the centrifugal plasma plane. 


i 
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The lifetime for sodium in Figure 1 is determined as follows. The 
adopted values for the electron temperature and density along the centrifugal 
equator are shown respectively in Figures 5 and 6 and reflect the currently 
best available values. Assuming a longitudinally symmetric plasma torus, a 
constant value along dipole field lines for the electron temperature and ion 
temperature (also shown in Figure 5) , and relative ion abundances given by un- 
published UVS and PLS results (Shemansky 1985; Bagenal 1985a), the electron 
density n e (p, z) in Figure 7 is determined by iteratively solving the diffusive 
equilibrium equations for a multispecies plasma (see Bagenal 1985b). The so- 
dium lifetime in Figure 1 then follows directly by using in equation (8) the 
electron density of Figure 7 and the electron temperature of Figure 5, as- 
suming T^ is constant along dipole field lines. 

A more complex description of the lifetime for the sodium cloud is ob- 
tained in Case C, where both the radial offset of the dipole magnetic field 
from the center of Jupiter and the east-west stretched shift of the plasma 
torus from the center of the magnetic field caused by an east-west electric 
field are included. The magnetic dipole offset, which is explicitly included 
in the packet formalism by the transformation vector D(t) in equation (1), 
makes the lifetime in a Jupiter-centered coordinate frame vary in System III 
magnetic longitude angle, even if the torus were assumed longitudinally sym- 
metric in the plasma torus coordinate frame. The east-west shift has three 
distinct effects on the plasma properties: (1) a geometric shift of the 

plasma torus center toward the east, (2) a decrease in the electron tempera- 
ture T^ of a plasma volume element as it moves from west to east of Jupiter, 
and (3) a decrease in the electron number density n e of a plasma volume ele- 
ment as it moves from west to east of Jupiter. The effect of the geometric 

shift is explicitly included in the packet formalism by the transformation 
— ► 

vector D(t) in (1). The second and third effects, however, make the electron 
density n e and temperature T^ depend, in addition to their spatial loca- 
tion R, on their east-west angular orientation (i.e., its "phase angle" $) 
about Jupiter. The collision frequency u and lifetime r for sodium in the 
plasma torus coordinate frame therefore have the functional form 

- - ~ ~ - ? ._i ? 

i/(R,$) - n (R.,$) K(T (R,$)) « r (R,$) (10) 

C 0 


so that equation (7) becomes 
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A -1 - - 

V - T - »/(r,$(t b +t' )) , (11) 

-► 

where r is given by equation (2) and $ is evaluated at the absolute time t^+t' 

corresponding to the position of the ensemble packet at flight time t' . 

For purposes of this paper, the electron density n e and temperature 

shall furthermore be assumed for Case C to be intrinsically symmetric in the 

plasma torus longitudinal coordinate for any fixed value of The life- 
time r in equation (10) will therefore depend only upon three independent va- 
riables, i.e., r(p, z,$). An explicit function for r(p, z,$) is developed in 
Appendix C. It is given by equation C.9 and is based upon assuming that n e 
and used in Figure 1 are valid at western elongation (i.e., for $ - 270°) 
near where the Voyager 1 PLS data were acquired. The lifetime for sodium in 
the plasma torus coordinate frame as determined by equation C.9 is illustrated 
in Figure 8. The motion of Io in the plasma torus is also shown in Figure 8 
as determined by the transformation (1). The complete radial domain of Io in 
the plasma torus coordinate frame is shown in Figure 9 as a function of System 
III magnetic longitude. The corresponding lifetime of sodium at Io's location 
in the plasma torus is also shown in Figure 10 as a function of Io System III 
magnetic longitude and Io phase angle. Note that the lifetime of sodium for 
Io at eastern elongation (90° Io phase angle) is significantly longer than at 
western elongation (270° Io phase angle). This will result in a brighter 
sodium cloud near the satellite when Io is east of Jupiter, as discussed in 
* the next section. 
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4. Model Results 

The new sodium cloud model developed in the previous section will now be 
used to explore two different topics: (1) the physical significances of a 

predominantly forward sodium cloud and its east-west orbital asymmetries (the 
latter discovered by Goldberg et al. 1978, 1980) and (2) the role of the 
plasma torus in explaining the east-west intensity asymmetry of the near Io 
cloud discovered by Bergstralh et al. (1975, 1977) and its extension into the 
forward cloud (more recently reported by Goldberg, Garneau, and. LaVoie 1984). 

For all model results to be presented, the sodium D 2 ’line emission 
brightness on the sky plane is computed for a satellite exobase of 2600 km in 
radius that is divided into 326 equal-area elements. An initially radial 
velocity vector is assumed for the trajectory of each ensemble packet ejected 
from each active area element. Active area elements are chosen to simulate a 
uniform flux of sodium ejected from the entire exobase or from a portion of 
the exobase such as a quadrant * hemisphere or band region. In computing the 
ensemble average for the sky plane D 2 brightness of the predominantly forward 
sodium cloud, an initial mean speed of 2.6 km sec"'** (i.e., a monoenergetic 
source, typical of the speed of the peak of the expected velocity distribution . 
for the forward cloud) is chosen for each active area element. For model cal- 
culations examining the east-west intensity asymmetry of the sodium cloud, 

initial speed distributions are in addition assumed for each active area ele- 

A 

ment. The ensemble packet trajectories r ( t ' , t^ |t*>,0) originating from each 
active area element are calculated including the effects of solar radiation 
pressure and the gravitational attraction of Io and Jupiter as discussed 
earlier. 

a. Forward Sodium Cloud 

The presence of a bright, well developed forward sodium cloud and 
the absence of an equally bright trailing cloud have been observed for Io for 
over a decade (Trafton, Parkinson and Macy 1974; Macy and Trafton 1975a, b; 

Munch and Bergstralh 1977; Trafton and Macy 1978; Murcray and Goody 1978; 
Murcray 1978; Matson et al. 1978; Goldberg et al. 1980; Morgan, Pilcher, and 
Fertel 1983; Trauger 1983, 1984; Pilcher et al. 1984; Morgan and Pilcher 1984; 
Goldberg, Garneau and LaVoie 1984; Trauger 1985). In the absence of knowledge 
about the spatial structure of the extremely radially asymmetric sodium sink 
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in Figure 2, pre-Voyager attempts to model the cloud interpreted the forward 
cloud by limiting a primarily low-velocity sodJtun source to a hemispherical 
region about the satellite and assuming a spatially uniform lifetime (Carlson, 
Matson and Johnson 1975; Smyth and McElroy 1977, 1978; Matson et al. 1978; 
Smyth 1979; Macy and Trafton 1980). The hemispherical source region, centered 
more or less on the inner satellite hemisphere, supplied sodium preferentially 
(because of orbital dynamics and radially ejected atoms) to the forward cloud. 
The adopted spatially uniform lifetime of about 20 hours both fixed the length 
of the forward cloud to that roughly observed and spatially limited sodium 
atoms that would have eventually (in 25-50 hours) overpopulated and over- 
extended the observed trailing cloud. 

As discussed in section 2, the first post-Voyager modeling of the cloud 
(Goldberg et al. 1980) incorporated the time-averaged radial structure? of the 
oscillating plasma torus to describe the sodium lifetime in a simple one- 
dimensional fashion. With this model, the same conclusion was drawn, namely, 
that the sodium source must still be hemispherical to reproduce properly the 
predominantly forward cloud. More recently, Smyth (1983) suggested that a 
real-time oscillating torus may enhance the radial asymmetry of the sodium 
lifetime so that the forward cloud might be consistently modeled with an iso- 
tropic sodium source. If this is true, as it shall be shown to be below, the 
predominantly forward cloud would then reflect the presence of the plasma 
torus, not the angular character of the sodium source. The new sodium cloud 
model will be used to explore and further develop this idea. 

To understand the impact of the plasma torus on the morphology of the 
sodium cloud, a number of monoenergetic model calculations are performed 
for a variety of active sources on the satellite exobase and for the three 
plasma torus descriptions (Cases A, B, C) discussed in section 3. Low-speed 
(—3 km sec"^) monoenergetic model calculations simulate reasonably well the 
main features of the forward cloud. This is possible, since earlier studies 
of spectral and spatial data have indicated that the most probable speed of 
the velocity dispersion at the exobase occurs at this low- speed value and that 
a significant percentage (>30-50%) of the velocity dispersion occurs in the 
speed interval (-2-5 km sec‘^) where atoms can dynamically populate the for- 
ward cloud. The velocity dispersion of sodium atoms at the exobase will be 
discussed in the next subsection. 
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Simple comparisons of monoenergetic model calculations for the forward 
cloud and observed D 2 -line emission images (Murcray 1978; Murcray and Goody 
1978; Matson et al. 1978; Goldberg, Garneau, and LaVoie 1984; Goldberg 1987) 
suggest a low-velocity (-2.6 km sec satellite source rate of -0.5-1 x 10 2 ** 
atoms sec' 1 . This source rate is smaller than the -2 x 10 26 atoms sec' 1 value 
for the forward cloud noted by Pilcher et al. (1984), because a lower and more 
accurate electron temperature is adopted here that increases the sodium life- 
time near lo by about a factor of 2. This forward cloud source rate of 
Pilcher et al. (1984) is similar to the —3 x 10 2 ** atom sec" 1 source rate 
determined in the earlier modeling by Goldberg et al. (1980) assuming a 
sputtering velocity dispersion and a simple one -dimensional description of the 

plasma torus sink, particularly when the higher velocity (-20 km sec' 1 ) source 
26 1 

of -1 x 10 atoms sec' 1 discussed by Pilcher et al. (1984) is included to 
simulate properly the fainter sodium directional features observed in the 
trailing cloud. Considering the larger velocity dispersion, a forward cloud 
source rate of -1 x 10 2 ** atoms sec' 1 is also consistent with the sodium D 2 
brightness of many tens of kiloRayleighs observed by Bergstralh et al. (1975, 
1977) in their 3x8 arc second slit centered on Io. This will be discussed 
in the next subsection. 

A forward cloud source rate of -1 x 10 26 atoms sec' 1 is, however, much 
larger than the earlier model calculations assuming a simple cutoff lifetime 
and a hemispherical source of only -2 x 10 25 atoms sec' 1 (Smyth and McElroy 
1978; Smyth 1983) and reflects the additional flux required to balance the 
plasma torus sink of Figure 1. For most forward- cloud model results to be 
presented, then, monoenergetic (2.6 km sec 1 ) ejection from the active exobase 
area with a flux of 2.3 x 10^ atoms cm 2 sec 1 (referenced to the surface area 
of Io and corresponding to a spherically symmetric source rate of -1 x 10 2 ** 
atoms sec 1 ) is therefore adopted as a standard source condition. An Io 
System III magnetic, longitude angle of 140 degrees is also adopted for con- 
venience to set the position for the plasma torus, since this is where the 
largest model-calculated value of the east-west intensity ratio is obtained, 
as discussed in the next subsection. Other values of the System III magnetic 
longitude angle of Io assumed in the model will modify the brightness of the 
sodium emission somewhat but will not alter the basic spatial shape of the 
predominantly forward cloud under examination here. 
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Monoenergetic calculations for the D 2 emission brightness of the forward 
cloud with Io at eastern elongation are presented in Figure 11. The standard 
exobase source conditions noted above are assumed for five different active 
exobase areas: symmetric ejection and ejection from each of its four separate 

quadrants. The simplest plasma torus description of Case A (i.e. , a rotating 
centrifugal plasma symmetry plane inclined 7 Q with respect to the satellite 
plane and centered on Jupiter) is also assumed in Figure 11. For symmetric 
ejection, the absence of an elongated trailing cloud and the presence and 
length of the predominantly forward cloud occur naturally in Figure 11 because 
of the spatially nonuniform nature of the sodium sink (see Figure 2) as anti- 
cipated by Smyth (1983). Comparisons for the four quadrant and symmetric 
sources indicate, as expected on dynamical grounds (Smyth and McElroy 1978), 
that the trailing inner quadrant is critical in populating the forward cloud 
and that the outer leading quadrant is most important in populating the 
trailing cloud. Model calculations for hemispherical sources are shown in 
Figure 12 and are obtained by properly adding model results for quadrant 
sources in Figure 11. Both trailing and inner hemispherical sources produce 
a well developed forward cloud as expected from the quadrant results of Figure 
11, but the trailing cloud for the inner hemisphere source is underdeveloped 
and somewhat less extended than typically observed in ground-based images for 
Io at this phase angle. 

The model results of Figures 11 and 12 suggest that the spatial morphol- 
ogy of the bright sodium cloud may be consistent with a symmetric source or 
other sources that appropriately mix sodium ejection from a portion of both 
the inner and outer exobase hemispheres. One such alternative source, con- 
sidered below, is a band source. Such an exobase source for high-velocity 
sodium (-20 km sec*'*') was suggested by Pilcher et al. (1984) to explain the 
fainter sodium directional features. The active exobase area for the low- 
velocity (2.6 km sec"'**) band source adopted here is -±40 degrees in width and 
centered on the great circle formed by the intersection of the exobase sphere 
and a plane that is perpendicular to the satellite orbit plane and that con- 
tains the centers of Io and Jupiter. Sodium is ejected uniformly from the 
band area, which is about 63% of the active exobase area of the symmetric 


source . 
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Model results for the sodium cloud with Io at both eastern and western 
elongations are shown in Figure 13 for a symmetric and band source. Compari- 
son of the model results at lower brightness levels for the symmetric source 
indicates a forward cloud that is more elongated at western elongation. This 
east-west orbital asymmetry follows directly from the combined effects of the 
orbital perturbation by solar radiation pressure of sodium atoms (Smyth 1979, 
1983) and the spatially nonuniform lifetime of sodium in the plasma torus. 

The typical east-west orbital asymmetry observed in ground-based images 
(Goldberg et al. 1980; Goldberg, Garneau, and LaVoie 1984) is, however, gene- 
rally more pronounced, and the forward cloud is not this extended in the 
west. For the band source, this preferential elongation of the west forward 
cloud is more enhanced and the west cloud is furthermore foreshortened. These 
additional effects occur because the band source restricts atom trajectories 
to those which are more strongly influenced by the acceleration of solar ra- 
diation pressure. Typically observed east-west orbital asymmetries are prob- 
ably best matched by a source that is somewhat intermediate to the band source 
and symmetric source. Comparison of model results for the symmetric and band 
sources in Figure 13 indicates spatial morphologies that are similar at the 1 
to 2 kiloRayleigh level, but are significantly brighter very near Io for the 
symmetric source. Near Io, this increased brightness occurs because the sym- 
metric source is ~1.6 times larger than the band source, while further from Io 
much of this larger source has since been ionized by the plasma torus and no 
longer contributes to the forward cloud. Model calculations using the band 
source and the plasma torus description with a magnetic dipole offset (Case B) 
are extremely similar to the results in Figure 13, as may be anticipated by 
the small lifetime differences for Case A and Case B shown in Figure 4. 

Model calculations for the band source with a plasma torus description 
including both the magnetic dipole offset and an east-west electric field off- 
set (Case C) are shown in Figure 14. The nominal value for the east-west 
stretched offset adopted in Figures 8 and 10 is assumed (e — 0.025 with 

— n; see the Appendices). For a plasma torus element 5.9 Rj west of the 
dipole center, this nominal value corresponds to an east- west L- shell change 
of -0.30 Rj (i.e., an effective displacement of -0.15 Rj). Comparison of the 
spatial morphology of the model calculations in Figure 14 indicates that at 
eastern elongation the cloud is significantly brighter (-25%) very near Io and 
has a brighter and hence more fully developed forward and trailing cloud with- 
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in a distance of -1/2 to 1 Rj from the satellite. This character of the 
sodium cloud is also evident in the observational image data published by 
Goldberg, Garneau, and LaVoie (1984) and reinforces the idea that the east- 
west electric field is present in the Jovian system. The brighter cloud for 
Io at eastern elongation occurs because of the larger lifetime of sodium that 
exist at eastern phase angles, as was illustrated in Figures 8 and 10. The 
significantly brighter sodium cloud for Io at eastern phase angles as il- 
lustrated in Figure 14 provides a natural explanation for the east-west inten- 
sity asymmetry discovered by Bergstralh et al. (1975, 1977) which is discussed 
below. 

b) East-West Intensity Asymmetry 

The presence of an east-west intensity asymmetry in the brightness 
of the D-line emissions from the sodium cloud very near Io was discovered by 
Bergstralh et al . (1975, 1977) from ground-based observations. These obser- 

vations, acquired for many different east and west geocentric phase angles of 
the satellite over two apparitions, measured the brightness of the cloud in a 
3x8 arc second slit centered on Io with the long dimension oriented perpen- 
dicular to the projection of the satellite orbit plane. As illustrated in 
Figure 15, this viewing slit samples primarily the sodium D-line emissions 
from the local (i.e., bound) atmosphere of the satellite which is within the 
Lagrange sphere. The Lagrange sphere is an effective sphere centered on Io 
inside of which the gravitational field of the satellite dominates the gravi- 
tational field of Jupiter. The slit measurements of Bergstralh et al. show 
that both the D^-line and D 2 -line emission brightnesses for Io east of Jupiter 
are consistently about 25% larger than for Io west of Jupiter. Furthermore, 
the slit- averaged D 2 brightness values for Io near eastern elongation were 
typically -40-45 kR, with a few measured values in the 50-60 kR range. 

Modeling results of Smyth (1979, 1983) suggest that a significant part of this 
east-west intensity asymmetry may be attributed to the acceleration effects of 
solar radiation pressure on sodium atoms if the velocity distribution at the 
exobase is sufficiently peaked at low subescape speeds. Results presented 
below for the sodium cloud model of section 3, however, suggest that the velo- 
city distribution may not be this highly peaked at these very low speeds and 
that the east-west asymmetries in the plasma torus induced by the east-west 
electric field are therefore primarily responsible for the observed east-west 
intensity asymmetry in the sodium cloud. 
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Model calculations presented below for the D-line emission brightness in 
the Bergstralh slit are computed for the Case C description of the plasma 
torus which explicitly includes both the magnetic dipole offset and the east- 
west offset resulting from a convective east-west directed electric field. 

For these model results, the east-west brightness ratio in the slit is mea- 
sured by dividing its calculated brightness for Io at eastern elongation (90 p 
geocentric phase angle) by its calculated brightness for Io at western elonga- 
tion (270° geocentric phase angle). Calculations are performed for the two Io 
System III magnetic longitude angles of 140° and 320° . For an Io System III 
magnetic longitude angle near 140°, the east-west brightness ratio in the slit 
has a maximum value. For an angle near 320°, the ratio has a minimum value. 
This dependence of sodium brightness and hence east-west ratio on the Io 
System III magnetic longitude angle occurs because of the changing space- time 
lifetime history experienced by sodium atoms near Io (see Figure 10) . This 
changing lifetime dependence is geometric in nature and occurs because Io's 
position in the plasma torus (see Figure 8) varies as the orientation of the 
magnetic dipole offset and the east-west offset cycle between alignment and 
near cancellation. This brightness modulation is exhibited in the data of 
Bergstralh et al. (1975, 1977), although it was not recognized by them because 
of an error in properly calculating the Io System III magnetic longitude 
angle. A detailed analysis of the Bergstralh et al. data augmented by addi- 
tional slit data acquired more recently by Goldberg (1987) is beyond the scope 
of this paper and will be published separately. Model computations presented 
here are selected to illustrate the primary character of the east-west bright- 
ness ratio and its relationship to the strength of the electric field and to 
the speed dispersion of atoms ejected from Io's exobase. 

For exobase ejection of sodium, the band source is selected for all model 
caluclated values of the east-west brightness ratio. This source represents a 
somewhat minimal active area exobase source (i.e., -63% of the symmetric 
source) that was shown in the previous subsection to be needed to simulate 
reasonably well the basic morphological features of the forward sodium cloud 
for speeds in the -2-5 km sec“^ range. Both lower speeds (with a likely pre- 
ference to more symmetric ejection) and higher speeds (that do not contribute 
significantly to the forward cloud) may also be present in the velocity dis- 
persion of atoms ejected from the exobase and could contribute significantly 
to the Bergstralh slit brightness. To explore these possibilities, the 
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character of monoenergetic model calculations is examined first. Model cal- 
culations are then performed for different velocity distributions in order to 
assess whether the observed slit brightnesses and east-west brightness ratio 
place constraints on the speed dispersion and source strength of sodium atoms 
ejected from the exobase. 

The model calculated dependence of the east-west brightness ratio on the 
plasma torus shift parameter e (see the Appendices) is shown in Figure 16 for 
monoenergetic (2.6 km sec~^) ejection of sodium from an exobase band source 
and for the two Io System III magnetic longitude angles of 140° and 320° . 

Such a sodium source is, of necessity, required, since the existence of the 
forward cloud source is well established. The parameter e, which is small, is 
defined by Barbosa and Kivelson (1983) as the ratio of the magnitude of the 
east-west electric field to the value of the corotational electric field 
(-0.113 volt/m) at Io's orbital distance from the magnetic dipole. For in- 
creasing values of the east-west electric field (i.e., increasing values of 
c), the east-west brightness ratio increases for both curves in Figure 16. 

For an c value of 0.025, the observed -25% enhancement of the east-west inten- 
sity ratio is achieved as an average between the maximum ratio value on the 
140° curve and the minimum ratio value on the 320° curve. This value of c 
implies an effective average east-west offset of -0.15 Rj for a plasma parcel 
at Io's distance from the magnetic dipole center and an east-west electric 
field of -2.8 mV/m. This value of e is very similar to the -0.03 value sug- 
gested by Barbosa and Kivelson (1983) to produce a -±10% modulation in the 
electron temperature that is required in explaining the dawn- dusk plasma torus 
asymmetry in the EUV S' 1 "* ion emission intensity measured in 1979 by the 
Voyager spacecraft (Sandel and Broadfoot 1982; Shemansky and Sandel 1982). 
Additional analysis of these data (Sandel 1985) has shown that the EUV asym- 
metry is time variable and corresponds to a time-averaged east-west offset of 
a plasma torus parcel at Io's distance from the magnetic dipole center of 
about 0.2 Rj, roughly toward dawn. Further confirmation of an east-west 
electric field of this magnitude is evident in east-west asymmetries in op- 
tical S + ion emissions from the plasma torus observed in the 1981-84 time 
frame (Morgan 1985a, b,c). Analysis of 1981 data by Morgan (1985a, b), when the 
plasma torus had changed in some important ways from the 1979 Voyager picture, 
led him to suggest a preferred value of 0.04 for e, with a value near 0.03 
still providing a reasonable fit to the intensity asymmetry but a less 
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preferred fit to differences between the east and west radial profiles of the 
S emissions. Since the Case C description for the plasma torus used in the 
sodium model is based upon the Voyager encounter picture, a representative 
value of 0.025 is adopted for e in subsequent calculations. 

The east and west D 2 brightness (per unit sodium atom flux) in the 
Bergstralh slit is shown in Figure 17 for monoenergetic model calculations 
with velocities ranging from 0 to 10 km sec' 1 and for Io System III longitude 
angles of 140° and 320° . The corresponding east-west brightness ratios are 
shown in Figure 18. The increase and then decrease of the east and west 
brightness and their ratio with ejection speed may be understood as follows. 
For ejection speeds less than about 1.7 to 1.8 km sec' 1 , the trajectories of 
sodium atoms are basically ballistic in nature (i.e., confined to the Lagrange 
sphere; see Figure 15). In this velocity range, the sodium abundance and 
brightness in the Bergstralh slit increase with increasing speed because the 
atoms fill an increasing portion of the slit viewing area. Because of the 
larger sodium lifetime near Io at eastern elongation, this increasing sodium 
abundance translates directly into a relative increase of the east intensity 
and hence to an increasing east-west intensity ratio. For ejection speeds 
larger than about 1.7 to 1.8 km sec 1 , sodium atoms may escape the satellite 
(i.e., the Lagrange sphere) with increasing ease. This reduces both the 
abundance, of sodium near Io and the residence time of sodium atoms in the 
Bergstralh slit, during which time the differing east-west lifetimes may alter 
the sodium abundance differently. These two factors cause the D 2 brightness 
and the east-west ratio to decrease with increasing velocity. In addition, 
for ejection velocities near 2.0 km sec 1 , the abundance of sodium near Io at 
eastern elongation is also preferentially enhanced by the action of solar ra- 
diation pressure (Smyth 1979, 1983), as is particularly evident in the 140° 
curve of Figure 18. 

For the 140° curve and 320° curve in Figure 18, the average east-west 
brightness ratio is equal to the observed 1.25 value for the ballistic ejec- 
tion speed of -1.05 km sec 1 and for the escape ejection speed of -2.6 km 
sec , the latter speed being consistent with the results of Figure 16. For 
an exobase ejection speed of —1.05 km sec 1 , the sodium atoms and hence D 2 
brightness are confined within a spherical shell about 1600 km above the sur- 
face of Io , while for the larger speed of 2.6 km sec 1 , the sodium escapes Io 
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and forms a forward cloud as discussed in the previous subsection. In 
reality, of course, there is a dispersion of sodium atom speeds at the exo- 
base. The slit brightness and corresponding east-west ratio then depend upon 
the convolution of the exobase flux distribution (as a function of speed) with 
' the monoenergetic model calculated brightness per unit flux results of Figure 
17. The nature of the speed dispersion of sodium near Io and its impact on 
the east-west intensity ratio are addressed in the remainder of this section. 

Various components of the speed dispersion of sodium atoms at the exobase 
of the satellite have been roughly characterized in past analyses of spectral 
and spatial data for the D-line emissions. Spectral line profile data 
(Trafton 1975; Trauger 1977; Trafton and Macy 1977, 1978; Carlson et al. 1978; 
Trafton 1983) have shown that the ejection speeds range from values comparable 
to the exobase escape velocity of the satellite (-2 km sec *") to those $s 
large as 15-20 km sec‘\ with data of Brown and Schneider (1981) and also 
more recent unpublished and high signal- to -noise data acquired at Oak Ridge 
Observatory extending this speed range to beyond 70 km sec’^. The typical 
line profile has two main characteristics: a low-velocity peak and a high- 

velocity asymmetric skirt. The peak, which is relatively narrow (FWHM of 
-±20-40 ma or -±1 to 2 km sec'^), is Doppler shifted -20-50 ma (i.e., an ap- 
parent speed of 1 - 2.5 km sec”'*') relative to the projected orbital speed 
of the satellite on the sky plane and may be associated with nonthermal, 
low- speed (2-3 km sec’*") atoms ejected from the exobase. The skirt, which is 
attached asymmetrically to the peak in the Doppler direction of velocities 
larger than the satellite orbit speed, decreases with increasing Doppler dis- 
placement from the peak and extends at low D 2 brightness levels to at 
least -1400 ma (-70 km sec'*-). The peak and skirt may be associated with non- 
thermal processes involving multicollisional neutral -neutral encounters 
as well as single ion-neutral elastic and charge exchange encounters that 
are driven by plasma torus interactions with Io's local atmosphere. Many of 
these speed components in the spectral data may also be identified in the spa- 
tial morphology data for the sodium cloud. The forward cloud may be iden- 
tified with the -2-5 km sec’*' speed range (Smyth and McElroy 1977, 1978; 

Matson et al . 1978; Smyth 1979; Goldberg et al. 1980; Macy and Trafton 1980; 
Smyth 1983), the fainter trailing-cloud directional features with the 
-10-25 km sec’* - speed range (Pilcher et al . 1984), and the very faint 
forward-cloud jets with the -30-80 km sec *" speed range (Trauger 1984, 1985). 
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To capture the essence of this speed dispersion and to probe the exis- 
tence of a more traditional thermal source, two different velocity distribu- 
tions (or, more correctly, flux distributions) are considered as a function of 
the initial ejection speed v of sodium atoms from the exobase. The first is 
what shall be called a modified- sputtering flux distribution, denoted by 
^( v ; v b , a, v M ) , to describe in an approximate manner the nonthermal source 
that is produced by the interaction of plasma torus ions with Io's atmosphere 
and exosphere. For the modified sputtering flux distribution, the velocity 
v b is Primarily responsible for determining the most probable speed of the 
distribution, the parameter a regulates the breadth of the dispersion re- 
sulting from the composite effects of single ion-neutral collisional ejection 
of sodium and sodium ejection from ion-neutral, neutral -neutral collisional 
cascade processes, and v^ is the maximum speed that a stationary atom can ac- 
quire by collision with a corotating plasma torus ion. The second dis- 
tribution is a Maxwell -Boltzmann flux distribution, denoted by ^(v; T) , to de- 
scribe a thermal source and its thermal corona, where T is the exobase temper- 
ature. The definitions and properties of these two flux distributions are 
presented in Appendix D. The two flux distributions are expressed for con- 
venience as a product of the total flux 4>q for all speeds and a normalized 
flux distribution <f> 


Hv; v b , a, v M ) - <jV , o * (v . v v 

E 

R 2 

T > - (r^) 4>q *(v; T) 
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(12. a) 
(12. b) 


where <f> is constructed so that it is equal to unity when integrated over the 

complete speed range and where is referenced to the satellite surface of 

radius R<, and not the exobase surface of radius R-.. 

E * 


The normalized modified sputtering distribution is illustrated in 
Figure 19 for a — 3 (classical sputtering) and a - 7/3 (Thomas -Fermi -like col- 
lisional ejection), where is selected so that the most probable speed is 3 

km sec’ 1 . The maximum speed of the flux distribution is 7v 2 - v 2 , where Vluf 

M b M 

* 6.62 km sec is determined for collision of corotating oxygen ions and 

stationary sodium atoms. For corotating sulfur ions and stationary sodium 
atoms, v M (* 66.16 km sec’ 1 ) is larger, but the shapes of the two flux distri- 
butions are otherwise very similar to those in Figure 19 and therefore need 
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not, in addition, be considered for our purposes. The significant difference 
between the two curves in Figure 19 is that the distribution for a - 3 is more 
populated near the most probable speed and consequently less populated at 
higher speeds ( > 12 km sec' 1 ) than the distribution for a - 7/3. As noted 
earlier by Sieveka and Johnson (1984) , the approximate two-to-one ratio in the 
sodium source for the forward cloud and directional feature that was deduced 
by Pilcher et al. (1984) is more appropriate for the Thomas -Fermi -like flux 
distribution. In Figure 19, this is illustrated simply by roughly assigning 
the forward cloud to the 2-5 km sec 1 speed range and the directional feature 
to the 12-25 km sec' 1 speed range. For the distribution with a - 3, the atom 
flux populations associated ■with the forward cloud and the directional feature 
are then 47.2% and 6.8%, respectively, which yield a sodium source ratio of 
-7. For the distribution with a - 7/3, these two populations are 34.6% and 
16.5%, respectively, and yield the desired sodium source ratio of -2. 

The normalized Maxwell -Boltzmann flux distribution is illustrated in 
Figure 20 for a moderately high exobase temperature of 1500°K and compared 
with the modified sputtering flux distributions of Figure 19. The Maxwell- 
Boltzmann flux distribution is highly peaked about the most probable speed of 
1.275 km sec' 1 . About 80% of the flux population is confined below the 
speed- of 1.8 km sec' 1 , while only 11.7% and 1.4% have speeds in excess of 
2.0 km sec ' 1 and 2.6 km sec' 1 , respectively. Such a flux distribution 
will produce a sodium corona confined almost entirely within the Lagrange 
sphere of Io and will consequently contribute to the Bergstralh slit bright- 
ness. This Maxwell -Boltzmann flux distribution is, however, clearly not ener- 
getic enough to be an appropriate source for supplying sodium to the forward 
cloud and to the directional features . 

Model results for the D 2 brightness per unit atom flux (i.e., normalized 
by * Q ) based upon the band source with a modified sputtering flux distribution 
are shown for the Io System III longitude angles of 140° and 320° in Figure 21 
for a - 3 and in Figure 22 for a - 7/3. The two corresponding east-west 
brightness ratios are shown in Figure 23. Both the brightness per unit flux 
and the ratio values are presented as a function of the most probable speed v m 
of the flux distribution, which is related to v b (see Appendix D) as indicated 
numerically in Figure 23. In Figure 23 the observed east-west ratio of -1.25 
is matched as an average value between the 140° and 320 curves for a most 
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probable speed of -3.2 km sec 1 with a - 3 and for a most probable speed of 
-2.9 km sec 1 with a « 7/3. Assuming a source rate of 2.3 x 10® atom cm"^ 
sec 1 for the forward cloud (i.e. f in the 2-5 km sec’^ speed range), the flux 

<(>q is then 5.05 x 10® atom cm ^ sec 1 and 6.53 x 10® atom cm"^ sec"^ for a — 3 

and a - 7/3, respectively. With these values for <f> Q and the most probable 
speeds of 3.2 km sec 1 and 2.9 km sec 1 , the corresponding D 2 brightnesses in 
the Bergstralh slit from Figures 21 and 22 vary from -28 kR to 21 kR for a - 3 

and from —32 kR to 24 kR for a — 7/3. For the favored flux distribution with 

a “ 7/3, these D 2 brightness values are similar to, but smaller than, the 
typical observed values of -40-45 kR and the less frequently observed values 
of -50-60 kR. If, however, a symmetric source had been used instead of the 
band source, the model calculated D 2 brightness values would be enhanced by 
about a factor of two, yielding D 2 brightness for a — 7/3 in the -50-60 kR 
range. A source with active exobase surface area somewhere between the band 
source and the symmetric source will therefore match both the observed D 2 
brightness and the east-west intensity ratio of -1.25. 

Model results for the D 2 brightness per unit atom flux (i.e., normalized 
by 4*q) based upon the band source with a Maxwell -Boltzmann flux distribution 
are shown in Figure 24 for the Io System III longitude angles of 140° and 
320°. The corresponding east-west intensity ratio is shown in Figure 25. 

In both figures, the model results are presented as a function of the most 
probable speed v m of the flux distribution, which is defined in Appendix D and 
is related to the exobase temperature T as indicated numerically in Figure 
25. The average of the 140° and 320° curves in Figure 25 matches the observed 
ratio of 1.25 for a low most probable speed of 0.71 km sec" 1 (i.e., T « 
460°K) and for a high and thermally unrealistic most probable speed of -3.65 
km sec 1 (i.e., T ~ 12,300°K). For exobase temperatures of 1000°K to 2000°K, 
which amply cover the range that has been suggested for Io (Kliore et al. 

1975; Kumar 1984, 1985; Summers 1985), the east-west ratio in Figure 25 as de- 
termined by the average between the 140° and 320° curves is significantly 
higher than the observed value of -1.25. This result follows directly from 
the small dispersion of the Maxwell -Boltzmann flux distribution illustrated in 
Figure 20 and the monoenergetic model results of Figures 17 and 18 which are 
significantly peaked near the ejection speeds of 1.6 km sec’ 1 and 1.7 km 
sec" 1 . 
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The fact that the east-west intensity ratio for a Maxwell-Boltzmann flux 

distribution with an exobase temperature of 1000°K to 2000°K is much larger 

than the observed value places a constraint on the flux of such a thermal 

o _ o - 1 

source. For a value of <f>Q no larger than 2 x 10 atom era’ sec , the D 2 
brightness in the Bergstralh slit as determined from Figure 24 varies between 
-27 kR and -18 kR. These brightnesses are comparable to those determined ear- 
lier for the modified sputtering flux distribution, and, when added to the 
contribution for the preferred modified sputtering flux distribution (a - 7/3, 
v m - 2.9 km sec* 1 , <f> Q - 6.53 x 10 8 atoms cm* 2 sec* 1 ), yield an average east- 
west intensity ratio of 1.31 and a maximum eastern brightness of 59 kR, both 
of which are already too large. A value of <f>Q - 1 x 10 cm’ sec is only 
slightly more acceptable, producing a maximum eastern D 2 brightness of -45 kR 
for the two flux distributions and an average east-west ratio of 1.29. Even 
for this value of <f>Q, however, the Maxwell-Boltzmann flux is already dominated 
by at least a factor of -7 at the adopted exobase radius (2600 km) by the non- 
thermal modified- sputtering flux. For a lower exobase temperature, the 

thermal flux is not so tightly constrained. With an exobase temperature of 

8 - 2 - 1 

460°K, the thermal flux may be as large as ^ « 2 x 10 atoms cm sec , 
since, when it is combined with the preferred modified sputtering flux distri- 
bution, a maximum eastern D 2 brightness of -45 kR and east-west brightness 
ratio of 1.25 are produced. This cooler thermal component of the corona is, 
however, very closely confined near the exobase ( < a few hundred kilo- 
meters). Even though the sodium data do not rule out the existence of such a 
cool thermal component, its presence is not required since, as noted earlier, 
the nonthermal flux distribution by itself is perfectly capable of matching 
both the observed brightness and the east-west ratio. These results suggest 
that the ion-neutral collisional processes produced by the interactions of the 
corotating plasma torus ions and Io's atmosphere/exosphere dominate over 
thermal processes in the exosphere and even bring into question the physical 
possibility that Io's atmosphere can support any appreciable classical thermal 
exosphere . 

The density of sodium at the exobase radius (2600 km) may be approximate- 
ly estimated from the <f>Q values for each flux distribution by using the ex- 
pressions derived in Appendix D. For the preferred modified sputtering flux 

distribution (a « 7/3. v m 2.9 km sec , ^ " 6.53 x 10 atoms cm sec ), 

u 3 

the nonthermal sodium density at the exobase is estimated to be —3.0 x 10 
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. 3 

cm . This value is very similar to the sodium density of (2.2 ± 1.1) x 10 3 

cm* 3 actually measured at this distance from Io by means of novel eclipse data 

for Io's sodium atmosphere obtained from ground-based observations (Schneider 

et al. 1987). Nontherraal model results to also match the radial density 

gradient of sodium, measured by Schneider et al., look very promising, but 

such an analysis is beyond the scope of this paper. For the Maxwell-Boltzmann 

flux distribution, the thermal sodium density is constrained to be < 1.7 x 10 3 

cm* 3 for a hot thermal exobase (T-1500°K, <f> 0 < 1 x 10 8 atoms cm* 2 sec* 1 ) 

3 3 u 

and < 6.0 x 10 cm* for a cooler thermal exobase (T - 460°K, <j>n < 2 x 10 8 

- 2 - 1 " - 
atoms cm sec ) . When compared with the measured sodium density , these 

model-constrained thermal densities suggest that cooler thermal components are 

not present at the adopted exobase radius, although they could, of course, 

exist at a somewhat lower altitude if the true collisional exobase were 

located there. As has been suggested by the relative ion compositions in the 

plasma torus (Gehrels and Stone 1983; Bagenal 1985b) and by the relative 

escape rates of atoms to supply the neutral clouds (Smyth and Combi 1984), the 

relative abundance of sodium to oxygen and sulfur is -10' 2 . If this is 

indicative of the mixing ratio at the model-adopted exobase radius, then the 

density of the major gas species at this distance from Io would be ~10 5 
-3 

cm . Such a density is, however, at least an order of magnitude smaller than 
the density required to form a classical collisional exobase, suggesting even 
more strongly that the true collisional exobase is likely somewhat below the 
model-adopted exobase radius of 2600 km. 
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5. Discussion 

The new model for Io's neutral gas cloud developed in Paper I has been 
applied to the Io sodium cloud. Three successively more complex descriptions 
of the plasma torus were explicitly considered for the sodium sink in addition 
to sodium sources at Io's exobase that were both thermal and nonthermal in na- 
ture. These sink and source signatures have been increasingly more evident in 
observational data for the sodium cloud in the past several years. A review 
of earlier sodium cloud models indicated that previous descriptions of the 
plasma torus and the sodium source have not adequately included these signa- 
tures. Model calculations to illustrate the importance of these sinks and 
sources in the sodium cloud were presented. These calculations were able both 
to reproduce and to provide a consistent physical explanation for the basic 
morphological characteristics and east-west intensity asymmetry of the cloud 
observed from ground-based telescopes. These model results and their physical 
implications are summarized below. 

Model results for the distribution of the D 2 “line emission brightness on 
the sky plane were used to analyze the general observed characteristics of the 
sodium cloud and to relate them to both the spatial structure of the plasma 
torus and certain source conditions on Io's exobase. The most basic observed 
characteristic of the sodium cloud — the absence of a bright, elongated 
trailing cloud and the presence of a predominantly bright forward cloud — 
was shown to occur naturally because of the extreme radial gradients in the 
density and temperature of the plasma torus inside Io's orbit. These grar 
dients produce a sharp increase in the electron impact ionization lifetime for 
sodium atoms (see Figure 1), so that the forward cloud (which dynamically oc- 
cupies the interior circumplanetary space) survives much longer than the 
trailing cloud (which dynamically occupies the external circumplanetary 
space). The overall shape and brightness of the basic forward and trailing 
clouds are consistent with an exobase source that contains radial sodium ejec- 
tion from at least a portion of the inner and outer hemispheres. Both a sym- 
metric source (-1.0 x 10 26 atoms sec' 1 ) and a band source (-0.65 x IO 26 atoms 
sec" 1 ) with a mean, nonthermal ejection speed of -2.6 km sec 1 produce a simi- 
lar and acceptable forward cloud for brightness levels of < 1-2 kR. In the 
model calculation, the east-west orbital asymmetry of the forward cloud, which 
is produced by D-line solar radiation pressure effects on sodium atoms, is, 
however, more enhanced for the band source than for the symmetric source. 
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When compared with typical image data (Goldberg et al 1980; Goldberg, Garneau, 
and LaVoie 1984) , a preferred source would appear to be somewhat intermediate 
to the band and symmetric sources. The generally brighter sodium cloud that 
has been observed in image data when Io is preferentially east of Jupiter 
(Goldberg, Garneau, and LaVoie 1984) occurs naturally in the model results for 
the most complex of the plasma torus descriptions (Case C) , where both an 
east-west electric field offset and a magnetic field dipole offset were 
included. 

Model results were also presented that quantitatively reproduce the ob- 
served east-west intensity asymmetry of the sodium cloud (Bergstralh et al. 
1975, 1977) and relate it to an east-west’ (i . e . , dawn-dusk) electric field in 
the Jupiter system and to the flux distribution for sodium atoms at the model - 
adopted exobase of Io. For a plasma torus description (Case C) with a mag- 
netic field dipole offset and an east-west electric field of -2.8 mV/m (i.e., 
an east-west offset of 0.15 Rj at Io's distance from the dipole center), mono- 
energetic (-2.6 km sec‘^) model calculations reproduced the basic structure of 
the forward cloud and also the average 25% east- to -west brightness enhancement 
observed in the measuring slit of Bergstralh et al. Larger and smaller values 
of the electric field produced enhancements that were, respectively, too large 
or too small when compared with the observed value. A similar value for the 
east-west electric field has been previously used to explain both the EUV 
dawn-dusk intensity asymmetry from S++ emission in the plasma torus (see 
Barbosa and Kivelson 1983; Ip and Goertz 1983) measured by the Voyager space- 
craft and the optical east-west intensity asymmetry from S + emission measured 
from ground-based facilities (see Morgan 1985a, b,c). The east-west asymmet- 
ries of the plasma torus and the sodium cloud are anticorrelated, as noted 
earlier by Smyth and Combi (1985 and 1987). To properly simulate the D-line 
brightnesses in the Bergstralh et al. slit, however, a flux distribution must 
be used to include the contributions from the speed of all sodium atoms in the 
exosphere/corona, since the slit is centered on Io and spatially samples the 
sodium brightness primarily in the Lagrange sphere of the satellite (see 
Figure 15). A modifed sputtering flux distribution, to describe the nontherm- 
al speeds evident in sodium cloud data, and a thermal Maxwell -Boltzmann flux 
distribution, to explore the importance of a classical thermal exosphere, were 
therefore adopted. 
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For the modifed sputtering flux distribution, two different physical 
situations were included to describe the nonthermal sodium source produced by 
interaction of the corotating plasma torus ions and the atmosphere/exosphere 
of Io. The first is for a sodium source generated by a complete ion-neutral, 
neutral -neutral cascade process, for which a classical sputtering flux distri- 
bution was chosen- The second is for a sodium source generated by a somewhat 
incomplete cascade process that exhibits at higher speeds ( > 10 km sec the 
character of single elastic collisional ejection of atoms by corotational 
ions, for which a Thomas -Fermi -like flux distribution was constructed. 

This second flux distribution more accurately reflects the relative speed 
populations of sodium atoms as roughly deduced from spatial and spectral data 
for the sodium cloud. Fast ( > 60 km sec’^) sodium atoms, which are produced 
by charge exchange with Na + and are evident in the very faint sodium jets 
ahead of Io (Trauger 1984,1985), were not, however, included along with the 
Thomas -Fermi -like flux distribution, because they are too fast and are too low 
in population to be important for our immediate purposes. From model calcula- 
tions, it was shown that the general spatial character of the sodium cloud and 
.both the average -25% east^to-west intensity enhancement and the absolute 
brightness measured by Bergstralh et al. were consistently reproduced by a 
Thomas -Fermi -like distribution with a most probable speed of -2.9 km sec ^ , an 

O _ O _ 1 

exobase source with flux of 6.5 x 10 atoms cm sec (referenced to Io's 
surface area) , and an active exobase area somewhere between the band source 
and the, symmetric source. This corresponds to a total sodium source between 
1.7 x 10 26 atoms sec' 1 and 2.7 x 10 26 atoms sec' 1 . The classical sputtering 
flux distribution for a symmetric source could just fit both the brightness 
and 25% east-to-west intensity enhancement, but it is much too deficient in 
flux in the 12-25 km sec"*** speed range to supply adequately the fainter sodium 
directional feature (Pilcher et al. 1984). 

Model results for the Maxwell -Boltzmann flux distribution indicated that 

for suggested thermal exobase temperatures in the range - 1000-2000°K, the 

calculated east-west intensity ratio is much higher than the observed value 

o 

of -1.25. This result constrains the thermal flux to be < 1 x 10 atoms 
cm' 2 sec' 1 , even when a nonthermal Thomas -Fermi -like flux distribution with 

O 

a minimum band source and flux of 6.5 x 10 is also assumed to contribute 
simultaneously to the slit D 2 brightness. The thermal flux distribution con- 
tributes primarily to the corona density within the Lagrange sphere and is far 
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too deficient in energetic sodium atoms to contribute any significant sodium 
to the forward cloud. The near order of magnitude dominance of the nonthermal 
flux over the thermal flux at the model-adopted exobase (2600 km radius) cor- 
responds to a thermal sodium density of <1.7 x 10® cm ~® and a nonthermal so- 
dium density of -3.0 x 10 cm . The nonthermal density compares very well 
with the measured value of (2.2 ± 1.1) x 10® cm ® at the model -adopted exobase 
recently reported by Schneider et al. (1987). A cooler exobase temperature of 
460°K is constrained by the sodium data to a thermal flux of < 2 x 10® atoms 
cm sec and a corresponding thermal sodium density of < 6.0 x 10 J cra"^ at 
the model -adopted exobase. Based upon these sodium densities, the classical 
collisional exobase would occur somewhat below the model -adopted exobase. 

Such heights for the collisional exobase are nearer the "minimum'’ SO 2 atmo- 
sphere for Io at the Voyager 1 encounter required by photochemical models of 
Kumar (1984). 

Based upon the sodium results, the escape of oxygen and sulfur from Io 
should also be characterized by a large, nonthermal flux distribution in the 
exosphere. In contrast to sodium, however, an equally large and very fast 
(-60-100 km sec ^) flux distribution will be produced for atomic oxygen and 
sulfur by charge exchange reactions because the elastic and charge exchange 
cross sections are comparable (Smyth and Combi 1984) . The very fast neutrals 
will, however, contribute minimally to the corona density, but the corres- 
ponding newly created hot oxygen and sulfur ions will provide at Io a major 
energy source for the plasma torus (Smyth and Combi 1984) . This energy source 
at Io, in contrast to the spatially distributed energy source for the torus 
produced by the extended neutral clouds, may provide the additional energy re- 
quired to properly balance the EUV energy radiated primarily by sulfur ions as 
noted recently by Shemansky (1987). The lower speed, nonthermal flux produced 
by elastic collisional processes will, however, support an extended nonthermal 
corona. The density of this nonthermal corona will be determined in a self- 
consistent manner by a balance between the outward- directed escape flux of 
c °H*- s l°nally energized neutrals and the upward supply rate from below of 
gases from Io's lower atmosphere and surface. From simulation of such a 
corona, McGrath and Johnson (1987) have reached conclusions similar to those 
of this paper, namely, that these collisional corona processes are the domi- 
nant gas escape mechanism for supplying mass to the plasma torus and that 
single collision ejection of neutrals is, for the "minimum” SO 2 atmosphere of 
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Io (Kumar 1984), the largest coronal ejection process. These nonthermal col- 
lisional processes will penetrate energetically somewhat below the classical 
collisional exobase level and may modify significantly the gas density and 
simple diffusive equilibrium (vertical stratification) assumptions used in 
photochemical models for Io (Kumar 1984, 1985). In this manner, the simple 
photochemical picture could be replaced in the upper atmosphere by a 
nonthermal heating layer and an effective global atmospheric blowoff process. 

Future studies of the sodium clouds of Io will be able to extend and re- 
fine significantly the physical insights for the satellite and the plasma 
torus gained from the investigations of this paper. Detailed model -data com- 
parisons of the brightness and the east-west brightness ratios near Io can be 
used to (1) calculate the sodium density profile inside and beyond the 
Lagrange sphere of Io, (2) estimate the amount of cooling that electrons may 
have undergone in the upper atmosphere/corona, (3) refine the shape and non- 
isotropic character of the nonthermal flux distribution of sodium atoms near 

the exobase and hence the character of the ion-neutral processes in the upper 
atmosphere/corona, and (4) investigate variations in the strength of the east- 
west electric field and/or the plasma torus properties over more than a. decade 
time span. Further model -data studies of time variations in the brightness 
morphologies of the forward sodium cloud and the directional features can be 
used to investigate inherent System III longitudinal dependence in the plasma 
torus and probe changes in the sodium source that occur on time scales which 
depend upbn the satellite phase angle, the System III longitude angle of Io, 
or longer secular changes in the Jovian system. Of particular pertinence to 
these morphological studies are the massive sodium cloud data set acquired by 
JPL (Goldberg, Garneau and LaVoie 1984) and the novel Doppler-resolved images 
for the sodium cloud acquired by Trauger (1984, 1985). 

Beyond the sodium cloud, valuable future studies for the potassium, oxy- 
gen and sulfur clouds can also be performed. In particular, by using the new 
gas cloud model of Paper I and cloud emission brightness data acquired for 
oxygen and sulfur with the aid of rocket (Durrance, Feldman and Weaver 1983; 
Skinner and Durrance 1986), the IUE satellite (Ballester et al. 1987), and 
ground-based telescopes (Brown 1981; Mekler 1986), a number of major investi- 
gations will be made possible. These investigations will be capable of cal- 
culating (1) the density in the upper atmosphere/corona of Io, (2) the loss 
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rate of gases from the atmosphere, and (3) the three-dimensional spatial 
character of the ion loading rate, plasma mass loading rate and ion energy 
input rate for the plasma torus supplied by different gas species lost by Io 
both very near the satellite and in the spatially extended neutral clouds . 
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6 . Summary 

A new general neutral cloud model (Paper I) has been applied to explain 
various observed phenomena regarding both the local and the extended regions 
of the neutral sodium cloud produced by Io. These include the predominant 
forward cloud, the east-west orbital asymmetry of the forward cloud, and the 
density of sodium in the near Io atmosphere/corona and its east-west intensity 
asymmetry. The model accounts for the spatially nonuniform and time -dependent 
sodium sink produced by a 7° inclined corotating plasma torus which in its 
most general form (as required by the sodium data) includes the physical and 
geometrical effects of the Jovian offset-dipole magnetic field and an east- 
west electric field of ~2.8mV/m whose presence and strength are implied by di- 
rect torus measurements. 

Sodium sources considered for Io include nonthermal atmospheric sputter- 
ing and incomplete collisional cascade flux distributions, both produced by 
collisions between torus ions and atmospheric atoms, as well as a thermal 
source. For the sodium cloud, a dominant nonthermal source for the incomplete 
collisional cascade flux distribution (peaked at -3 km sec"^ with a strength 

nr -i nr I 

between 1.7 x 10* D atoms sec"' 1 ' and 2.7 x 10* atoms sec ) is required from an 
active exobase area which is basically spherically symmetric with the most 
forward and backward elements absent. 

This first application of the model to the sodium cloud has uncovered 
a number of important new results regarding the nature of the local Io 
atmosphere/corona and its interactions with the plasma torus, the escape of 
material from Io, and the nature of the plasma torus itself. Exciting future 
applications of the new model to the neutral clouds of potassium, oxygen, and 
sulfur, as well as detailed time -dependent studies of the wealth of sodium 
data, were also discussed. The most important of these future studies are 
(1) the space and velocity distribution of gases in the local atmosphere as it 
relates to the escaping gas and to torus electron cooling, (2) variations in 
the east-west electric field, (3) system III longitudinal asymmetries in the 
plasma torus, and (4) the ion, mass, and energy input rates of the local atmo- 
sphere and the extended neutral clouds to the plasma torus. 
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APPENDIX A 

PLASMA TORUS COORDINATE FRAME TRANSFORMATIONS 

The transformation from a location R - (x,y,z) in the magnetic (or, more 
generally, the plasma torus) frame to a location R — (x,y,z) in the satellite 
fixed coordinate frame at an absolute time t is represented by 


R - T(t) • R + D(t) 


A. 1 


in expression (1) of the text. In matrix form, the transformation A.l is 




r x 1 

y 

- [ T ( t ) ] 

y 

. Z . 


. Z . 


+ [D(t)] 

where the tensor matrix is given by 


[T(t>] 


cos 

■sin(^-S) cos^ 
-sin(^-0) sin^ 


and the column matrix is given by 


sin(^-0) 0 

cos (i>-6) cos 4> -sin^ 

cos sin<f> cos <f> J 


A. 2 


A. 3 


[D(t)] - 



O 

O 


r d cosS^sinC/^-p^) + e (x 2 +y 2 ) sin 

- 

0 cos<f> -sin<0 


-d cos5^cos(/x N -p d ) + e (x 2 +y 2 ) li cos (^-V» s ) 


^ 0 sin<£ co s<f> J 


^-d sin5 J 


A. 4 


The tensor transformation A. 3 includes the tilt angle <f> of the plasma torus 
plane and both the angular motion ^ - iKt) of the plasma torus and the angular 
motion 6 - 0(t) of Io about Jupiter. The vector transformation A. 4 contains 
the effects of the magnetic dipole offset distance d and the stretched-offset 
distance c(x 2 +y 2 ) of the plasma torus proposed by Barbosa and Kivelson (1983) 
because of the presence of an east-west electric field in the Jupiter system. 
The definitions of the plasma and satellite coordinate frames as well as a 
more detailed description of the physical transformation A. 2 and its related 
parameters and variables are discussed below. 


The plasma torus is centered at the origin of the (x,y,z) frame with the 
z-axis normal to the centrifugal plasma plane and inclined by an angle <f> with 
respect to the direction of the Jupiter spin axis. The Jupiter spin axis is 
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normal to the equatorial plane of the planet in which Io's orbital motion is 
assumed to be confined to a circular orbit about the planet. The projection 
of the positive y-axis onto the equatorial plane of the planet occurs at the 
System III (1965) magnetic longitude of the projected north magnetic polar 
axis. The magnetic field is described by an offset tilted dipole. The 
offset-dipole distance d has a direction vector in System III defined by a 
latitude angle S ^ and a longitude angle /i^. The spin axis of the plasma torus 
is parallel to the Jupiter spin axis but is displaced from the center of the 
torus by the dipole offset vector. The spin period of the plasma torus is the 
System III rotation period denoted by T IXI . The orientation of the rotating 
plasma torus at an absolute time t is determined by specifying the geocentric 
phase angle 


*(t) 


^1965 (t l* 


2tt 


^III 


(t - 


V 


A. 5 


of the System III longitude (/^ + *r/2) , where g 5 ( ) *- s P^ ase angle as 

determined in the System III (1965) Jupiter longitude convention for a con- 
veniently chosen absolute reference time t^. The System III central median is 
therefore given by (/i N - n/2 + ij>) . 

The spin axis of the plasma torus is furthermore displaced in the 

equatorial plane from the center of Jupiter by an east-west electric field 

(Barbosa and Kivelson 1983; Ip and Goertz 1983) assumed to be in the equa- 

torial planet plane. The stretched displacement is given by e(x + y ) , 

where c is the ratio of the magnitudes of the east-west electric field and 

the corotation electric field as defined by Barbosa and Kivelson (1983) and 
2 2 h 

(x + y ) is the radial displacement from Jupiter in the equatorial plane 
as measured in the satellite fixed coordinate frame. For generality, the 
orientation of the electric field is assumed not to be exactly east- to -west as 
defined by an earth observer. Its orientation is determined by the geocentric 
phase angle if> s of a symmetry line (perhaps solar wind in origin) that is drawn 
in the equator plane, normal to the electric field and directed toward the 
sunlit portion of the planetary equator. The orbital motion of Io about the 
center of Jupiter is in the same direction as the spin of the plasma torus but 
has a longer rotation period Tj q . The satellite fixed coordinate frame 
(x,y,z) is centered on Jupiter with the (x,y) -plane in the equator plane and 
Io on the positive x-axis. The angular orientation of Io about Jupiter at an 
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absolute time t is determined by the geocentric Io phase angle 8 referenced to 
the time t^ as in A. 5: 

*(t) - 0(t x ) + Y~ (t - tp . A. 6 

Io 

For model calculations in this paper, the following values for the 
transformation parameters are chosen: <f> - 7* , d - 0.131, (5 d , - 

(-8°, 149°), - 202° and - 180°. The values for d, $ d , p d and are 

for the 0^ offset tilted dipole model of the Jovian magnetic field (see Acuna, 
Behannon and Connerney 1983) . Different values for the parameter e are chosen 
as discussed in section 4. 
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APPENDIX B 

OBSERVER'S COORDINATE FRAME TRANSFORMATION 


The transformation from a location R - (x ,y , z ) in the observer's 

^ o o o o 

coordinate frame to a location R - (x,y,z) in the satellite fixed coordinate 
frame at an absolute time t is represented by 


R o * V*) * R 


B. 1 


and used in expression (46) of Paper I. In matrix form, the transformation 
B. 1 is 


> 

^ X 
o o 

J 

- [T o (t)] 

r X 

y 

^ z > 


^ z 


o 


where the tensor matrix is given by 


[T o (t>] - 


cos <f> cos 9 
sin 9 

^ sin <f> cos 9 


-cos <f> sin 9 -sin <f> ^ 
cos 6 0 

-sin^ sin 9 cos <f> 


B.3 


The observer's frame is centered on Jupiter with x Q -axis directed along 
the Earth- Jupiter line and the Earth located on the negative x Q axis. The 
angle <f> is the apparent tilt angle of the orbit plane of Io as seen by an 
Earth observer and 9 — 9( t) is the geocentric phase angle for Io defined by 
equation A. 6 . 
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APPENDIX C 

ELECTRON IMPACT IONIZATION DESCRIPTION FOR SODIUM 


The electron impact ionization description of the lifetime of sodium in 
the plasma torus coordinate frame (x,y,z) is developed below for Case A (a 
rotating, Jupiter-centered and tilted dipole), Case B (a rotating, dipole- 
offset-centered and tilted dipole) , and Case C (a rotating, dipole-offset and 
tilted dipole in the presence of an electric field in the planetary equator) . 
These cases are adopted in section 3 and used in section 4. The descriptions 
developed below are appropriate to the encounter time of the Voyager 1 
spacecraft with Jupiter. 

Voyager 1 in situ plasma data have been reduced (Sittler 1986; Sittler 
and Strobel 1986; Bagenal 1985a, b; Bagenal et al. 1985; Bagenal and Sullivan 
1981) for each measured position along the spacecraft trajectory to L- shell 
coordinates , adopting an offset tilted dipole description for the magnetic 
field. Assuming that the electron and ion temperatures are constant along the 
magnetic field lines and that the electron and ion density distributions along 
the field lines are determined by solving the diffusive equilibrium equations 
for a multispecies plasma (see Bagenal 1985b) , the electron temperature and 
density can be determined in the centrifugal plasma torus equator plane at 
specific radial and longitudinal locations for each spacecraft measurement. 
Because of the lack of sufficient longitudinal plasma information, the simpli- 
fying assumption has generally been made. that these densities and temperatures 
in the centrifugal plane are independent of magnetic longitude. With this 
assumption, radial profiles for the electron temperature and density in the 
centrifugal plasma equator plane follow and are shown, respectively, in 
Figures 5 and 6. This simplifying assumption, of course, ignores the skewed 
(i.e., nonnormal) orientation of the field lines to the centrifugal equator 
plane which varies with system III longitude and causes the real plasma torus 
description above and below the centrifugal plane to be asymmetric except at 
the node longitudes (-/i^±90 o ). This assumption also ignores any intrinsic 
system III magnetic longitude dependence which has been seen in more recent 
plasma torus data obtained from ground based observations (Pilcher and Morgan 
1985; Pilcher, Fertel, and Morgan 1985; Pilcher and Morgan 1980; Trafton 
1980a; Trauger, Munch, and Roesler 1980). Adopting this simplification, a 
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two-dimensional description for the electron density in the plasma torus is 
chosen at the node longitude's projection (an average condition) where the 
density above and below the centrifugal plasma plane is symmetric. This 
longitudinally symmetric description for the electron density is shown in 
Figure 7, and the corresponding lifetime description for sodium based on 
electron impact ionization is shown in Figure 1 . In Case A the radial 
coordinate in Figure 1, which is actually centered on the offset-dipole 
location, is assumed simply to be centered on Jupiter. In Case B the radial 
coordinate in Figure 1 is centered on the dipole -offset location. 

In Case C the additional presence of an electric field in the equator 

plane destroys the longitudinal symmetry of the electron temperature and 

density as depicted in Figures 5-7. As a plasma volume element moves from 

— .o o u 

west to east of Jupiter, its radial displacement L - (x + y ) in the cen- 
trifugal plane from the dipole centered coordinate frame (x,y,z) increases, 
resulting not only in a new center for the essentially circular path of the 
plasma volume element (Barbosa and Kivelson 1983) but also in an expansion in 
the density and a cooling in the temperature of the plasma volume element. 
Ignoring the approximate 7° tilt of the centrifugal equator, which because of 
the rotation of the plasma torus about Jupiter causes the component of the 
equatorial electric field parallel to the radial coordinate L to be modulated 
very slightly, the radial displacement of the plasma volume element is given 
(to first order in e) by 

7 1 + €[1 + sin($ - (i/> - *))] . C.l 

Li S 

W 

Here $ is the phase angle of the plasma volume element about Jupiter ($ - n/2 

at eastern elongation, $ - 3 tt/ 2 at western elongation, and - it for an east- 

to-west directed electric field), is the radial displacement of the plasma 

volume element at $ — 3?r/2 + (^ g - n ) , and c is the ratio of the uniform east- 

west convection and corotation electric fields as defined by Barbosa and 

o o h 

Kivelson (1983) . The radial coordinate L - (x + y ) in the plasma torus 
centered coordinate frame (x,y,z) is related (to first order in e) to the 
dipole centered radial coordinate L by 

L - L[1 - e sin($ - (^ - *))] 

s 


C . 2 



so that (to first order in e) L is, as required, independent of $ - (^ g - tt) 
for a particular plasma volume element 


L - (1+e )L . 

w 


C . 3 


To describe the plasma torus in the L centered coordinate, the descrip- 
tions of the electron density in Figure 7 and the electron temperature in 

Figure 5, both deduced from Voyager 1 measurements made west of Jupiter, are 

w — 

adopted as the description of the electron density n e (L w ,z) and electron 

temperature T^(L w ) at that phase angle near western elongation $ - ^ + (if> s -n) 

where the plasma volume element is in the Jupiter meridian plane containing 

the east-west electric field vector. The electron density n and tempera - 

- _ e 

ture T in the plasma- torus centered coordinate system (L,z,$) are then given 
e 

by 


n (L,E,$) - 
e 




(L ,z) 
w 


C . 4 


T (L,5,$) 
e 




(L ) , 

w 


C . 5 


where the expansion factor is 



C . 6 


so that to first order in e 


f - 3c [1 + sin ($ - (^ s *w))] V C.7 

The expansion factor C.6 for the electron temperature follows from the approx- 
imate integral of the motion for the diamagnetic moment of a charged particle 
in a slowly (space and time) varying magnetic field. This adiabatic invariant 
also requires that the magnetic flux through a flux tube containing the plasma 
volume element be invariant for the west- to -east radial displacement. For a 
dipole field, the cross-sectional area in the symmetry plane of the flux tube 
is proportional to L^. Assuming a constant height for the plasma volume 
element as it is displaced radially, the expansion factor C.6 is therefore 
also applicable for the density C.4. 
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The lifetime ? for sodium atoms in the plasma torus centered coordinate 
frame is given by the reciprocal of the collision frequency, 

»>(L,z,$) - n (£,z,$) K(T (£,$)) , c . 8 

“ G 

where K(T^) is the electron impact ionization rate coefficient. Combining the 
expressions C.4-C.8, expanding the rate coefficient K(T e > about T -T in a 
Taylor series in c, and retaining terms to first order only in e^the lifetime 
for sodium atoms r — is given explicitly by 

f(L,z,$) - r W (L z) [1 + 3c (1 + sin ($ - (if> - *))) (1 + 
w s 

T^(L w ) K'(T^(L w )) / K (T"(L w )))] C. 9 

where L — (1+0 L and z * z. Here, 
w 

r W (L w ,I) - [n^(L w ,I) K (T*^))]’ 1 C.10 

3 7T 

is the lifetime at $ ■= ^ + *■) as given by Figure 1, T^(L w ) is given by 

the adopted curve in Figure 5, and K(T) and K' (T) are obtained from Figure 3, 
where 

K>(T)-4££1. 


C.ll 



47 


APPENDIX D 

THERMAL AND NONTHERMAL FLUX DISTRIBUTION 


The atom flux ejected from the exobase to describe the primary population 
of the cloud is given by 

A A (primary) 

4>(v,n,t h ) - w • O f(w,x b ,t b ) D. 1 

(primary) 

in equation (11. b) of Paper I, where f(td,x b ,t b ) is the one-particle distribu- 
tion function. The rate at which atoms are ejected at time t b from a differ- 
ential surface area element dfl (- R 2 sin $ d$ d9) with an outward- directed 

A E 

normal unit vector Q - (0,$) located on the exobase of radius R £ at - (R E> 8,$) 
and with initial velocity vectors in the range w to w + dtj is then 

^(w,Q,t b ) dO dw. D. 2 

For purposes of this paper, the atom flux is assumed to be uniform on the 
active portion of the exobase, zero on the inactive portion of the exobase, 
and independent of the ejection time t b . The one-particle distribution func- 
tion on .the active exosphere is also assumed to depend only upon the speed 
— ► 

v - |w| , so that dropping the '(primary)' superscript, 

- w • 0 f(|w|), D. 3 

where f(|w[) is normalized so that the number density n of this gas at the 
exobase is given by 

n- J f(|w|) d 3 w, D.4 

all space 

3-+ -+ 2 

where d w - \w\ dfl'dv is the volume element in velocity space and dfl' is the 
solid angle element in velocity space. An upward- directed flux distribution, 

^(v) , is defined as a function of the speed v by 
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*(v) 



J dO' « 7rv^f(v) 

upward 

hemisphere 


D . 5 


so that the uniform flux of atoms ejected from the active area of the exobase 
in the speed range v to v+dv is ^(v) dv. In this paper, the upward- directed 
flux distribution is adopted and used in the model calculations, where it is 
simply assumed that the ejection velocity vectors of all atoms are initially 
normal to the exosphere (i.e., radial ejection). 


For the Maxwell -Boltzmann flux distribution adopted in this paper, the 
one -particle distribution function depends parametrically on the exobase 
temperature T and is given by 


where 


f(v;T) 




-<v/v T ) 2 


D . 6 



D . 7 


is the most probable speed of the velocity distribution and where m is the 
mass of the gas atom. The Maxwell -Boltzmann flux distribution then follows 
from equation D.5 and is given by 


*(v;T) 


^ <-> 3 
h V T 


-(v/v T ) 2 


D. 8 


The most probable speed for the flux distribution D.8 is denoted by v m and is 
given by 


v 

m 



D. 9 


The total upward flux at the exobase may then be related to the flux <f>Q , 
referenced to the surface area of Io, by the expression 

® nv nv R 2 

f *(v;T) dv z--=- 5T D - 10 

0 2 U 76 x R E 

where Rg is the radius of the satellite. The normalized flux distribution 
^(v;T) is then defined by 
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*<v;T) 


nv, 


*(v;T) / 


T 

2jn 


-(v/v T )‘ 


D.ll 


For the modified sputtering flux distribution adopted in this paper, the 
equivalent one -particle distribution function is given as a function of three 
parameters (a,v b ,v M ) by 


f(v;a,v b ,v M ) 


2a- 3 


*C(a,v M /v b ) (V 2 +V 2) a 



D.12 


where the dimensionless function C(a,v M /v b ) is given by the integral 


C(a,v M /v b ) 


7(v M /v b )2 . i 

4 ; — ^ — 
0 (1 + **)“ 


1 - — ^ Jl + x 2 
V M 


dx 


D. 13 


which is not easily expressed in an analytical form for arbitrary a. The 
dimensionless parameter a regulates the breadth of the speed dispersion, and 
describes the composite effects at the exobase of single-collisional en- 
counters and multicollisional encounters (i.e., cascade processes or 
sputtering) of neutral sodium atoms in the atmosphere/corona with corotating 
ions in the plasma torus. Classical sputtering, which describes a collisional 
cascade process that has proceeded to completion, is obtained for a-3 (see 
Sieveka and Johnson 1984). As is indicated by the directional feature data 
for the sodium cloud (Pilcher et al. 1984; Sieveka and Johnson 1984), incomp- 
lete collisional cascade processes for the higher speed (>10 km sec' 1 ) sodium 
atoms are evident and are more appropriately represented by its leading term 
in the cascade chain, which is single elastic collisional ejection. For 
single elastic collisional ejection of sodium atoms by plasma ions (Johnson 
1985), the atom flux distribution is directly proportional to the Thomas-Fermi 
differential cross section (see Johnson 1982; Sieveka and Johnson 1984), which 
for the energy range of interest in the plasma torus and for the mathematical 
form adopted for the modified sputtering flux distribution yields an a value 
of 7/3. The velocity parameter v b in the distribution function D.12 is pri- 
marily responsible for determining the most probable speed of the velocity 
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distribution. In the theory of sputtering from surfaces, the speed v b is 
determined by the surface binding energy (mv b 2 /2) for the ejected atom and 
also ensures that the distribution is well behaved at low values of the speed 
v. The velocity parameter Vjj is the maximum speed that a stationary sodium 
atom of mass m can acquire in an elastic collision with a plasma torus ion of 
mass M and relative speed v^: 


M 


2M 

(M + m) V R* 


D. 14 


The velocity Vj| establishes a maximum speed for the velocity distribution 
which is given by - v b 2 . The form of the term [1 - ,/v 2 + v b 2 /v M ] in 

equation D.12 was introduced by Thompson (1968) to limit the distribution to 
this maximum speed. For model calculations in this paper, v^ is chosen to 
have the nominal corotational speed of 56.8 km sec' 1 relative to Io, but this 
speed is likely to be much larger (>100 km sec' 1 ) as determined from 
simulations of Io's interaction with the plasma torus (Wolf-Gladrow, Neubauer, 
and Lussem 1987) . 

The modified sputtering flux distribution then follows from equation D.5 
and is given by 


^(v;a,v b ,v M ) 


n v_ 3 v b “ 

C(<*,v M /v b ) v b v 2 +v b 2 


v 2+ V 
V 


D. 15 


The most probable speed, v ffl , for the flux distribution D.15 is determined by 
solving the nonlinear expression 


m 


where 


/F 7 ? 


D. 16 


e - i 


i - [v M /(7v m 2 + v b 2 ) . lr 1 . 


D. 17 


The total upward flux at the exobase is then related to the flux ^q, refer- 

enced to the surface area of Io, by the expression 


7V-V 

S ^< v ;«.v b ,v M ) 
o 



D ( a ’ v m/ v h> _ 
c («.v M /v b ) R e 


2 


4 


0 


D . 18 
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where the dimensionless function D(a,v^/v^) is given by 

AW 2 - 1 

D(o,v M /v b ) - ^ [1 - ^ yi + x 2 ] dx D. 19 

J (1 + x 2 ) a V M 
0 

and can easily be integrated analytically. For the case discussed in this 
paper (a - 7/3, v m - 2.9 km sec*^, v b - 2.2159 km sec*^, - 46.62 km sec'^), 
the dimensionless coefficient in equation D.18 is evaluated to be 


D(a, yv 

C(Q -W 


0.4823. 


D. 20 


The normalized, modified- sputtering flux distribution, ^(v;a, v b > v ^) , is 
finally defined by 


- . , *(v;q,v b ,v M ) 

0(v,or,v b ,v M ) - [nv b (D(a,v M /v b )/C(a,v M /v b ))] 


(— ) < 


v b D(a > W v b 


v, z a 
- b — ) 


V 2 + v b 2 



[1 


]. D. 21 
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Table 2 

Reactions for Sodium Atoms in the Io Plasma Torus 


Reaction 


Reaction Time (hr) 


Rl. Na + hu -» Na + + e 400 

R2 . e + Na -+ Na + + 2e >1 

R3 . S + + Na -► S + Na + 

R4 . 0 + + Na -+ 0 + Na + 

R5 . Na + + Na -► Na + Na + 50* 

R6 . S ++ + Na -*■ S + + Na + > 40* 

R7 . 0++ + Na - 0 + + Na + > 300* 

R8. X + + Na -*■ X + + Na (elastic collision) > 30* 


- 3 + 

Assumed an ion density of 50 cm for Na 


.-3 


500 cm for S , 100 cm' for 
0 4 * 4 ", and 1000 cm"^ for X + with an ion-neutral velocity of 60 km sec"'*'; 
cross-sectional information for electron impact ionization is from 
McFarland (1965), McFarland and Kinney (1965) and Zapesochnyi and 
Aleksakhin (1968); cross sections for S 4 ", 0 + , S 4-4 " and O 4 " 4 ” are from 
Johnson and Strobel (1982), for Na 4 ” from Rutherford et al. (1972) and 
for X 4 * - 0 + or S + from new estimates by Johnson (1985). 
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Figure 1 


Figure 2 


Figure 3. 


Figure 4. 


Figure 5. 


FIGURE CAPTIONS 


Electron Impact Ionization Lifetime of Sodium. The lifetime in 
hours is shown as calculated in the plasma torus for plasma 
conditions appropriate to the Voyager 1 encounter epoch. The 
radial and vertical distances are given in units of the planetary 
radius (Rj). 

Interaction of the Sodium Cloud and the Plasma Torus . The 
spatially projected overlap in the satellite plane of the Io 
plasma torus with the solar radiation perturbed cloud shape 
(solid boundary with heavy shading) and with the unperturbed cloud 
shape (short dashed boundary) is indicated for the diametrically 
opposite satellite phase angles as published by Smyth (1983). 

The cloud shape illustrates the envelope of sodium after 20 hours 
of flight time. 

Electron Impact Ionization Rate Coefficient for Sodium. The rate 
coefficient is shown as a function of the electron temperature and 
has been calculated assuming a Maxwell -Boltzmann distribution and 
cross-sectional information from McFarland (1965), McFarland and 
Kinney (1965) and Zapesochnyi and Aleksakhin (1968) . 

Electron Impact Ionization Lifetime of Sodium at Io. The lifetime 
as a function of the System III magnetic longitude angle of the 
satellite is shown at Io's position for a 7° tilted plasma torus 
that is centered on Jupiter (light line, Case A) and centered on 
the offset magnetic dipole location (heavy line, Case B) . 

Temperature in the Io Plasma Torus. The electron temperature as a 
function of radial distance (L-shell) in units of planetary radii 
(Rj) is a representative average from inbound and outbound PLS 
electron data of Voyager 1 (Sittler and Strobel 1986; Sittler 
1986) . The ion temperature is also shown as obtained from the 
corrected analysis of the PLS ion data of Voyager 1 (Bagenal et 
al. 1985; Bagenal 1985b). 
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Figure 6 . 


Figure 7 . 


Figure 8. 


Figure 9 . 


Electron Density in the Plasma Torus. The electron density as a 
function of radial distance (L-shell) in units of planetary radii 
(Rj) is shown in the centrifugal equator as deduced from the PLS 
electron data of Voyager 1 (Sittler and Strobel 1986; Sittler 
1986) and from Voyager UVS remote sensing data (Shemansky 1985). 

Electron Density in the Plasma Torus. The density in units of 

_ o 

cm is shown in a cross-sectional cut of the plasma torus and is 
appropriate to the Voyager 1 encounter epoch. See the text for 
more details. 

Electron Impact Lifetime of Sodium in the Plasma Torus. The life- 
time in units of hours is shown at eastern and western elongation 
in the coordinate frame of the plasma torus (Case C) that is 
tilted by 7°, has a magnetic dipole offset and is in the presence 
of an east-west electric field of -2.8 mV/m. The lifetime at Io's 
location in the plasma torus for these two elongation pictures 
varies periodically as a function of the System III magnetic 
longitude of the satellite and is shown by the heavy closed line. 

Location of Io in the Plasma Torus. The shaded area shows the 
domain of Io's location in the plasma torus as a function of the 
System III magnetic longitude angle and the "effective L-shell" 
(denoted by L' ) for a 10-degree tilted and rotating magnetic 
dipole that is offset from the center of Jupiter and that is in 
the presence of an east-west electric field of -2.8 mV/m. The 
dotted line in the center of the area is the location of Io if the 
east-west electric field were zero. 
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Figure 10. 


Figure 11. 


Figure 12. 


Figure 13. 


Electron Impact Ionization Lifetime of Sodium at Io's Location. 

The lifetime in units of hours at Io's location for the plasma 
torus description of Figure 8 is shown as a function of the System 
III magnetic longitude angle of Io and the geocentric phase angle 
1° f° r Jupiter at opposition. The highest lifetime contour of 
15 hours is centered on an Io System III longitude angle of -160° 
and an Io geocentric phase angle of 90° . 

d 2 tightness of the Sodium Cloud on the Sky Plane. Monoenergetic 

(2.6 km sec model calculations for the D 2 emission brightness 

are shown with Io at eastern elongation for the four labeled 

quadrant sources and for a symmetric source . A sodium source flux 
8 2 1 

of 2.3 x 10 atoms cm sec (referenced to the satellite surface 
area) and the simplest plasma torus description (Case A) are 
assumed. The D 2 contour brightness levels, in units of kilo- 
Rayleighs , are, from outside to inside, 0.2, 0.5, 1.0, 2.0, 5.0, 
10.0, 20.0 and 50.0. See the text for further discussion. 

^2 Brightness of the Sodium Cloud on the Sky Plane. Monoenergetic 
(2.6 km sec model calculations for the D 2 brightness are shown 
with Io at eastern elongation for the four labeled hemispherical 
sources. The sodium flux, plasma torus description and contour 
levels are the same as in Figure 11. 

East- West Orbital Asymmetry in the Sodium Cloud. Monoenergetic 
(2.6 km sec‘l) model calculations for the D 2 brightness of the 
cloud on the sky plane with Io at eastern and western elongation 
are compared for the symmetric source and the band source. The 
sodium flux, plasma torus description and contour levels are the 
same as in Figure 11 with the exception that an outer 0.1 
kiloRayleigh contour has been added. 
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Figure 14. East-West Orbital Asymmetry in the Sodium Cloud. Monoenergetic 
(2.6 km sec''*') model calculations for the brightness of the 
cloud on the sky plane with Io at eastern and western elongation 
are compared for the band model. The sodium flux and contour 
levels are the same as in Figure 11. The plasma torus description 
is for Case C. See the text for further discussion. 

Figure 15. Bergstralh Slit for the Sodium Cloud. The size of the slit 
(shaded rectangle) used by Bergstralh et al. (1975, 1977) in 
measuring the sodium cloud brightness near Io is shown for Jupiter 
at opposition in relation to the satellite disk (dark circle) and 
the Lagrange sphere (large circle) of the satellite. 

Figure 16. East-West Intensity Ratio. The east-west intensity ratio in 
the Bergstralh et al. (1975, 1977) slit is calculated as a 
function of the plasma torus shift parameter e for monoenergetic 
(2.6 km sec"^) ejection from a band source and for the plasma 
torus description Case C. The ratio is shown for a System III 
magnetic longitude angle for Io of 140° (maximum ratio) and 320° 
(minimum ratio) . The dashed line indicates the observed average 
ratio. 

Figure 17. D 2 Brightness in the Bergstralh Slit. The monoenergetic model 
calculated D 2 brightness in the observing slit expressed in 
normalized flux units for Io at eastern (east) and western (west) 
elongation is shown as a function of the sodium ejection speed 
from a band source. The plasma torus description (Case C) is 
described in the text and the System III magnetic longitude angle 
of Io is 140° in (a) and 320° in (b) . 

Figure 18. East-West Intensity Ratio. The east-west intensity ratio, as 
determined from the model results of Figures 17a and 17b, is 
shown. The dashed line indicates the observed average ratio. 
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Figure 19. Modified- Sputtering Flux Distribution. The normalized modified- 

sputtering flux distribution, v, ,v K# ) , is shown as a function 

b M 

of the sodium ejection speed v for a - 3 (classical sputtering) 
and a - 7/3 (Thomas -Fermi -like collisional ejection), where the 
most probable speed for both distributions is 3.0 km sec’^. See 
the text for further discussion. 

Figure 20. Maxwell -Boltzmann Flux Distribution. The normalized Maxwell- 

Boltzmann flux distribution, ^(v;T), is shown as a function of the 
sodium ejection speed v for an exobase temperature of 1500°K. The 
two flux distributions of Figure 19 are also shown for comparison. 

Figure 21. D 2 Brightness in the Bergs tralh Slit. The model calculated 

brightness in the observing slit expressed in normalized flux 
units for Io at eastern (east) and western (west) elongation is 
shown as a function of the most probable speed of the modified- 
sputtering flux distribution for a — 3. The band source and the 
plasma torus description. (Case C) assumed in the calculation are 
described in the text. The System III magnetic longitude angle of 
Io is 140° in (a) and 320° in (b) . 

Figure 22. D 2 Brightness in the Bergstralh Slit. Same description as Figure 
21, except a - 7/3. 

Figure 23. East-West Intensity Ratio. The east-west intensity ratio is shown 
in (a), as determined from the model results (a » 3) in Figure 21, 
and in (b) , as determined from the model results (a - 7/3) in 
Figure 22. The horizontal dashed line indicates the observed 
ratio. The vertical dashed line indicates the most probable speed 
where the average value between the 140° and 320° curves is equal 
to the observed average ratio. 
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Figure 24. D 2 Brightness in the Bergs tralh Slit. The model calculated 

brightness in the observing slit expressed in normalized flux 
units for Io at eastern (east) and western (west) elongation is 
shown as a function of the most probable speed of the Maxwell - 
Boltzmann flux distribution. The band source and the plasma torus 
description (Case C) assumed in the calculation are described in 
the text. The System III magnetic longitude angle of Io is 140° 
in (a) and 320° in (b) . 

Figure 25. East-West Intensity Ratio. The east-west intensity ratio as 

determined from the model results in Figures 24 for a Maxwell - 
Boltzmann flux distribution is shown as a function of the most 
probable speed (lower scale) and corresponding exobase temperature 
(upper scale) . The range of exobase temperatures from 1000°K to 
2000°K is indicated by the shaded area. The dashed line indicates 
the observed average ratio. 
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